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This  report  was  prepared  by  the  Industrial  and 
Experiment  Station,  Department  of  Mechanical  Engineer¬ 
ing,  University  of  Florida,  Gainesville,  Florida,  under 
Contract  No.  F08635-71-C-0073  with  the  Air  Force  Arma¬ 
ment  Laboratory,  Eglin  Air  Force  Base,  Florida,  during 
the  period  from  9  December  1970  to  9  December  1971. 
Lieutenant  Robert  J.  Karner  (DLWG)  monitored  the  pro¬ 
ject  for  the  Armament  Laboratory. 

The  principal  investigator  for  the  contractor  was 
Dr.  J.  Mahig. 

This  report  consists  of  two  volumes.  Volume  I  is 
devoted  to  the  Six-Degree-of-Freedom  Simulation  while 
Volume  II  is  concerned  with  the  Three-Degree-of-Freedom 
Simulation.  This  is  Volume  I. 


This  technical  report has  been  reviewed  and  is 
approved.  ^  / 


'HEYWVRD  H.  STRONo  [Y  / 

Actina  Chief,  Atr-to-Surlare  Guided  Weaoons  Div. 


ABSTRACT 


This  report  describes  a  six-degree-of-f reedom  pro¬ 
gram  which  can  be  used  to  determine  the  trajectory  and 
miss  distance  of  a  missile  system.  The  options  for  the 
program  are  such  as  to  permit  variation  of  the  aerody¬ 
namics,  seeker,  autopilot,  actuator,  and  missile  motor 
performance  for  the  purpose  of  accurately  simulating  a 
given  missile  design  and  evaluating  the  effects  of  any 
changes  in  system  parameters.  Sufficient  detail  has 
been  included  in  the  text  in  order  to  minimize  the  users’ 
effort  needed  to  know  how  to  update  or  modify  the  program 
for  his  purposes. 
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SECTION  I 


INTRODUCTION 

The  purpose  of  this  report  is  to  provide  a  ref¬ 
erence  which  will  enable  ready  access  to  the  use  of  a 
six-degree-of-f reedom  program  which  is  capable  of  ac¬ 
curately  determining  the  trajectory  and  miss  distance 
of  a  semi-active  or  passive  guided  missile.  The  pro¬ 
gram  is  divided  into  convenient  blocks,  called  modules 
or  subroutines,  which  do  specific  tasks:  e.g.,  deter¬ 
mine  aerodynamic  forces,  seeker  output,  state  of  auto¬ 
pilot,  current  value  of  thrust,  etc.  As  a  result,  the 
user  will  be  able  to  easily  locate  the  section  of  the 
program  where  specific  calculations  are  performed  and 
modify  them  or,  if  necessary,  to  add  modules  to  achieve 
other  purposes. 

This  program  has  been  derived  from  a  program  in 
the  library  of  the  North  American  Rockwell  Company, 
Columbus  Division,  and  is  described  in  NR  70H-232-1 
and  -2.  The  purpose  of  the  program  was  to  determine 
trajectory  and  miss  distance  of  an  air-to-air  or  air- 
to-surface  missile.  This  manual  goes  into  somewhat 
greater  detail  in  identifying  the  variables  and  defin¬ 
ing  coordinate  svstems  than  heretofore.  This  has  been 
possible  because  of  the  extensive  work  carried  out  with 
the  program  by  the  author  in  satisfying  the  requirements 
of  this  contract  and  information  supplied  by  Mr.  R.  D. 
Ehrich  and  P.  D.  Capcara  of  the  North  American  Rockwell 
(NAR)  Corporation.  The  program  described  below  has  been 
modified  from  the  original  version  supplied  to  USAF  by 
North  American  Rockwell  Corportion  by  personnel  at  the 
Air  Force  Armament  Laboratory  to  permit  the  considera¬ 
tion  of  the  effect  of  a  random  spot  motion  on  the  miss 
distance  of  a  laser  guided  missile.  Incorporated  into 
the  version  presented  in  this  report  are  additional 
capabilities  which  provide  an  accurate  simulation  of 
the  quadrant  detector,  range  closure,  proportional  lead 
guidance,  simplified  program  reset  mechanism  for  multi¬ 
ple  runs,  greater  target  maneuverability  in  air-to-air 
simulations,  and  a  more  general  high  frequency  actuator 
routine  which  will  accept  either  experimental  or  theo¬ 
retically  derived  transfer  functions. 
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SECTION  II 


PROGRAM  DESCRIPTION 
2 . 1  Subroutines ,  Modules,  and  Tables 

A  complete  listing  of  this  program  appears  in 
Section  V,  The  program  consists  of  three  types  of 
subprograms : 

(a)  Tables  of  aerodynamic  coefficients  in  block 

data  form. 

(b)  Modules  describing  missile  subsystems. 

(c)  Executive  subroutines  and  the  main  program. 

The  block  data  subroutines  must  be  physically  loca¬ 
ted  at  the  front  of  the  program  deck  after  the  main  pro¬ 
gram  for  proper  operation.  Data  is  extracted  from  these 
tables  in  the  module  Al  which  makes  use  of  the  table 
look-up  subroutines  TABLl ,  TABL2,  and  TABL3  which  form 
a  part  of  the  executive  routines. 

For  each  module  (e.g.,  Al ,  C4)  the  programmer  has 
the  option  of  using  an  associated  initialization  module 
(e.g.,  C4I) .  These  initialization  modules  may  be  used 
to  compute  initial  conditions  from  input  data  or  add  to 
the  list  of  state  variables  to  be  integrated.  The  ini¬ 
tialization  modules  are  executed  only  once  at  the  start 
of  each  simulated  mission.  It  is  in  the  modules  them¬ 
selves  (e.g.,  C4)  that  the  derivatives  of  the  state 
variables  are  computed.  Time  is  incremented  by  a  fixed 
amount  (At)  after  every  other  pass  since  a  predictor- 
corrector  integration  algorithm  is  used. 

A  large  block  common  array,  called  C,  allows  the 
communication  of  certain  variables  between  modules  and 
subroutines  for  input/output,  integration,  and  control 
purposes . 

The  mathematical  relationship  of  various  modules 
and  subroutines  are  shown  in  Figure  1,  and  a  corres¬ 
ponding  list  of  the  modules  is  given  in  Table  I. 
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Figure  1.  Flow  Chart  for  State  Variable  Calculations 
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TABLE  I.  SUBROUTINE  AND  MODULE  LIST  SIX-DEGREE- 
OF-FREEDOM  DIGITAL  PROGRAM 


I  GEOPHYSICAL  AND  EXTERNAL  ENVIRONMENT 

I 

G2  -  Steady  winds 

G3  -  Air  data  -  including  dynamic  pressure, 
density,  speed  of  sound 

G4  -  Terminal  geometry  -  computes  miss  distance 

G5  -  Transformations  of  position  and  velocity 
between  various  coordinate  systems 

II  SENSORS 

CIO  -  Spot  motion  -  including  boresight  error, 
aiming  error,  hotspot  motion,  etc. 

£1  -  Seeker  -  Seeker  performance  and  platform 

motion 

QUADET  -  yuadrant  detector  simulator 

III  COMPUTERS 

Cl  -  Autopilot  -  computes  steering  commands 
from  seeker  output 

C4  -  Actuators  -  includes  flap  motion  and 
limits 

IV  AIRFRAME 

A1  -  Aerodynamics  coefficients  -  table  look-up 

A2  -  Aerodynamic  forces  and  moments  -  in  wind 

axis,  includes  forces  and  moments  on. lugs 
while  missile  is  on  rail 

A3  -  Engine  -  computes  thrust  forces  as  well  as 
c.g.  shifts  and  mass  changes. 

V  DYNAMICS 

Dl  -  Translation  dynamics  of  missile  -  accelera¬ 
tions  in  body  axes  are  transformed  into 
earth  coordinates  and  integrated  into 
velocities  and  positions. 

D2  -  Rotational  dynamics  of  missile  -  computes 
rotational  accelerations  and  velocities 
referred  to  missile  body  axes. 
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2.2 


A2  -  Aero  Forces  and  Moments 


Figure  2  shows  the  relationship  between  the  body  axis 
and  wind  axis  coordinate  system.  In  addition,  the  coor¬ 
dinate  directions  are  shown  for  the  positive  direction  of 
the  dimensionless  aerodynamics'  coefficients  in  both  the 
body  axis  system  and  the  wind  or  primed  axis  system.  The 
body  axis  system  and  the  wind  axis  system  are  related  by 
the  following  system  of  equations: 


XB 

■ 

X*B 

II 

4>' 

y’b 

•  ZB 

. 

,z'b 

The  aerodynamic  coefficients  are  functions  of  the 
aerodynamic  roll  angle  { <J> 1 )  and  angle  of  attack  o'.  It 
can  be  seen  that  the  angles  <f>‘  and  o'  locate  the  wind 
vector  in  much  the  same  way  that  a  magnitude  r  and  angle 
6  locate  a  vector  in  polar  coordinates.  With  reference 
to  Figure  2,  it  is  apparent  that  the  plane  containing  the 
wind  vector  is  obtained  by  rotating  the  XBZB_  plane  through 
4> *  about  the  missile  centerline  (XB  axis).  The  wind  vec¬ 
tor  is  located  in  this  plane  by  the  angle  o'  measured  from 
the  XB  axis.  The  angles  <p '  and  a'  are  related  to  the  angle 
of  attack  a  and  sideslip  8  by  the  following  equations: 

Cos  a*  =  Cosa  Cos8 

Sin  4>  *  =  Sin  8/(l-Cos2a  Cos28)^ 

where  if  a  and  8  are  small,  one  finds 

a ' 2  =  a2  +  82  . 

Since  if  a  and  8  arc  small,  a'  will  similarly  be 
small  and  it  will  be  found  that 

6  =  a'Sin^' 

a  =  a’Cos$' 

tan$'  =  0/a 

a'  =  Va2+  62.' 
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♦'a  is  the  aerodynamic  roll  angle  referenced  to  2ero 
with  the  missile  flying  in  the  +  configuration.  If  the 
missile  is  intended  to  fly  in  the  X  configuration,  $'A 
equals  45°  with  8  =  0.  Thus,  4>'=4>'A-450.  M 

It  will  be  found  that  the  following  relationships 
hold  with  respect  to  $  and  <t>*A: 

Cos  4<j>'A  =  -Cos  4$' 

Sin  4<j>'A  =  -Sin  4$' 

Cos  ($A  -  45°)  =  Cos 

Sin  ($A  -  45°)  =  Sin 

Some  of  the  above  relations  can  be  experienced  in 
terms  of  the  angle 

In  order  to  facilitate  application  to  the  program. 
Table  II  lists  the  correspondence  between  variable  names, 
commonly  used  aero  symbols,  and  their  COMMON  location  in 
the  program. 
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TABLE  II.  CORRESPONDENCE  BETWEEN  VARIABLE  NAMES, 
AERO  SYMBOLS,  AND  THEIR  COMMON  LOCATION 


Name 

Symbol 

Common 

Location 

CLDRP 

Cn’/6r 

1225 

CNQ 

V/Sq 

1226 

CLD 

Ci'/fip 

1227 

CLMP 

c  * 

m 

1228 

CLNP 

<V 

1229 

BDEFL 

l«l 

1230 

CDCM 

<V  W) 

1231 

DDL 

6P 

1232 

BDM 

fiq 

1233 

BDN 

6r 

1234 

CDCN 

CN '  (4,*) 

1235 

CL2 

C*'  (0) 

1240 

CL3 

C£,(b)  lug 

1241 

CNPU 

CN '  ($'  ,a') 

1244 

CYPU 

CY '  (<M 

1245 

CMP 

Cm’  (a '  ,  4>  * ) 

1247 

CNP 

Cn '  (($>*) 

1248 

CLR 

cz'  (o') 

1249 

CZP 

CN ' 

1250 

CYP 

cY. 

1251 

2.3.1  D1 ,D2  -  Translational  and  Rotational  Dynamics  Module 


The  following  list  of  symbols  applies  to  the  equa¬ 
tions  of  motion  which  are  developed  in  paragraph  2.3.3 
for  modules  D1  and  D2. 


2.3.2  List  of  Symbols 


m 

P 

q 


r 


T 

^Z 


u 


V  *  ua 
W  =  Ug 

XB' YB' ZB 

Xe,Te , Ze 
a 
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F (t)  instantaneous  mass  (slugs) 

Rolling  velocity  =  angular  velocity  along 
X  axis  (rad/sec) 

Pitching  velocity  =  angular  velocity  along 
Y  axis  (rad/sec) 

Yawing  velocity  =  angular  velocity  along 
Z  axis  (rad/sec) 

Moment  of  inertia  about  X  axis  (slug-ft) 

Moment  of  inertia  about  Y  axis  (slug-ft) 

Moment  of  inertia  about  Z  axis  (slug-ft) 

(I2  =  Iy  for  a  perfectly  symmetrical  missile) 

Linear  velocity  along  the  X  body  (xE)  axis 
(ft/sec) 

Linear  velocity  along  the  Y  body  (Yg)  axis 
(ft/sec) 

Linear  velocity  along  the  Z  body  (Zg)  axis 
(ft/sec) 

Airframe  axis  system  that  moves  with 
airframes 

Earth  coordinates 

Angle  of  attack  =  angle  between  a  fuselage 
reference  line  and  the  relative  wind  in  the 
X3,ZB  plane  (rad) 

Tan  a  =  W/U ; a  =  W/U 

Angle  of  sideslip  (rad) 


Tan  &  =  V/  U2  +  W2 ;  S  *  V/U 
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Euler  angles 
<P,d ,  and  <p 


ip  is  the  rotation  about  ZB,  6  is  the 
rotation  about  YB,  and  <t>  is  the  rota¬ 
tion  about  XB  in  that  order  (rad) 

g  -  Acceleration  due  to  gravity  (ft/sec2) 

T  -  Thrust  along  XB 

cN,cy,cc,cR  Dimensionless  aerodynamics  coefficients 

Cr,/Cn  (body  axes) 

C'niC'yiC'i;  -  Dimensionless  aerodynamics  coefficients 
C'RfC'ntfC'n  (primed  axes  system  -  Figure  2) 

Density  (slug/ft3) 

q0  -  Dynamic  pressure  =  ^pu2 (lb/ft2) 

5  -  Body  reference  cross  sectional  area  (ft2) 

l  -  Reference  body  length  (ft) 

AX  -  Shift  of  center  of  gravity  from  a  refer¬ 
ence  point  along  the  XB  axis  (ft)  - 
negative  aft 

6  -  Control  surface  deflection  (rad) 

6  -  Control  surface  deflection  to  give  pitch¬ 

ing  motion  (rad) 

6p  -  Control  surface  deflection  to  give  rolling 
motion  (rad) 

<5r  -  Control  surface  deflection  to  give  yawing 

motion  (rad) 

C^j  -  Dimensionless  aerodynamic  derivatives 

d'  -  Aerodynamic  or  wind  angles  of  attack 
(rad)  -  Figure  2 

4>'  -  Aerodynamic  roll  angle  (rad)  -  Figure  2 

*  a  "  Aerodynamic  roll  angle  (rad)  referenced  to 
zero  when  flying  in  the  +  configuration. 

<t>  •  A  =  4>  *  +  45° 
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2.3.3  Equations  of  Motion 


The  six-degree-of-f reedom  equations  of  motion  im¬ 
plemented  in  the  computer  program  in  terms  of  the  body 
axes  are  given  below*.  (See  Figure  3  for  coordinate 
system  orientation.) 

Longitudinal  Force 

=  m[U  +  Uaq  -  U^r]  =  q0SltCc  -  mgSinQ  +  T 
Lateral  Force 

ZFy  =  m (d/dt  (U^ )  +  Ur  -  Uap]  =  qoS^  [Cost)  'C  ’  Y  - 
Sin4>*C'j.,]  +  mgCosBSintj) 

Vertical  Force 

EF  =  m[d/dt(Ua)  +  Uep  -  Uq)  =  -q0 [Cos$ ' C ’ N  + 
Simj'C'y]  +  mgCos0Cos£ 

Rolling  Moment 

EMx  =  !xP  =  q«SwUC4  +  £/2UCipP) 

Pitching  Moment 

ZMy  =  IYq  +  (Ix  -  I2)pr  -  q0S^[Cos*'C’m  +  Sin*'C'n  + 
£/2U  Cmqq  -  AX/MSin^'Cy  +  Cos<t>*C'N)] 

Yawing  Moment 

IMZ  =  Izr  +  (IY  -  Ix)pq  =  qvSTtS-lCos<J',C'n  -  Sin^'C^  + 
£/2U  C_  r  -  AX/£  (Cos^’C’v  -  Sin<J>’C\,)] 

Uj1  1  »■* 

*  UMj  velocity  in  X  direction 

VMJ  velocity  in  Y  direction 

WMJg  velocity  in  Z  direction 
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Missile  c.g. 


Figure  3.  Definitions  of  Angles  and  Coordinate  Systems 


The  Euler  transformations  of 


'C  =  p  =  rCos<l>/Cos0  +  qSin<j>/Cos0 

6  =  q  =  qCos$>  -  rSm4> 

<p  =  r  =  p  +  (rCos$  +  qSiniJ>)  tan© 

The  velocity,  in  terms  of  the  earth  axes,  can  be  obtained 
as : 


X  =  UCos6CosT  +  Ug  (Sin$Sin6CostJj  -  Cos<f>Sin^) 

+  Ua  (Cos0Sin0Cosip  +  Sin^Sinj)) 

Y  =  UCosQSiniJi  +  Ug  (Cas4>Cos^  +  Sin^Sin9Sin<{) ) 

+  Ua  (Cos<f>Sin8SiniJj  -  Sin<J>CosiM 

2  =  -  USin0  +  UgSin4>Cos0  +  UaCosQCos<t> 

Tne  Euler  angles,  shown  in  Figure  4,  and  the  posi¬ 
tion  of  the  missile  in  earth  coordinates  can  be  obtained 
through  an  integration  of  the  above  equations. 

The  block  diagram  of  the  implementation  of  the  equa¬ 
tions  of  motion  and  the  Euler  transformations  are  shown  in 
Figure  5. 

2.4  Subroutine  G2 


This  subroutine  is  called  the  wind  and  gust  module. 

This  module  determines  the  velocity  and  direction  of  the 
wind.  The  module  assumes  that  there  is  no  wind  above  an 
altitude  RHW.  Below  that  altitude  the  wind  direction  and 
magnitude  are  assumed  to  be  constant  throughout  a  layer 
RWINC  in  depth.  (It  should  be  noted  that  RWINC  is  measured 
along  the  line  of  sight.  Since  most  missiles  fly  with  only 
small  deviation  from  the  original  line  of  sight,  the  alti¬ 
tude  increments,  if  needed,  may  be  readily  estimated.)  Two 
random  variables  are  associated  with  the  wind  in  each  layer: 
the  magnitude  and  angular  orientation  which  are  considered 
constant  in  each  layer.  The  mean  value  of  the  wind  magni¬ 
tude  is  VWTE,  and  its  standard  deviation  is  given  as  SW.  The 
mean  value  of  the  angular  orientation  of  the  wind  in  a  layer 
is  3PSIW,  and  the  standard  deviation  of  the  angular  variation 
is  SWl .  The  current  value  of  the  wind  magnitude  and  direction 
is  given  by  VWTEV  and  BPS1WV,  respectively.  The  relationship 
between  these  mean  values  and  the  inertial  coordinate  system 
is  shown  in  Figure  6, 
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(Earth  Acad 


FIGURE  5.  SIX-DEGREE-OF-FREEDOM  EQUATIONS  OF  MOTION 


2.5  Subroutine  CIO 


This  subroutine  determines  the  ground  plane  coordi¬ 
nates  of  that  point  in  the  area  illuminated  by  the  laser 
beam  which  the  missile  seeker  physics  causes  the  autopilot 
to  consider  the  designated  target.  This  distinction  is 
necessary  since  some  seekers  are  hotspot  trackers  while 
others  are  cer.troidal  trackers.  The  procedure  used  to  de¬ 
velop  this  apparent  target  is  accomplished  first  through 
the  designation  of  the  coordinates  of  the  illuminator 
(XIL,  HILL)  which  may  either  be  on  the  forward  air  con¬ 
troller  or  on  the  launch  aircraft.  The  maximum  errors 
generated  on  the  ground  are  considered  to  be  made  up  of 
three  parts:  the  maximum  boresight  error  (BORE) ,  the  maxi¬ 
mum  pointing  error  (WAND) ,  and  the  maximum  deviation  of  the 
hotspot  from  the  resulting  beam  centroid,  which  is  desig¬ 
nated  as  (RADIUS) .  Each  of  these  variables  is  considered 
a  random  variable  with  a  uniform  distribution.  The  re¬ 
sulting  random  variables  generated  are,  respectively, 

BOREF,  WANDF,  and  SPWID.-,  The  variables  are  considered  to 
vary  independently  in  the  XE  and  YE  direction  and  are  then 
appropriately  summed  in  order  to  determine  the  apparent 
target  location.  The  coordinates  of  this  point  are  desig¬ 
nated  as  the  variables  ZLASR  and  YLASR.  The  location  in 
earth  coordinates  may  be  found  by  equating  ZLASR  to  XE  and 
YLASR  to  YE  and  setting  ZE  equal  to  zero. 

2.6  Subroutine  QUADET 


Subroutine  QUADET  is  called  by  Si  for  the  determina¬ 
tion  of  the  signal  generated  to  the  autopilot  by  the  quad¬ 
rant  detector  (Figure  7) .  The  quadrant  detector  is  orient¬ 
ed  such  that  the  dead  zone  is  in  the  same  direction  as  the 
fins,  assuming  the  missile  flies  in  the  X  configuration. 

The  subroutine  determines  the  current  size  of  the  circular 
image  through  the  assumption  that  the  image  size  increases 
inversely  proportional  to  the  range  of  the  missile  from  the 
laser  spot.  RTl  is  the  variable  designating  the  ratio  of 
the  size  of  the  current  spot  to  its  size  at  infinity.  The 
laser  image  on  the  detector  is  assumed  to  be  circular.  In 
order  to  determine  the  portion  of  the  area  of  each  quadrant 
covered  by  the  laser  image,  the  area  of  the  detector  is  sub¬ 
divided  into  LT  segments.  (In  the  current  program  LT  is  set 
equal  to  16.)  In  order  to  effect  a  dead  zone,  an  area  round 
the  axis  of  the  coordinate  system  equal  to  half  the  segment 
width  is  not  included  in  the  area  of  the  image  which  cover 
these  segments.  If  a  portion  of  the  laser  image  falls  off 
the  assumed  circular  detector's  surface,  it  is  not  consid¬ 
ered.  The  variable  DETRID  is  half  the  instantaneous  field 
of  view  of  the  detector  in  degrees.  DEFICS  is  half  the 
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instantaneous  field  of  view  intercepted  by  the  image  of 
the  laser  spot  on  the  detector.  The  subroutine  will  de¬ 
termine  if  the  following  blind  range  and  breaklock  cri¬ 
teria  are  met  and  print  this  information  on  the  line 
printer.  The  breaklock  criteria  is  met  if  there  is  no 
portion  of  the  laser  image  on  any  of  the  four  quadrants. 
The  blind  range  criteria  is  met  if  the  image  of  the  laser 
spot  on  the  detector  exceeds  90  percent  of  the  total  area 
of  the  detector. 

2.7  Subroutine  SI  (Module) 


The  purpose  of  the  SI  Module  is  to  simulate  the  re¬ 
sponse  of  several  types  of  seeker  models  and  to  generate 
the  commai  is  which  are  transmitted  to  the  autopilot. 

The  subroutine  will  simulate  the  seeker  response  to 
either  a  continuous  information  source  or  a  sampled  data 
source.  This  is  accomplished  hy  setting  the  variable 
OPTKR  either  to  zero  or  one,  respectively.  If  operating 
from  a  continuous  information  source,  the  seeker  is  as¬ 
sumed  by  the  module  to  be  a  proportional  seeker.  In  the 
sampled  data  mode  the  seeker  can  be  programmed  as  either 
a  proportional  or  a  bang-bang  seeker  by  the  choice  of  the 
magnitude  of  the  variable  DEFOCS.  If  this  variable  is 
chosen  so  that  it  is  equivalent  to  DELF  [detector  radius/LT 
(in  current  program) ] ,  the  seeker  will  simulate  a  bang-bang 
laser  seeker;  whereas,  if  this  variable  is  chosen  so  that 
it  is  larger  than  DELF,  it  will  produce  a  proportional 
laser  processor. 

In  the  sample  data  mode  the  seeker  will  simulate 
either  a  pursuit  or  a  proportional  navigation  system  by 
setting  the  variable  CAGE  equal  to  zero  or  one,  respective¬ 
ly.  In  the  continuous  information  mode,  corresponding 
changes  in  the  guidance  law  will  occur.  In  either  mode 
of  operation  the  PLG  option  may  be  implemented.  This  is 
done  by  removing  the  C  from  the  two  cards  following  the 
PLG  OPTION  card. 

The  mode  of  operation  of  the  subroutine  in  either 
mode  is  to  initially  determine  the  true  location  of  the 
target  in  the  gimbal  axis  coordinate  system  (RXG,  RYG ,  and 
RZG)  and  then  determine  the  angles  the  lines  of  sight  make 
with  the  RYG,  RXG  plane  and  the  RZG,  RXG  planes  (BEPSY  and 
BEPSZ,  respectively,  shown  in  Figure  8  and  Figure  9). 
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TARGET 


Figure  8.  Coordinate  Relations  Between  Body 
and  Gimbal  Axis  System  and  the  Line  of  Sight 


Gyroscope  Platform  Gimbal  Angles 
0g  outer  gimbal  -  pitch 
inner  gimbal  -  yaw 

Figure  9.  Schematic  Diagram  of  Platform  Gimbal  Angles 
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The  rate  of  pulse  loss  is  determined  by  the  value, 
between  zero  and  one,  initially  given  the  variable  VLAZRP. 
This  is  done  by  comparing  a  uniformly  distributed  random 
variable  [C(103)]  whose  range  is  also  between  zero  and  one 
with  VLAZRP.  If  it  is  greater,  it  is  assumed  that  the  in¬ 
formation  in  the  pulse  is  lost.  If  pulse  loss  has  not 
occurred,  then  the  apparent  location  of  the  target  is  de¬ 
termined  in  the  gimbal  axis  coordinate  system  which  has 
resulted  from  boresight  errors,  wander,  etc.  Subroutine 
3UADET  is  then  called  to  determine  the  output  of  the  quad¬ 
rant  detector.  This  output  is  used  to  generate  the  re¬ 
quired  gimbal  rate  and  autopilot  commands.  If  pulse  loss 
has  occurred,  previously  generated  commands  (e.g.,  gimbal 
rate,  autopilot  signals)  are  maintained. 

In  addition,  Appendix  I  shows  the  mechanics  of  the 
coordinate  transformation  from  the  body  axis  to  the  gimbal 
axis  system  for  easy  reference. 

2 . 8  Cl  Autopilot  Module 

The  following  high  and  low  frequency  autopilot  block 
diagrams  are  suitable  representations  for  an  autopilot 
that  would  prove  to  be  consistent  with  either  a  propor¬ 
tional  or  bang-bang  seeker.  The  block  diagrams  for  each 
of  those  autopilots  are  given  in  Figures  10,  11,  12,  and 

13.  These  systems  correspond  to  those  mechanized  in  the 
program  listing  found  in  Section  V  for  the  low  frequency 
autopilot  and  in  Appendix  III  for  the  high  frequency  auto¬ 
pilot. 

2 . 9  C4  -  Actuator  Module 

The  actuator  module  simulates  the  action  ■  f  the 
actuator  up  to  a  third  order  system,  as  shown  in  Figure 

14,  which  corresponds  to  a  high  frequency  actuator.  Under 
these  conditions  it  is  capable  of  simulating  the  dynamics 
of  either  a  torque  balance  system  whose  block  diagram  is 
shown  in  Figure  15,  or  that  of  the  position  loop-controlled 
actuator  shown  in  Figure  16.  It  will  also  simulate  the 
dynamics  of  an  actuator  whose  transfer  function  has  been 
determined  from  hardware  test  data  up  to  the  third  order. 

The  transfer  function,  given  in  general  form  as  ex¬ 
pressed  in  this  module,  is  shown  below: 

£  =  K _ _ 

<5C  7il  *  S3  +  A2  *  Sz  +  A3  *  S  +  A4 
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FMH(I) 


FIGURE  14.  HIGH  FREQUENCY  ACTUATOR 


aerodynamic  load 


Actuator  Torque  balance  System 


figure  16.  Actuator  Position  Loop  Block  l) 


The  transfer  function  for  either  the  position  feed¬ 
back  system  or  the  torque  balance  system  can  be  brought 
into  the  following  form: 

6  _  KVA*  £ (57.3) 

Tc  Kl  [s  (KlJS  (TS+l)  +  Kr(tS+1)  +  (A  4)Vkl)  + 

K£K657. 3 (tS+1) ] 

and  similarly  for  the  position  loop  block  diagram. 

If  either  the  torque  balance  system  or  the  position 
loop  control  system  is  to  be  activated,  then  CKACT  should 
be  either  set  equal  to  one  or  zero,  depending  on  whether 
the  aerodynamic  tables  for  FMHl,  FMH2,  FMII3,  and  FMH4  are 
included  in  the  data  tables.  The  variable  BDMAX  limits 
the  maximum  amplitude  of  the  fin  motion.  The  low  fre¬ 
quency  actuator  equations  are  developed  below. 


Low  Frequency  Actuator 


BDELT(l)  = 
BDELT (2 )  = 
BDELT (3)  = 
BDELT (4)  = 
6,  = 

*3  = 

= 

where 


BDELT (1)  - 
BDELT  (2)  - 
BDELT (3)  + 
BDELT (4)  + 
BDELT  (1) 
BDELT (2) 
BDELT ( 3 ) 
BDELT (4 ) 

DELTPB 

DELTQB 

DELTRB 


6p 

+ 

6q  - 

6  r 

6p 

+ 

6q  + 

6r 

6p 

■f 

6q  - 

6r 

6p 

+ 

6q  + 

6r 

The  mechanization  of  these  equations  may  be  found 
in  the  Program  Listing  (Section  V)  for  the  low  frequency 
actuator.  The  high  frequency  actuator  program  listing 
may  be  found  in  Appendix  II. 
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2.9.1  Fin  Deflection 


A  positive  pitch  rate  (motion  up)  is  obtained  with 
a  negative  6q,  where 


A  positive  roll  rate  (motion  clockwise  about  the  X 
body  axis)  is  obtained  with  a  positive  6p,  where 


A  positive  yaw  rate  (motion  clockwise  about  the  2 
body  axis)  is  obtained  with  a  negative  6r. 

A  positive  surface  deflexion  is  defined  as  a  trail¬ 
ing  edge  down.  The  surfaces  are  labeled  by  looking  at  the 
missile  tail-on,  with  6 t  being  the  upper  right  surface,  62 
the  lower  right  surface,  6  3  the  lower  left  surface,  and 
the  upper  left  surface,  as  shown  in  Figure  l7. 

It  is  assumed  that  the  surface  effectiveness  will  be 
given  in  terms  of  6q,  6p,  and  6r  as  a  function  of  a'  and  4> '  . 
These  terms  will  be  considered  as  a  part  of  the  aerodynamic 
coefficients  given  in  the  primed  axis  system. 


2.10  A3  -  Engine  Module 


As  a  result  of  various  sources  of  error  occurring  in 
the  manufacture  and  assembly  of  a  solid  propellant  motor, 
the  thrust  alignment  is  not  perfect.  The  coordinate  system 
used  in  determining  the  misalignment  the  user  wishes  to 
simulate  is  shown  in  Figure  18. 


RFXCG  -  X  Component  of  thrust  vector  with  respect  to 
body  axis  in  the  X  direction 

RFYCG  -  Y  Component  of  thrust  vector  with  respect  to 
body  axis  in  the  Y  direction 

RFZCG  -  Z  Component  of  thrust  vector  with  respect  to 
body  axis  in  the  Z  direction 

FTHRST  -  Missile  Thrust 

Figure  18.  Offset  Thrust  Vector  Coordinate  System 


SECTION  III 


VARIABLE  LOCATIONS 


3, 1  Variable  Names,  Block  Locations,  and  Definitions 

Since  the  proper  use  of  this  program  requires  that 
the  definition  of  upwards  of  five  hundred  singly  dimen¬ 
sioned  variables  as  well  as  many  multidimensioned  be 
made,  it  is  clear  that  some  order  must  be  maintained  in 
the  allocation  of  storage  or  serious  programming  diffi¬ 
culties  could  arise.  Therefore,  blocks  of  common  loca¬ 
tion  have  been  allocated  to  specific  subroutines  as 
shown  in  Table  III.  This  procedure  should  be  continued 
in  the  event  it  is  necessary  to  add  variables  as  a  re¬ 
sult  of  program  modifications. 

Of  the  large  number  of  variables  actually  listed  by 
the  program,  only  two  hundred  and  fifty  appear  to  be 
significant  in  the  preparation  of  the  input  or  of  an 
aid  in  understanding  the  output.  Therefore,  it  was  felt 
that  they  should  be  separately  defined.  This  is  done  in 
Table  IV.  It  should  be  noted  that  the  units  of  the  vari¬ 
ables  listed  in  that  table  should  be  considered  to  be  in 
feet,  seconds,  pounds,  or  degrees  unless  otherwise  speci¬ 
fied. 


TABLE  III.  BLANK  COMMON  ASSIGNMENTS 


Array  Index 

C(  ) 

Module 

Name 

Module  Description 

0-50 

CIO 

Unsteady  Illuminator 

50  -  102 

G2 

Steady  Winds,  Variable  Winds 

200  -  299 

G3 

Air  Data 

350  -  399 

G5 

Coordinate  Conversion 

400  -  499 

SI,  SlI 

Seeker  -  Platform 

800  -  899 

Cl,  C1I 

Autopilot 

1100  -  1149 

C4 ,  C4I 

Actuators 

1200  -  1299 

Al 

Aero  Table  Look-Up 

1300  -  1399 

A2 

Forces  and  Moments 

1400  -  1499 

A3,  A3I 

Engine 

1600  -  1699 

Dl,  Dll 

Translational  Dynamics 

1700  -  1799 

D2,  D2I 

Rotational  Dynamics 

Note:  Locations  1950  -  4310  are  reserved  for  Executive 
Subroutines,  Initial  Conditions,  and  Input-Output  Instruc¬ 
tions. 
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TABLE  IV.  COMMON  LOCATION ,  VARIABLE  NAME,  AND  DEFINITION 


Common 

Location 

Variable 

Name 

Definition 

C  (1) 

BORE 

Maximum  boresight  error 

C  (2) 

WAND 

Maximum  pointing  error 

C  ( 3) 

RADIUS 

Maximum  deviation  of  hotspot 
from  beam  centroid 

C  (4 ) 

HILL 

Height  of  illuminator 

C(5) 

RILL 

Ground  range  of  illuminator 

C  (6) 

AISPOT 

{ 

0.  -  Centroid  tracker 

1 .  -  Hotspot  tracker 

C  (7) 

AILL 

1 

0.  -  Stationary  illuminator* 

1.  -  Moving  illuminator 

C  (8) 

SPOTMO 

0.  -  No  spot  motion 
[1,  -  Spot  motion 

C  (9) 

XSPOT 

X  -  Coordinate  of  centroid  or 
hotspot 

C  UO) 

YSPOT 

Y  -  Coordinate  of  centroid  or 
hotspot 

C(ll) 

A I  FAC 

0.  -  Tracker  on  ground  or  on 
launch  aircraft 

1.  -  Tracker  on  separate  air- 
>*  craft 

CU2) 

VI LM 

Velocity  of  illuminator,  Mach 
number 

C  (18) 

XILL 

Ground  range  in  X  direction  of 
illuminator 

^Must  give  HILL,  XILL  for  input 
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TABLE  IV  (Continued) 


Common 

Location 

Variable 

Name 

Definition 

C  { 51) 

BPSIW 

Mean  angle  of  wind  velocity 
vector 

C  (52) 

VWTE 

Mean  wind  velocity 

C  (53) 

RHW 

Altitude  above  which  all  the 
winds  are  zero 

C  (55) 

sw 

Standard  deviation  from  mean 
angle  of  wind  velocity  vector 
BPSIW 

|C (56) 

1 

RWINC 

Shear  layer  of  wind.  Depth  of 
wind  layers  at  which  wind  velo¬ 
city  and  angle  remain  constant. 

C  (56) 

SW1 

Standard  deviation  from  mean  wind 
velocity  VKTE 

C(100) 

VWXE 

Wind  velocity  (X  component  with 
reference  to  the  earth-fixed 
coordinate  system) 

C(101) 

VWYE 

Wind  velocity  (Y  component  with 
reference  to  the  earth-fixed 
coordinate  system) 

C  (102) 

VWZE 

Wind  velocity  (Z  component  with 
reference  to  the  earth-fixed 
coordinate  system) 

C  (203) 

PDYNMC 

Dynamic  pressure 

C  (204) 

VMACH 

Mach  number 

C(205) 

DRHO 

Air  density 

C  (206) 

VSOUND 

Speed  of  sound 

C  (207) 

VAIPSP 

Missile  velocity  v/ith  respect  to 
air  mass  in  earth  axes 
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1ABLE  IV  (Continued) 


Common 

Location 

Variable 

Name 

De  f inition 

C (208) 

RhZRO 

Initial  altitude  of  the  missile 

C (200) 

R1I 

Present  altitude  of  the  missile 

2  v  3  5  0  } 

BTHT 

Pitch  angle  in  degrees  -  G 

2 \ 35 1 ) 

LPSI 

Roll  angle  in  degrees  -  if 

C (332) 

bPHI 

Yaw  angle  in  degrees  -  <f 

C ( 356) 

VTOTE 

Tc  il  missile  velocity 

C (357) 

BGAMH 

Horizontal  proportional  naviga¬ 
tion  angle  (degrees) 

C (35S) 

LGAMV 

Vertical  proportional  naviga¬ 
tion  angle  (degrees) 

C (367) 

i 

IS  ALPHA 

Vertical  component  of  angle  of 
attack 

'iC  (  368) 

BALPHY 

Horizontal  component  of  angle 
of  attack 

:  C  (  3  C  9 ) 

1 

1 

1 

BALPHP 

(ci'  =  VbALPHA2  +  BALPHY2) 
total  angle  of  attack 

!  2  (  3 70 ) 

bPilIF 

v’  orientation  of  wind  vector 
in  roll  axis 

|c (371) 

RANGE 

Range 

I  2  (  372) 

11XBA 

Range  (X  component  in  body  co¬ 
ordinate  system) 

2  ,3*3) 

ki'BA 

Range  (Y  component  in  body  co¬ 
ordinate  system) 

1 

i  2 ( 37  4) 

K2BA 

Range  (2  component  in  body  co¬ 
ordinate  system) 
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TABLE  IV  (Continued) 


Common 

Lee  ci  t  ion 

Variable 

Name 

Definition 

C ( 377 ) 

BALPD 

d (BALPHA) /dt 

C (378) 

BALYD 

d ( BALPHY ) / dt 

C(379) 

BALPPD 

d (BALPHP) /dt 

C ( 380) 

RAN  GO 

The  distance  of  missile  from  the 
laun_h  point 

C (403) 

E2 

Seeker  output  to  autopilot  (pitch) 

C (4  07) 

EY 

Seeker  output  to  autopilot  (yaw) 

C  (427 ) 

1 

13THTG 

Platfc  -m  position  (9g) 

C ( 4  31 ) 

BPSIG 

Platform  position  yaw  gimbal 
angle  (4*g) 

C  (4  32 ) 

RXG 

Range  X  in  gimbal  axes 

C  (433) 

RYG 

Range  Y  in  gimbal  axes 

C (4  34) 

RZG 

Ranee  Z  in  gimbal  axes 

!  C  (435) 

BEPSZ 

Angular  position  of  the  line  of 
sight  in  gimbal  axes  (see 

C  (435) 

BEVSY  \ 

Figure  8) 

C  ( 4  3  7 1 

WZ 

j  Missive  body  rate 

C  (4  38 ) 

1 

l  wy 

i 

iirssile  body  rate  Wy 

C  (441) 

i  SZGBtS 

Pitch  girbal  torque  bias  (deg/sec) 

jc  (442) 

1 

1  SYCBIS 

j  Yaw  gimbal  torque  bias  (deg/sec) 

j  C  ( 4  4  3  ) 

{  OPT NR 

>  Optics  routine 

1 

C  (4  44 ) 

1 

OI’TBKL 

Optical  breaklock 

C  ( 44  5 ) 

UT 

Time  at  which  next  pulse  expected 

C (44  6 ) 

CDT 

Pulse  rate,  sampling  period 

TABLE  IV  (Continued) 


Common 

Location 

Variable 

Name 

Definition 

C  (450) 

CSK 

Acquisition  gain  seeker  constant 

j  C  (  4  5 i) 

I  I 

CAGE 

[(>0)  uncaged  gimbals]  [ ( <  0 ) 
remain  in  caged  position] 

’!  c  (  4  5  2 ) 

1 

Q  BREAK 

(Breakiock  has  occurred  due  to 
loss  of  signal) (automatically 
parameterized) 

C  ( 453) 

RBKLOK 

Range  at  breakiock  (maximum  range 

at  which  lock-on  can  take  place) 

jC(454) 

LEDGE 

Half  the  field  of  view 

!  C  ( 4  5  5) 

WEPSMX 

Breakiock  drift  rate 

|  C  (  4 56 ) 
j 

CKEKR 

Seeker  gain 

jc  (457) 

1 

Cl'OSTP 

Pitch  to  yaw  friction  coupling 

o 

.fcw 

in 

CO 

CROSPT 

Vav  to  pitch  friction  coupling 

j  C (4  6  0) 

CHIDE 

(=1  missile  guidance  system  in 
effect)  (=0  missile  guidance 
system  not  in  effect) 

jc  (4  61) 

i 

i 

i 

i 

1 

I 

1 

SAMP 

Preprogrammed  guidance  trajectory 
(cutoff  check  automatically  para¬ 
meterized)  (0  -  missile  uses  pre¬ 
programmed  flight  path)  (1  -  mis¬ 
sile  uses  preprogrammed  flight 
path  until  seeker  acquires  target) 

C  \  4  6  4  ) 

CGAMVS 

Vertical  trajectory  programming 
constant 

j c  ^4  o  5  ) 

i 

CGA.'lilS 

Horizontal  trajectory  program], ling 
constant 

C  (4  06) 

ZLAZR 

Location  of  laser  spot  on  target  in 
X  direction  due  to  ground  or  air¬ 
borne  FAC 

* 
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TABLE  IV  (Continued) 


Common 

Location 

Variable 

Name 

Definition 

C  (467) 

YLAZR 

Location  of  laser  spot  on  target 
in  Y  direction  due  to  ground  or 
airborne  FAC 

C (468) 

DEFOCS 

Half  angle  in  degrees  of  angle 
intercepted  by  image  of  laser 
spot  on  detector  surface 

C  (469) 

DETRAD 

Half  angle  in  degrees  of  angle 
intercepted  by  quadrant  detector 

C  (472) 

CKSKl 

Seeker  driver  constant 

C  (473) 

VLAZRP 

Used  in  pulse  loss  calculation 

C  (850) 

HLIMO 

Limit  on  6C  from  pitch  and  yaw 
plane  (deg)  (fins  1  and  3) 

C  ( 851) 

HLIME 

Limit  on  6C  from  pitch  and  yaw 
plane  (deg)  (fins  2  and  4) 

C (852) 

i 

| 

QBIAS 

Pitch  body  rate  bias  (deg/sec) 
(used  as  "g"  bias) 

i C (853) 

RBI  AS 

Yaw  body  rate  bias  (deg/sec) 

C  (855) 

GZ 

Navigation  ratio  for  pitch  plane 

C (856) 

GY 

Navigation  ratio  for  yaw  plane 

C  (863) 

TAUZ 

Pitch  guidance  lag  filter 
(rad/sec) 

C (864 ) 

TAUY 

Yaw  guidance  lag  filter  (rad/sec) 

C  (865) 

TDY1 

Rate  loop  gain  switch  1  (sec) 

C ( 866 ) 

TDY2 

Rate  loop  gain  switch  2  (sec) 

C (  87  7 ) 

TAUL 

Guidance  lead  filter  (rad/sec) 

C  (888) 

CKSK2 

Seeker  gain  constant 

TABLE  IV  (Continued) 


Common 

Location 

Variable 

Name 

Definition 

CU117) 

BSURF  1 

*,] 

C (11 18 ) 

BSURF  2 

fin  deflection  output 
for  fins  1,  2,  2,  4 

C  (1119) 

BSURF  3 

C  (1120) 

BSURF  4 

J 

,  C  (1160) 

l 

DELTPB 

=  (-  6,-6a  +  6,  +  «„)/4 

C  (1161) 

DELTQB 

= 

=  (6  j  +  6  2  +  o3  +  5v)/4 

C  (1162) 

DELTRB 

6r  = 

=  (-5j  -6  2  +  6  j  +  6  4)  /4 

C  (1260 ) 

CXERR 

Drag  coefficient  error 

C (1261) 

CZERR 

Normal  force  (C  *  2 )  coefficient 
error 

C (1262) 

CYERR 

Side  force  (C*  )  coefficient 
error 

C  (1263) 

CLERR 

Ro3 l  moment  (C'L)  coefficient 
error 

j  C (126  4) 

CMERR 

Pitch  moment  (C'M)  coefficient 

error 

2  ( i  2  6  5 ) 

CNERR 

Yaw  moment  ( C  ' -^ )  coefficient 
error  1 

C ( 1 300) 

FXBA 

The  X  component  of  aero  force  in 
body  coordinate  system 

C  (1301) 

FYBA 

The  Y  component  of  aero  force  ir, 
body  coordinate  system 

C ( 1302 ) 

FZBA 

The  Z  component  of  aero  force  in 
body  coordinate  system 

C  (1  303) 

FMXBA 

The  X  component  of  aero  moment  in 
body  coordinate  syste 

C ( 1 304) 

I- MY  BA 

The 

i  component  of  aero  moment  in 

1 

body  coordinate  system  | 
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TABLE  IV  (Continued) 


Definition 


The  2  component  of  aero  moment 
in  body  coordinate  system 

Missile  reference  area  (ft2) 

Missile  reference  length  (ft) 

Centsr  of  gravity  shift  (ft) 


Hinge  moments 


RFXCG 


RFYCG 


FMXTH 


FMYTII 


FMZTH 


FMXLUG 


FMYLUG 


FMZLUG 


Thrust  vector  displacements  (ft) 


Distance  between  lugs  (ft) 

Rail  length  (ft)  (between  rear 
of  front  lug  and  end  of  rail) 

X  component  of  moment  caused  by 
thrust  misalignments 

Y  component  of  moment  caused  by 
thrust  misalignments 

Z  component  of  moment  caused  by 
thrust  misalignments 

X  component  of  moment  due  to  lugs 

Y  component  of  moment  due  to  lugs 

2  component  of  moment  due  to  lugs 


TABLE  IV  (Continued) 


Common 

Location 

Variable 

Name 

Definition 

C (1401) 

C (1402) 

DALPHT  ' 

LPHIT 

The  angles  as  indicated  in 

Figure  18 

C (1403) 

QNALGN 

(>  0;  include  thrust  misalignment 
angles) 

C (1404) 

FCFTH 

Fractional  increase  in  total 
thrust 

C (1405) 

QbUR'-l 

Parameterized  by  program 

C (1410) 

FTHRST 

Missile  thrust 

C (1411) 

FT  MX 

X  component  of  missile  thrust 

C (1412) 

FTHV 

Y  component  of  missile  thrust 

C  (1413) 

FT  HZ 

Z  component  of  missile  thrust 

C  ( 14 1 4  ) 

CISP 

Specific  impulse  (lb/sec) 

C  (1415) 

DWT 

Total  missile  plus  propellent 
wt  (lb)  initial 

C  (1416) 

DWP 

Propellent  weight  (lb) 

C  (1417) 

PJDCGO 

Initial  value  of  c.g.  shift  (ft) 

C  (1418) 

HDCGF 

Burnout  value  of  c.g.  shift  (ft) 

C  (1413) 

FMIXO 

Initial  value  of  moment  of  in¬ 
ertia  about  the  roll  axis 
(slugs  ft) 

C (1420) 

FMIYO 

Initial  value  of  moment  of  in¬ 
ertia  about  the  pitch  axis 
(slugs  ft) 

C  (1421  ) 

FLCGA 

Distance  between  launch  c.g.  and 
rear  lug  (ft) 

C  (1422) 

RLCG 

Present  position  of  c.g.  of  missile 
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TABLE  IV  (Continued) 


Common 

Variable 

Location 

Name 

Definition 

C (1496) 

UIMPD 

Thrust  of  the  motor 

C (1499) 

UIMP 

Impulse  of  the  motor 

C (1603 ) 

VXE 

X,Y,Z  coordinates  of  missile 

C (1607) 

VYE 

velocity  with  respect  to  the 
earth  fixed  coordinate  system 

C (1611) 

VZE  J 

C (1615) 

RXE 

X,Y,Z  coordinates  of  missile 

C (1619) 

RYE 

| 

c.g.  with  respect  to  the  earth 
fixed  coordinate  system 

C  ( 162 3) 

RZE 

1 

C (1624) 

AX  BA 

X  component  of  acceleration  in 

body  coordinate  axis 

C (1625) 

AYbA 

Y  component  of  acceleration  in 

body  coordinate  axis 

C ( 1626) 

A  2  BA 

Z  component  of  acceleration  in 

body  coordinate  axis 

C  (1627) 

AGRAV 

Gravitational  constant 

C ( 1628) 

DtlASS 

Current  mass 

C  (1629 ) 

ATHRST 

Target  thrust 

C( 1630) 

ATURNT 

Maximum  transverse  acceleration  of 

target  in  terms  of  g 

i 

C  (1632) 

VDELX 

Relative  velocity  c-f  missile  to 

target  in  X  direction 

C (1633) 

VDELY 

Relative  velocity  of  missile  to 

target  in  Y  direction 

C (1634) 

VDELZ 

Relative  velocity  of  missile  to 

1 

1 

target  in  Z  direction 
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TABLE  IV  (Continued) 


— 

Common 

Location 

Variable 

Name 

- 1 

Def inition 

C (1635) 

RDELX 

Range  difference  between  target 
and  missile  in  the  X  direction  of 
the  earth  fixed  coordinate  system 

C  (1636) 

RDELY 

Range  difference  between  target 
and  missile  in  the  Y  direction  of 
the  earth  fixed  coordinate  system 

C (1637) 

RDELZ 

Range  difference  between  target 
and  missile  in  the  Z  direction  of 
the  earth  fixed  coordinate  system 

'  C  (1644  ) 

ATARG 

Acceleration  of  the  target 

C (1647) 

VTARG 

Velocity  of  the  target 

C ( 1648 ) 

1 

RTXED 

The  X  component  of  the  velocity  of 
the  target  in  the  earth  fixed  co¬ 
ordinate  system 

C ( 1651) 

P.TXE 

The  x  coordinate  of  the  position  of 
the  target  in  the  earth  fixed  co-  | 
ordinate  system 

C ( 1652} 

RTYED 

The  Y  component  of  the  velocity  of 
the  target  in  the  earth  fixed  co¬ 
ordinate  system 

' C  (1655) 

RTYE 

The  Y  coordinate  of  the  position  of 
the  target  in  the  earth  fixed  co¬ 
ordinate  system 

C (1656) 

RTZED 

The  Z  component  of  the  velocity  of 
the  target  in  the  earth  fixed  co¬ 
ordinate  system 

C (1659) 

RTZE 

The  Z  coordinate  of  the  position  of 
the  target  in  the  earth  fixed  co¬ 
ordinate  system 

0(1660) 

VTXL 

The  X  component  of  the  velocity  of 
the  target  in  the  earth  fixed  co¬ 
ordinate  system 
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TABLE  IV  (Continued) 


Common 

Location 

Variable 

Name 

Definition 

C (1661) 

VTYE 

The  Y  component  of  the  velocity  of 
the  target  in  the  earth  fixed  co¬ 
ordinate  system 

C (1662) 

VTZE 

The  Z  component  of  the  velocity  of 
the  target  in  the  earth  fixed  co¬ 
ordinate  system 

C (1663) 

VDXB 

The  X  component  of  the  acceleration 
of  missile  in  the  body  coordinate 
axes 

C ( 1664 ) 

VDYB 

The  Y  component  of  the  acceleration 
of  missile  in  the  body  coordinate 
axes 

C  (1665) 

VDZB 

The  2  component  of  the  acceleration 
of  missile  in  the  body  coordinate 
axes 

C  (1666) 

BDIVE 

Initial  pitch  orientation  of  the 
aircraft  (missile  assumed  oriented 
parallel  to  aircraft) 

C  (1667) 

RSLANT 

Initial  slant  range 

C  (1668) 

RXO 

The  X  component  of  the  original 
launch  point  of  the  missile  in  the 
earth  fixed  coordinate  system 

C  (1669) 

RYO 

The  Y  component  of  the  original 
launch  point  of  the  missile  in  the 
earth  fixed  coordinate  system. 

C  (1670) 

RZO 

The  Z  component  of  the  original 
launch  point  of  the  missile  in  the 
earth  fixed  coordinate  system 

C  (1672) 

BPSITD 

The  angular  rate  of  turn  of  target 

C ( 16  75) 

bpsit 

The  total  angle  through  v/hich  the 
target  has  turned  in  degrees 
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TABLE  IV  (Continued) 


Common 

Location 


C  (1676) 


C  (1677) 


C  (16  78 ) 


Variable 

Name 

Definition 

ANGX 


ANGY 


ANGZ 


The  X  component  of  the  accelera¬ 
tion  of  the  missile  in  terms  of 
g  with  respect  to  body  axes 

The  Y  component  of  the  accelera¬ 
tion  of  the  missile  in  terms  of 
g  with  respect  to  body  axes 

The  Z  component  of  the  accelera¬ 
tion  of  the  missile  in  terms  of 
g  with  respect  to  body  axes 


0(1700) 

CFA11U 

Derivative  of  CFA11 

C  (1703) 

CF  All 

Cosip  Cos0 

C  (1704) 

CFA12D 

Derivative  of  CFA12 

C (1707) 

CFA12 

Sinip  Cos0 

C  (170a) 

CFA13D 

Derivative  of  CFA13 

C  (1711) 

CFA13 

-Sin0 

C (1712) 

CFA21D 

Derivative  of  CFA21 

C  (1715) 

CFA21 

Sinip  Cos9  +  Cosip  Siri9  Sinip 

C (1716) 

CFA22D 

Derivative  of  CFA22 

j  C { 1719) 

CFA22 

Costp  Cos<P  +  Sin ip  5in0  Simp 

C ( 1720) 

CFA2 3D 

Derivative  of  CFA  23 

C (1723) 

CFA2  3 

Cos6  Sinip 

C  (1724) 

CFA31D 

Derivative  of  CFA31 

C  Q727) 

CFA31 

Cost)'  SinQ  Cosip  +  Simp  Sin4» 

C (1728) 

CFA32D 

Derivative  of  CFA32 

TABLE  IV  (Concluded) 


Common 

Variable 

Location 

Name 

Definition 

C (1731) 

CFA32 

Sin^  Sin0  Cos<>  -  Cosiji  Sin<p 

C (1732) 

CFA33D 

Derivative  of  CFA33 

C ( 1735) 

CFA33 

Cos6  Cosiji 

C  ( 1736 ) 

WPD 

d (WP) /dt 

C (1739) 

WP 

Roll  rate  of  missile 

C  (1740) 

WQD 

d  (WQ)  /dt 

C ( 174  3) 

WQ 

Pitch  rate  cf  missile 

C (1744) 

WRD 

d  (WR)/dt 

C (1747) 

WR 

Yaw  rate  of  missile 

C (1748 ) 

FMIX  ' 

Missile  moments  of  inertia  about 

C (1749) 

FMIY 

► 

the  X,Y,  and  Z  missile  body 
axes  in  flug-feet2 

C (1750) 

FMIZ 

md 

C (1751) 

CRAD 

Conversion  factor  (from  radians 

to  degrees) 

C (1752) 

BPHIO 

Initial  roll  angle  of  missile 

C(1753) 

BTHTO 

Initial  pitch  angle  of  missile 

C (1754) 

BPSIO 

Initial  yaw  angle  cf  missile 

3 . 2  Subroutine  Call  Sequence 

The  subroutine  call  sequence  is  determined  by  the 
order  in  which  these  subroutines  are  identified  in  the 
data  card  assembly.  A  data  card  is  identified  by  the 
program  as  a  subroutine  call  by  the  number  2  located  in 
column  2.  The  iaentif ication  number  of  the  subroutine 
may  be  called  as  either  (MODNO (NOMOD) ]  or  [ XMODNO (NOMOD) ] 
by  the  program.  This  integer  must  be  right  adjusted  to 
column  25  on  the  card.  Table  V  shows  the  identification 
numoer  and  the  subroutines  called  in  the  example  problem. 
If  other  routines  are  required,  they  will  be  found  listed 
with  their  identification  numbers  in  subroutine  AUXSUB. 

3.3  State  Variables 


The  state  variables  within  this  six-degree-of-f ree- 
dom  simulation  program  are  defined  in  the  initialization 
subroutines  (modules) .  These  variables  are  identified 
through  the  IPL  table  which  also  defines  the  location  of 
the  state  variables.  Only  these  variables  are  integrated 
by  the  integration  routine  AMRK.  Other  variables  found 
in  the  program  which  are  derivatives  are  not  state  vari¬ 
ables  by  this  definition.  A  listing  of  the  sequence  num¬ 
ber,  IPL  numbers,  and  variable  names  are  found  in  Table 
VI.  The  listing  is  for  the  program  when  it  contains  the 
high  frequency  autopilot  and  actuator. 

In  the  event  a  location  is  defined  as  a  state  vari¬ 
able,  the  following  convention  must  be  observed: 

C(J  +  3)  State  variable 


then 


C ( J)  is  the  derivative  of  that  state  variable. 
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TABLE  VI.  STATE  VARIABLES  AND  DERIVATIVES 
NAMES  AND  LOCATION  CODES 


Sequence 

No. 

N 

IPL(N) 

C  <IPL  CO  ) 

Module  Location 

1. 

1 

424 

BTHTGD 

SI,  S1I 

2. 

427 

BTHTG 

SURFACE  AND  PLAT¬ 
FORM 

3  „ 

2 

428 

BPSIGD 

4  . 

431 

BPSIG 

5. 

3 

800 

BPHISD 

Cl,  C1I 

6. 

803 

BPHIS 

HIGH  FREQUENCY  AUTO¬ 
PILOT  INITIALIZATION 

7. 

4 

804 

WQSDD 

MODULE 

8. 

807 

WQSP 

9. 

5 

808 

WQSD 

10. 

811 

WQS 

11. 

6 

812 

WRSDD 

12. 

815 

WRSP 

13. 

7 

816 

WRSD 

14. 

819 

WRS 

15. 

8 

820 

ESUMOD 

16. 

823 

ESUMO 

17, 

9 

824 

ESUMED 

18. 

827 

ESUME 

19. 

10 

828 

EZSDD 

20. 

831 

EZSP 

21. 

11 

832 

EZSD 

22. 

835 

EZS 

TABLE  VI  (Continued) 


Sequence 

NO. 

N 

IPL (N) 

C (IPL (N) ) 

Module  Location 

23. 

12 

836 

EYSDD 

24. 

839 

EYSP 

Ci,  C1I 

25  . 

13 

840 

EYSD 

HIGH  FREQUENCY 

AUTOPILOT 

26. 

84  3 

EYS 

INITIALIZATION 

MODULE 

27. 

14 

880 

EZSSD 

28. 

883 

EZSS 

29. 

15 

8  84 

EYSSD 

30. 

887 

EYSS 

-  —  — 

31. 

16 

1100 

BDELTD(l) 

- - - — - 

32. 

1103 

BDELT(l) 

C4 ,  C4I 

33. 

17 

1104 

BDELTD ( 2 ) 

.  HIGH  FREQUENCY 

MODULE  (ACTUATORS) 

34. 

1107 

BDELT (2) 

35. 

18 

1108 

BDELTD (3) 

36. 

1111 

BDELT (3) 

37. 

19 

1112 

BDELTD (4) 

38. 

1115 

BDELT (4) 

i 

39. 

20 

1124 

BDLTDD  (1) 

j  40. 

1127 

BDELTP ( 1) 

1  41. 

21 

1128 

BDLTDD (2) 

42. 

1131 

BDELTP  (2) 

43. 

22 

1132 

BDLTDD  (3) 

44. 

1135 

BDELTP (3) 

45. 

23 

1136 

BDLTDD (4) 

46. 

1139 

BDELTP (4) 
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iAMX  VI  (Continued) 


Sequence 

A 

t«W  • 

_ 1! _ L 

ITL  i  N ) 

r 

C  ( i PL l N  J  ) 

Module  Location 

47. 

n 

1496 

LIMPD 

A3I 

48. 

■ 

1499 

LIMP 

ENGINE 

49. 

25 

1600 

VXED 

Dl,  Dll 

50. 

1603 

VXE 

TRANSLATIONAL 

1 

DY  NAT-1  ICS 

51. 

26 

1604 

VYED 

52. 

1607 

VYE 

53. 

27 

1612 

RXED 

54  . 

1615 

RXE 

55. 

28 

1616 

RYED 

56  . 

1619 

RYE 

57  . 

29 

1620 

RZED 

5a  . 

1623 

RZE 

59. 

30 

1640 

VTARGD 

60. 

1643 

VTARG 

61. 

31 

1644 

BPSITD 

62. 

1647 

BPSIT 

63. 

32 

1648 

RTXED 

' 

64  . 

1651 

RTXE 

C  5 . 

33 

1652 

RTYED 

C6  • 

1655 

RTYE 

67. 

34 

1656 

RTZED 

68  . 

1659 

RTZE 
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TAULL  VI  (Concluded) 


■HHI 

C  (IPL(N)  ) 

Module  Location 

1700 

CFAIID 

D2,  D2I 

1703 

CFAII 

ROTATIONAL  DYNAMICS 

1704 

CFA12D 

1707 

CFA12 

1708 

CFA13D 

1711 

CFAl  3 

1712 

CFA14D 

1715 

CFAl  4 

1716 

CFA22D 

1719 

CFA22 

1720 

CFA23D 

1723 

CFA23 

1724 

CFA310 

1727 

CFA31 

1728 

CFA32D 

1731 

CFA32 

1732 

CFA33D 

1735 

CFA33 

1736 

WPD 

1739 

V.T 

1740 

WQD 

1743 

WQ 

1744 

WRD 

1757 

V.'R 
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SUCTION  IV 


INPUT  KL JU I REMUNTS 


•i .  1  Initial  Conoitions 

In  order  to  simplify  the  input  data,  the  options 
which  had  existed  in  the  original  program  have  been 
eliminated.  It  is  assumed  by  the  program  that  all 
variables  not  initialized  are  automatically  set  equal 
to  zero.  Input  data  and  initial  conditions  are  entered 
into  the  program  by  entering  a  number  3  in  column  2  of 
the  data  card  which  identifies  the  type  of  information. 
The  name  of  the  variable  may  be  entered  in  column  3  to 
20.  Common  location  of  the  variable  must  be  entered 
right  adjusted  in  columns  21  to  25  and  the  numerical 
data  in  columns  31  to  45.  Figure  19  shows  the  position 
of  the  data  card  in  the  completed  program  deck  which  is 
ready  for  submission ,  as  well  as  the  actual  data  card 
format . 


In  addition,  since  the  seeker  is  generally  assumed 
to  lock  on  before  launch,  the  gimbal  angles  are  automa¬ 
tically  initialized  to  this  position.  However,  in  the 
cases  where  gimbal  angles  must  be  chosen  in  any  other 
position,  the  transformations  and  angular  displacements 
between  gimbal  axes  and  body  axes  coordinate  systems  are 
given  in  Appendix  I  and  Figure  T-1 ,  respectively . 

Initial  position  and  velocity  can  only  be  specified 
in  one  manner  for  simplicity.  They  are  specified  in 
terms  of  the  following  variables: 


BDIVL 

(in  degrees,  negative 
low  horizontal) 

nsuA::"’ 

l  - 

:  (fu-t) 

i.ALPl:A 

(367) 

(degrees) 

BALPHY 

»  j58) 

(degrees) 

VMACH 

(204) 

(Mach  number) 

when  orientation  be¬ 


lt  should  be  noted  that  the  program  when  used  to  simu¬ 
late  many  mis.sions  requires  only  that  subsequent  changes  in 
data  be  added  since  the  program  will  only  update  the  last 
data  set  for  the  next  run.  (Sec  "Program  Description.") 


5G 


Executive  Rout 


Description 
olumn  No. 
Format 
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SECTION  V 


PROGRAM  LISTING 

5.1  Complete  Six-Degree-of-Freedom  Program  Listing  with 
Example 


c 

€•§•••01  FOCI  10  IB  US  10  WIN  FORTRAN  ANRx  INTEGRATION  ROUTINE 


CCPPCN  C  (4110)  , GRAPH 
EQUIVALENCE  ICI 266 J.'HNIN  l« 
C  ICI  256W,N  It 

u  (cuoooi.r  u 

C  (Cl  2012 ULSTER  It 

C  ICI  I02JJ.0P01NI  I» 

C  (C(  lUTi.NCOUT  It 

C  ICI  2665WEU  It 

EQUIVALENCE 

EQUIVALENCE  ICI 29 7  )  I  A tt  It 
QIFENSICN  CRAPH|?,2»,TIHE(2I 

DIMENSION  VI  ABLE! 2*15)  « 

C I FEM  (CN  VAR  I  10)1  t 

EQUIVALENCE  ICI  1 98C  1  *RN  I 
ECU  I  VALENCE  ICI |96l  I, RUT  I 
EQUIVALENCE  (C  (  198  2  J  ,PL0TN4  I 
EQUIVALENCE  ICI  1983UPLOTN2I 
EQUIVALENCE  ICI 1984)»nPL0T  I 
IMEGtR  cpoint 
Integer  plotno 
c  integer  CPf 

external  Auxsua 
DC  22I*l#43IO 
22  Cll)*0. 

SP1*2 

CALL  CCUNTV 
CALL  2CR0 

1000  tCMlNUE 

1001  If  (HCTNC.LEtO.  IGGTO? 

If IREPPLI-GT*0. ICOTOT 

C 

c  rspplt 

c 

c 

If  IREfPlT.GT.-i.OI  NOOUT  -  U 
NPlCT.O 
7  CALL  ClNPIl 
1C  AS  E  *0 

If IRKUTTAtGT.O.OI  NPT-I 
LSTEP  •  STEP 
NPlCU-PLCThA 
NPLCr2*fL0TN2 

ncplct-plotno 

1002  call  suali 

1001  CALL  AUXI 

1004  call  suti2 

100>  CC  *0  I  *  2 »M 

Jr-  -II 

ELI  1°  I J ♦ 1 1 
EUI  !“1  InCI J<2I 

variu  •  ci jo i 

60  DEMI)  R  CIJl 
VARIU  *  T 

1006  CALL  AUXSUA 

1007  NJ*N-1 

CALL  AFRKIAUXSUAl 

1008  00  SO  I  •  2,  N 

j  •  iPLd-ii 

90  C(J4l)  •  VAR  I  I ) 


ICI2663) *HMAK  ICI2664I*0ER  I 

ICI2S62)*IPL  I C I  2  969  I  *  VAR  I 

<Ct20li)  ,K$TfP  ),  1C  I 201Q I • 1 Tf P  I 
1C  I  2008  I  *  PLOT NQ I v  ICI 2C09I ,NC*LOT I 
ICI  2^251  ,  TlHE  I*  1C  (2125)  tVUALfl 
1C  I  2022  I »OPtNlO) •  (CI2006I .REPPLT) 
ICI276S) *  EL  )»  1020071. PTLCSSI 
(Cll971ltRIU  It  I C  I  1972  I  t*M»T  1  At 
ICI  1974  1  *NJ  It  ICU97»1*NPT  I 

1PU  100)  •  OEM  101 1 

ELI  100 )  •  Eul 1001 


•  0,  USE  NEV  NO. 4,7  IOISCARO  OLO) 
1.  USE  OLO  PLUS  THOSE  AUDEO 
•1*  USE  NEkt  NO.  7  IOISCARO  CIO) 
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t  •  VAfttl  » 

1009  call  sueo 

\f  t  KST£P  .50.  n  CO  10  1C07 
DC  155  JV-*,N 
115  VAft<JVUO. 

CHI  SB 
CALL  PRCCES 
CALL  «€S£T 

If  (LSI  EP.EQ.5. Oft. L STEP, 6Q.7 .OR. NOPLOT .EQ.QIG0T05 
CALL  T  tMVi  OILT  » 

I  »f  <6*9MC£LT 

96  fC**At(LH  *  17HST  ART  PLOTTING  ATF14.71 

lESSPt -PTLES5 

CPCIM •OPCIM-LESSPT 

Call  PICT4  ( GRAPH, CPQ[HT*VLA8LE»TlME »NPL0T4 I 
CALL  PLCT2  IHPICT2) 

CALL  PLCTN  IaOPlOT  ) 

CALL  MhEVICELT) 

WR  ITCI6.91  JCEtT 

97  F  C  &  *  A I  I  ih  ,  i  8hplCT  T  ING  ENDED  AT  F14.7) 

IF  1 (RSI .Cl -0. 1  .ANC.  9n.EQ.RnM.AnO.  lSTEP.EQ. 2  ICO  10  70 

5  CC  TC  < lOCO.lO0l.lCO2. 10QJ,  1004. 1005, 1006, 1007, 1008, 1009, tOlOl, 
l  LSTEP 

1010  If  (CPTMO.GT.O.  > 

lkaRfT£(6,6l(  l,C(  ll,C(t+l)»C(  M2  )  »C  (  l ♦ 3) *CT I ♦* 1 »C 1 1 ♦! ) ,C 1 1 ♦* I  % l *1 , 3$ 

•  10,7  I 

6  fCPPAfl  1HI/U5,  IP7C15.7)) 

70  CALL  S8 

CALI  EXIT 
EKC 

BLCCA  CAfA 

COKh  /KCXCP/NCXI6I 

•  /CXABG/ALPI 7»rAhI8) 

•  /CXf IN/CXI 5fc  ) 

c •  i  a  *cx/ 1 ,e, o.o, o.o/ 

CATA  Al P/0.  .4. ,6., 8., 10., 15. .20./ 

CAT*  AK/0.,.7,.9, J  .0  25,  1.2, 1.5,2. 0,2. 3/ 

C  A  1  A  CX/ 

1  0.0, 0.0, .18, .12, .<6,  1.56, 2. 66, 

2  0.0. 0.0, ,18, .12, ,66,1.56*2,66, 

)  O.C.C.C. .221.15, .71, 1.7,3. 27, 

4  0.0. 0.0, .23a. 3), .79, 1.8, 3. 2* 

1  0,0.0. C». 21,. 14, ,76, 1.65,3.2*, 

6  0.0. 0.01.13, .25, .51,1. 11,2.06, 

7  0.0,0.0,.;0,.2l,.!e,.75,l.67, 

8  0.0, 0.0,  .09k. 19, .  35, .82,1.57/ 

EKC 

81  CCA  CATA 
COPPCN  /NC2/hCNl4l 

•  /C2 AftC/Al p|  e  1, ANI 9) 

4  /C2 PUh/CN I  72  I 

CCPPC6/6CC2/NCCI4I 

•  /CCARC/AiPCm,  APCI6) 

•  /CC2*t*/DtNl 3*1 
CCPPC6/^C^/nCh<  A  1 

•  /CPARG/  ALfl8),Bhl9  ) 

•  /CKFUK/CMI 72 1 
CCPPCN/NCC*/NOM( 41 

«  /CCr.AftC/AlfOI6),ftKD(  1 1 

•  /(C»«fu?!/CCHI  3*1 
CATA  KCN/8, 0,0.0/ 

OAT  A  ALP/O. ,2..4.,6.,8., 10., 15., 20./ 
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OAT*  **/0  •«•?,. 9,1  .05, l.i,  1.5,  1.7, 2.0,2. 1/ 
CAT  A  CM/ 

•  0.0, .65, 1.36 ,2. 21*3. 12*4. 13, 7.06, 10.71, 

•  0. 0,. 66,1.41 .2. 2 5, 3. 17, A. 19, T. 18.1O.0A, 


• 

0.0, 

.69, 

1.46, 

2.32, 

3.27. 

4.31, 

7.  36, 

11.  lit 

• 

0.0, 

.77, 

1.62, 

2.56, 

3.59, 

4.72, 

7.98, 

12.02, 

• 

0.0, 

.77, 

1.62, 

2.67, 

3.6C, 

4.74, 

8.C2, 

12.11, 

•  0.0, .71,  1,5.  2.30,  3.36,  6.47,  T.6  ,  11.60, 

•  0.0,  .7,  1*49,  2.37,  3.34,  4.40,  7.59,  11.61, 

•  0.0, .6Q,  1,47,  2.34,  3.30,  4.36,  7,56,  11,61, 

«  0.0, .67,  1*43,  2.27,  3.21,4,26,  7.42,  11*47/ 

OATA  HCC/6, 6^0,0/ 

C  A  r  A  ilPC/O.I  4.,  e.,  12..  16.,  20./ 

OA  f  A  A»C/0.0».7«  .4,  1.1,  1.5,  2.3/ 

CAfA  CCN/ 


0  .0, 

•  06, 

.30, 

.  75, 

1.40, 

2.  20, 

0.0, 

.06, 

*  30, 

.73, 

1.40, 

3.20, 

0.0, 

.06, 

*  30, 

.79, 

1.45 

.  2,2 6. 

0.0, 

.06, 

*30. 

.79, 

1.45,  2.26, 

0.0, 

.06  , 

.30, 

.77, 

l.*3, 

2.2*, 

0.0, 

•  06. 

.30, 

.it. 

1.41, 

2.22/ 

C  A  ?  A  MO'/8,9V0,0/ 

CATA  Alf/0.»?«,  4 • ,  6.,  8.,  10. »  1$.,  20./ 

CAT#  8H/0.0,  .7,  ,9,  1.05,  1.1,  1.5,  1..,  2.0,  2.3/ 


CATA  CM/ 

•  0.0*  -1.16.  -2.48, 

-3.97.  -5.62, 

-7.51, 

-12*83,  -19.28, 

0.0, 

-1,19,  -2.54, 

-4.06,  -5,76, 

-7.68, 

-13.10.  -19.64, 

c.o. 

-1.25,  -2.66, 

-4.23,  -5.9e, 

-7.97. 

-13.54,  -20.26, 

0.0. 

-1.591  -2,35, 

-5.29,  -7.40,- 

9.75,  - 

16.24,  -23.67, 

0.0, 

-1.66,  -3.48, 

-5.40,  -7.66, 

-10.09 

,  -16.75,  -24.56V 

0.0, 

-1,07,  -2.31, 

-3.74,  -5.35, 

-7.22. 

-12.48,  -18.91, 

o.o. 

-1.114  -2.39, 

-3.65,  -5.50, 

-7.40. 

-12.70,  -19.12, 

o.o, 

-1.15,  -2.47, 

-3.97.  -5.65, 

-7.57. 

12.90,  -19.29. 

0.0, 

-1.14,  -2.46, 

-3.96,  -5.63, 

-7.54, 

-12.81,  -19.09/ 

CATA  NC>»/6,6VO.O / 

CATA  A  L  F  0/0 • V  4.,  0.,  12.,  16.,  20./ 
CATA  0MC/O.,  .7,  .9,  1,1,  1.5,  2.3/ 
CAT  A  CCM/ 


0.0, 

-.13, 

-.64, 

-A. 82. 

-3.45. 

-5.40, 

0.0, 

-.13. 

-.64, 

-1.82, 

-3.45, 

-5.40, 

0.0, 

-.13, 

-.73. 

-1.93, 

-3.65, 

-5.65, 

0.0, 

-.14, 

-.77, 

-1.49. 

-3.70. 

-6.00, 

0.0, 

-.13, 

-.73, 

-1.93, 

-3.65, 

-5.65, 

0.0, 

-.13. 

-.73. 

-1.91, 

-3.65, 

-5.65/ 

BMC 

6LCCK  CATA 

CCMMCN  /NCXO/NCXOI 21 

•  /CxCAAC/AMI  8  I 

o  /CX0TUN/CX01 6 1 

OITA  hCXC/8,1/ 

CATA  AM  /O. », 7, .9, 1.025, 1.2,1. 5. 2.0, 2.3/ 
CATA  CXG/.24, .24, .2 3,. 6 6,. 76, .91, 1.0, .94/ 
CMC 

BUCK  CATA 
COPPCh/  MCN2/NCN2141 

•  /  CN2AAG/ALPI 6),  AMI  3 1 

•  /  CN2fUS/CN2l 10 ) 

DATA  6CM2/6, 3,0,0/ 

CATA  AlP/O.,  4.,  8.,  12. ,  U.»  20./ 

CATA  AM  /O..  .7.  2.3/ 

OATA  CM2/ 

•  0.0,  ,17,  168V  1.50,  2.60.  3.90, 
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a  0.3.  .»».  I'll) 

a  0.0.  .17.  l*50'  1’e0*  1 

CSC 

ti.CC*  C»T» 

rrrrO,  /kCYJ'SCYJIO 
,  /CY2AKG/ai  *<  4  •  • 

4  /CY2FUN/CY2U8) 

cat*  SCY2/6. 1.0.0/ 

CAT*  AlP  /0..  9*>  ,J**  ' 


CAtA  AH  /Q .Q  .  i 2 . 3/ 

CATA  CV2  / 

0  •  0  «  **0T  t  *.30, 

-t6S  •  1  *3  * 

O.o*  **»0T  ♦  3C* 

-.65.  •M» 

C.<5»  -•<>/»  -.30. 

->■65.  **Ul» 

e*c 

81 CC*  C  al  A 

CCK-Cs/  NCL  2/  -lCL2l*l 
■  /Ci.2»sG/*Cpl  61.**!** 

a  /Cl2H>/Cl2t  2*' 

C»T*  *CL2/6.a.O.O>  a>/ 

data  alp/0. Oi  A. I  8  •  * 

Cara  av/0..  .2.  .0. 

cat*  Ct  2  /  ,, 

•  0.0.  .02.  .ot.  .1*.  ***’ 

a  0.0.  .02.  .0?.  .!»•  •*’' 

a  0.0. .03.  -12.  .273.  **'5* 
a  0.0,  .02.  .07.  .1*.  .13.  •»« 

CSC 

8i.CC*  C*T* 

CC**C*a*Cci7NCC»t*3 

,  /CcJaBC/ACPI 4 I. 

<  /CL3CU*/CCl( 1*1 

oali  8CC3/6 .5.0.0/ 
cata  acc/0..  •«>  '*•»  1 

cat*  >*/0.,  >7.  2.3/ 

.C*0*o“c22..0*>.  -08.  •{**•  •“** 

.  0.0  .022.  .0*3,  .0*.  .»*»>  •***? 

a  0.0. .022. .0*3  .  -08.  .1*1.  .2*’/  \ 

t*C 

8I.CC*  cal* 

<cc**C*/*U0/*CN(l»Jjn.cMjoJ)(tMjoluCHa.,«|.CN»Oil.  Ch*H2». 

*  C.nM  ?2> »  0  8  (  J2> 

«  /CfCARC/  ALP(  8  ItOdH  A  l  .AM  8) 

0ata  KCi./o.aa*. 0.0.0/  „  , 

cata  Cff  / •  20 • .  *10;.  ‘6*: 

oai a  ak/0.  ..7.  *Y ,  1.08.  3.2.  »«■*•  a.'1**' 

CATA  CM/  .  .  204.-.  194.-.181* 

•  -.222.-.218.*.210"*.  •  *  110. -.181  . 

a  ..2CJ.-.20A.-.203.-.201.  •  ,J08,  <i07, 

:  -j?r.  :I25: :»!:  *.»»:  -»5'  -”7' 

.e!:;^“2Ui-.2r -.m.-.jo.  :;7o:;::r9a.;ae»; 

a  -.206. -.20*.-. 203.  .203.  •*  .206.  .207, 

:  :%r.  :\\l: :»?: :«!: 

Cata  C*  3/  2J2, -.  220. -.207. 
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J 


r/ 

i 


« 

.  236  , 

.237*  .233, 

•  2  3*« 

.240 

•  240, 

.240,  .236, 

• 

•  252, 

.255  ,  .256. 

•  260, 

.262 

•  265» 

.268,  .269/ 

0 A 1  A  CK4/ 

4 

-.276, 

-.273, -.269, 

-.264,- 

.260 

-.255,- 

.241,-. 226, 

• 

-.260, 

-.259, -,257, 

-.254,- 

.251 

-.2*.9,- 

.2)9, -.229, 

■» 

.260  , 

,2611  .262, 

.26), 

•  264 

.264, 

•264,  .261, 

4 

.276  , 

•279V  .252, 

.283, 

•  288 

.285, 

.294,  .295/ 

C  A  I  A  ChS/ 

-.246, 

-.243, -,240, 

-.236,- 

.232 

-.227.- 

.215, -.202, 

-.230. 

-.229,-. 227, 

-.225,- 

.222 

-.220,- 

.212, -.20), 

■  230  . 

•232.  ,2 32, 

.234, 

•  234 

.2  34, 

.234,  .2)1, 

.246, 

.249,  ,252, 

,255, 

•  257 

.258, 

.262,  .263/ 

CATA  CN6/ 

-.166, 

-.1651-. 163, 

-.161,- 

•  158 

-.155,- 

. 146, -.138, 

-.152. 

-.  15U-.  149. 

-.149,- 

.147 

-.145,- 

.139, -.1)4, 

.  152, 

.152.  ,153. 

.153, 

•  154 

.154, 

•154,  .15), 

.168, 

.170%  .172. 

.173, 

.175 

.176, 

.179,  .180/ 

DATA  Chit 

-.131, 

129,-.  127, 

-.125,- 

.123 

-.121,- 

.115, -.108, 

-.115, 

-.114, -.113, 

-.112.- 

•  111 

-.109,- 

. 106,-. 102, 

•  115, 

.116*  .116, 

•  116, 

.116 

.117, 

.116,  .115, 

.131  , 

•132.  .134, 

•  135. 

.136 

.137, 

,139,  .140/ 

OATA  CK»/ 

-.123, 

-.122*-.  120, 

-.118,- 

.116 

-.114,- 

.108, -.106, 

-.107. 

-.106*-. 106, 

-.105,- 

.104 

-.103.- 

.099,-. 095, 

.107, 

.108,  -108, 

.109.  . 

109, 

.109,  . 

109,  .107, 

.12-3, 

•125,  .126, 

•  127. 

128, 

•  129,  • 

131,  .132/ 

CKO 

bicck  cata 

CCrPQ*  /KClC/NCl  t  6 } 
n  /ClCA*G/AiFl6),  (j£F|4»,  AM  131 

•  /ClCPuN/tl lt/4 I,  Cl2t24>,  amo 
C  A  t  A  NCL  /&,4 , 3.0,0,0/ 

DATA  Al?  /Q.,  A. ,  8 . •  12.,  16.,  20./ 

CATA  C&F  / -  20 . «  -10.,  10.,  20./ 

CATA  AP /O .0 •  1 7 , 2. 3/ 

CATA  CLl  / 

•  1**  * -.1*9.-. Ml,-. 138,-. 138, -.13*, 
•-.131 .-.128,-. 125,-. 12*,-. 12 3, -.121, 

«  .131.  .123.  .123,  .124,  .123,  .121, 


•  .144, 

i  .149.  .141, 

.138,  .136,  .134/ 

DATA 

Cl  2  / 

•-.144, 

,-.149, -.141, 

-.  1  38, -.  1  )8,-.  134, 

•-.1)1, 

,-,  1  20  .-.  125, 

-.I24.-.123.-.121* 

«  .131, 

,  .128,  .125, 

.124,  .123.  .121, 

•  .144, 

,  .149,  .141, 

.130,  .138,  .134/ 

CAIA  Cl 3  / 

•-.148, -.143, -.143, -.140, -.138, -.136, 

•-. 144, -.138, -.136, -.13 3, -.130, -.130, 

•  .146,  .130,  .136,  .133,  .130,  ,130, 

■  .U8,  .US,  .143,  .140,  .138,  .136/ 

€NC 

eicc*  CATA 
CC**'C4  /fiCHC/Ne>M6) 

•  /C»CAftC/ALP(8l,CEFf4lpAlM8) 

•  /CrCfUS/CMl ( 3  2  I, CM 2 1 32  I, CM3 1 3  2  I ,CH 4 1 32 ] »CH5 ( 32 1 ,C*6 1 32 ) ,C*7 1 32 1 , 

•  C'8 < 32  » 

OATA  KCP  /8 .4,8^0, 0,0/ 

LATA  AlP/O.  ,2., 4. ,6..  8., 10. *15,  ,20./ 

CATA  CEP  / • 20 ..  -10.,  10.,  20./ 

CATA  AP/O., .7,. 0,1.0 5, 1.2, 1.5, 2 .,2.3/ 

OATA  t*l/ 
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••i.h  it.iollHW' 


t 

.935 

1937) 

.937, 

.938, 

.930, 

.939, 

.944, 

.940, 

2 

.  664 

*  84  4k 

*064, 

.365, 

•  064  , 

.  e«4, 

.670, 

.874, 

3 

-.864 

*  *  B  5  4  ,  — 

.804, 

-.e65. 

-.864, 

-.864,- 

.870,- 

*874. 

4 

-.938 

-•937*.- 

.937, 

-.9)0  , 

-.930, 

-.9)9,- 

.944 

.948/ 

CATA  C*2/ 

1 

-.933 

-.937,- 

.  9  37  , 

-.938, 

—.9)0, 

-.939,- 

.944,- 

*948, 

2 

-.864 

-.664)- 

.864, 

-.665. 

- . 064, 

-.865.- 

.670  ,- 

.874  , 

3 

.864 

.06*., 

.  864  , 

.665, 

.  e64. 

,065, 

.  670. 

.  874  , 

4 

.938 

.9)7, 

.  9  3  7  , 

.  9  50  . 

.9)0, 

.939, 

.944  , 

.948/ 

CATA  03/ 

1  -1.123.-1.  1)7,- 1.137,  -l.|  J  9,- 1.1 38.-1.  14.-  1.1 44,-1.  152, 

2  -1.0ft4,-l.C64, -I.C64,-!. 0fc4,-1.064.-l. 066, -1.072,-1. 077. 

3  1,064,  L.064,  I.C64,  1.064.  I.C65,  W072,  1.077, 

4  1.12).  1.07.  1.07.  1.120,  1.13e.  1.14,  1.144,  1.152/ 
cata  c>4/ 


t 

-1.277 

,-1.277,-1.277,-1.277,-1.277,-1 

.275, 

- 

1.285, 

-1.29, 

2 

-  1.204 

,-i; 203, -1.20 3,-1.203. -1.204,-1 

.205  , 

- 

1.211 

-1.216 

3 

t  -20  4 

,  1.203,  1.203,  1  .  2 C  3 »  1.204,  l 

.  205  , 

1.211, 

1.216 

4 

1.277 

,  1.277,  1.277,  1.277,  1.277,  1 

.279  , 

1.285, 

1.29/ 

CATS  Cl*  5/ 

1 

-  .  796  , 

787  , 

- 

.  /fl  6  * 

2 

-.7/0, 

-.Tu.-./is.-.m.-.m.-.m.- 

.  711  , 

- 

.711  , 

3 

.  720  , 

) 7 1 8 1  .715,  .715,  .713,  .712, 

.711  , 

.711, 

4 

.  796. 

.793.  .791,  .789,  .7.4,  .747, 

.78  7  , 

.  7  86/ 

C  A  I  A  C *6/ 

l 

-.621. 

-.620  k -.6 19, -.6  IT, -.6  15, -.6 14.- 

.612, 

- 

.61), 

2 

-.545, 

-.4**4-. 4*3. -.4*1, -.4*0. -.4 39.- 

.53  7  , 

- 

.5)4, 

3 

.545. 

.S'.*.  .5*5,  .4*1.  .5*0,  .439. 

.  53  7  , 

.43*. 

4 

.621. 

,620k  .619.  .617,  .615,  .614, 

.6t2, 

.613/ 

c * r a  07/ 

1 

-.566  , 

-.56 4% -.53  3, -.582  ,  -  .  580 ,-.5  79 »- 

.477. 

- 

.578, 

2 

-.610, 

-.40B.-.  40  7,-.->0  6,-.*03.-.5C*.- 

.502, 

- 

.50), 

3 

.510. 

.409.  .407.  .406,  .404,  .4C*. 

.502, 

.50), 

4 

.586, 

.48*4  .493.  .581,  .480,  .479, 

.577  , 

.570/ 

CATA  09/ 

1 

-.586. 

-.5044-. 583,-. 501,-. 590, -.379»“ 

.  378  , 

- 

.575, 

2 

-.510, 

-.506. -.50 7, -.50 6, -.505,-. 504,- 

.  502  , 

- 

.SC3, 

3 

.510, 

.506,  .507,  .504,  .505,  .504, 

.  502  , 

.503  , 

4 

.566, 

.46*4  .483.  .481,  .580,  .479, 

.578, 

.574/ 

CSC 

buck  cata 
CC^^C6  /nTi-/nTM21 
«  /TMA«G/ThA< 20> 

•  /  jHFuVThM  20) 

CAT!  Mh/13,1/ 

CAT*  7hA> 

•  0*.  .05,  ,1,  ,2,  .4,  « 0 ,  1,3*  1*S.  1.6,  1.0,  2,0*  2,0*  100,/ 

E  A  7  A  T  Y.t  J 

«  1.,  5C0Q  .  ,  4000.,  3800.,  3650.,  3000,,  3600.,  2000.  .HOO.t 

•  700.,  250. v  0.0,  0.0/ 

£KC 

BLOCK  CATA 

CO^CN  /MCLP/NLIAI 

•  /CL  P  ARC/ At.P  I  8  )  .  AKP<  7  ) 

•  /CL^FUK/CLP( 56  ) 

CGPKCS  /6O»C/S0l*l 

•  /C^QARG/ALC  l  6  ).AXQ{4 ) 

•  /C*CPuN/C*Qt  24) 

CATA  ML/fl. 7,0,0/ 

CATA  AlF/O. ,2.l4.*6.,6.,10,fl),»2O,/ 

OAT  A  AKP/o. r. 7, .9,  1.05, 1,5, 2. 0,2.3/ 

CATA  CL»/ 
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t 

c 

c 

c 

c 

c 


.  •  •  _  a  *»  d  —  50 1  '*  1 0 •  5 6  i  •  3^  * 

.  -e.OT.'».*Ol-t.OT.-T.T».  ?  2!  _,.yi.-io.*4.-n.»i  • 

»A,  43  ,  “4 .  11  1 

CAT*  hC/t.'.tl.O/ 

CAT* 

CAT*  Al»C/0»  .  A 1  •  3.  ** 

CATA  Cl<i/  _  .. 

«  „  ->*  -.4  I  ,  - 

Sua^CUt  |N£  C«*OfTCt€fMC.0*f.CS.Mi.i.CC.OO.Wll..in« 

Ct^CKClOlO ) 
ecuivaLtNcetct ) 

t  irfcSSICMCf 20  )  ***  20)  *0' 
nPTftfcC-Hauf  FQ*  °F  CMfcCTOft 

cUm->u*»cv  Of  CffOtus  spot 

AA-CCfct®  L£ft  OUACPNI 
te-LCUt®  »(  CU4C»A\t 

CC-OPPE®  L6ff  QUACPNf 

?uX  otfotus  Spot  »  iHMHm 

CEfCCS-CEf ICS/St-J«*T» 

RT 1 • I  1 . ALOO . /®  r  T 

ACf «efPSf/i 7.29SI# 

YCf .HlfSA/t 1 «JS5t» 

UCf«'.f07l*ITOf-I<Of  > 

*Cf •OLF 

CETRAC-CET«lC/57.3SSf» 

ICf.SCMUQf  •■i»YOP«»J» 
tev«Ar.-cEi*>t»oefocs*»H 

iFiiCf.ce  ,c6i®ao>ociow 
C2T«-EEIPfO»*3 

Off RAC-OetOC4«Rtl 

02f RAC-CEF»AC»»3 
Ch.lt/3 

OELf*CfT®AC/T.N 

CACN’Ofir/A. 

VC  i -e  t » l 
CCVOt'l.ACl 
Al 1 )»0. 

15  ei 1 1“0. 

COOI  •  >  «icl 
CIK"C£T®AC“I  l"ll*Otl-F 

;j  COTOJO 

ix-lCAl ICJFRsO-ICIPAtOf )**3> 

CQ20Nt-l*2 

tetrt  i j.Lf*0*)dCT035 


CFIill  >  .tf.O.lM  tl-O. 

B< I l*C( 2  > *OZCN 

■I  F  C  6  <  I  l.GT.0.191  1 1  ■  0  • 

ccrc  jo 

15  A( 1 1«OI2>~02CN 

IF  (At  |  1.11.0*1 A(  n-o. 
tin<cii  i*ckn 
ifiei i  i.cr.o.iei f i-o. 
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23 

24 


25 

10 


40 


SI 

42 


50 

60 


C 

C 

c 


CC  10  10 

IflCI  ll«C(2H23.24*2S 
ec  t i--aes(0( 1 1 —ci 2 1 i 
4(11-0. 

CCTCBO 

4(11*0. 

B  (  I)  *0 . 

ccrc  30 

44 (  )  •  A  6  S  (  C 1  11-0(21  I 

B( 1  1-0. 

ccs  r  imjC 

44-0. 

B9-0. 

cc-o. 

cc*o. 

4<lS»U*0. 

e(tN»n%o. 

OC40I-1.1N 

aa.abSI  <AIII»»I  l*l>»o.B»0£tM»A* 

ea.ACii(eili*0(t«ll )•#» 

^ C - *es  ( (  a <  I  •  L«  i  ».* 1 1  *ln *  1 1  (•.3«0£1F  l*cc 

CC>ASSIielI«I.NI*8It  .  5«CEtMAQD 

ACE I-l. l.l».«C2TR«.» 

Tcrii»a*>c:  ‘CC'tc 

IF  lIOML.tt  •*OEU.«lf£<».*Uftl 

FCRMATI IBM  etlNO  8<N0£  *IP£12»** 

ccrc.o 

HS(rt(6.50l«T 

C8»EAii*-l. 

AA«0. 

B8«0. 

CC«Q. 

FGRfaMIJM  eatlKLCCH  »I  HANGS  EOOAl  T0 
RE tuHn 


t\c 

5<JA«CuUF.er.AUSStS,AM.V! 

CGFfCnCIAJIOI 

ECuIVAUEf.Ce  (Cl  103l.m» 

ti.uivAi.SNce  iciioooi.ti 

J-tFE  RCCui’.SC  STAN CARO  DEVIATION 
AF-IS  IFE  REt'JlREC  FEAN 

Ft  f  rwPtiTF  r  NORMALLY  0ISIRI6UTE0 


*•0. 


8-1. 

CC  50  IMil! 

VFURCNGlBI 


50  A«A«VFL 

V* ( A-6 I «S«A" 

rsiora 

EKC 

subroutine  cioi 

CCFFCN  CINJIOI 

ECUlVALENCElCITI.AH.t.1 


.1PE12.M 


RANDOM  NUMBER 
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CCUIVUEHCEICI  lll.Birjci 
E.LHvll  tNCE  let  1  2  ) ,  vu*> 
f  curv*(.ff.cElc,  20U  lV„At(.. 

€CUIV<lfNrj,C,  ,aiiJ  ' 

fCUIViU^CtlCI  16)5), 

Evtl  V AL  COCE  (Cl  6  l.HIil  1 
ECU(V<cf5CE(c,  lt2J|>(ue| 

ECU|V«1.  ENCf  let  5  )  ,1  111  ) 

ECuivilenceici 1  ll.yiif) 

EQUIVALENCE ICIltJ.ull 
Eci/i vi c f 6ce ( e 1 2000 1  n 
tCUIVAlEl.CEiei  1666  1,801  Vfj 

EQUIVALENCE (Cl  1667  l,«UAItf  > 

?'“!*‘l£NCC|E|a).SEcriiO) 

ir  t §,oC*€C*°*°’  "ETuaN 

IMMll.fC'.OI  OC  1C  10 
J  '•"‘C.ci.n  co  loo 

•  VI-ACN 

*{.*■*  au'f,r  •  ccsciBomei 
J  ■  Ji?NT  *  SINtfBOlVE) 

Kite  •  iBSIHIccl 
CC  1C  1 
10  CCM  thUE 
mi  •  o.o 
o  continue 

Jiii  *  SC"!'*'f?  «  XIU..2! 

VSCc*c  •  *0. Cd392»hlLL  ♦1117  % 

liH  fOO«»r  l2l.9En.4i 

‘"“•"“MU,  T  .HUu.zE.8lu.VSOw,0iy,u 

€^c 

SUQ^OUTfNC  CIO 
CCH^CNCU  MO  \ 

eeuiv4i6-.cE  ici  m ,vuii 

fcul  V41  EUCE  ICI  le  I,  »H  ( 

eeeiv4iE6CE(ci?652),H«|Nj 

ECU  I  VALENCE (Cl  5  I ,R  III  ) 

ecUfV4tthCE«r.l*|,HlLll 
equivalence. cm, ecuet 
equ  vale-, ceicm, wand; 

6CUlV4ie.-.CE(C(i),B4CI0S» 

EQUIVALENCE (C(6],AlS8or) 

ECUIVAlEnCEIC.VI.XSPOI) 

eeuivi-  .sceicuoi.vspoM 

ecu  ivin  r,ce  i  ci  i «)  ,or  i 

EcutvieEKCdcuovoi.r* 

ECUIVlLEr,Cf(C(466),Cxrrl 
ECUIVUEiiCf  tcut51l,«i*{> 

ECU  VALENCE  (Cu655  l,»rv£j 
6  Ctrl  VALENCE  1C(8  I  .SPCINIJ) 

<  ECIPC.EC.,.0,  sETu»N 

p:SS%:5?ai!r!ssi 

IfU.CC.g,  CO  10  4  °  * 

i  ccMinue 

1«0 

Ul  •  Uf  »  CHDT 
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IF! AILL .EC.CVOR.VlLF .EQ.0.0  I  GO  TO  10 
RILL  »  SORT  lXfL«*2*MLL**2l 
10  CCMlSUf 

ANG  ILL  «  AS  IKOIHIL/MLU 
8C  R  6  F  •  BCRF  *H  Hi  /  ICCO.O 
WANCF  ■  WANC*R  ILL /  10C0. 0 
SPwIO  »  RAC  1US  *R  ILL  /  loco  .0 
SPLENG  m  SPwIC/SlSCi AN G ILL  ) 

OTwC  »  RCi.OGi  ftCREF  » 

CTwCP  »  RLSLC I GCRtF  I 
C  T  if*-  A  €  E  *  RCMTitwASCF  ) 

CTr-REP  »  RCNCGIW.NCF) 

OffLft  »  -SKL ESC*A l SPOT 
CFCLRP  »  -S  Pn l C*  A  I  SPOT 
CMVF  aRCNOI  SPl  E N C < 2 . 0  I • A l SPOT 
CM  V  €  P  -  RLSCl  Y  PN  IC*2.0)*AISPQT 
2l  ASft»C7uC*CHAE‘ *C FOUR* OF IV£ 

YL  A  SR  -  C  T*CP  ♦  C  r  l-.rt  €  €  P*CFCLRP  *CFI  VEP 

3  CCMIMjE 

1234  f  CRPAK2X,  10E13.4  I 

»fi!TEl6«l?)4)xl| ,V  ILf.PH  IN.  RILL, HILL,  ANGIU.aOREP.eOAE ,  HAF-OF  AND 
na 1  ret fe. 12  34  1SPW 1C .R AO  I  US. SPL ENG, SPMID.OTWO.OTHOP ,0  THREE ,OTmR£P, 
ICFCuA.CFCuRP 

wfi  1 IE  (  6,1234  )CMv6,  CF  tvtP.  X SPOT  .RTXE.V SPOT  ,RTY£,XU,XIl»XU,T 
RETURN 

2  XSFCT  •  0.0 
YSPCT  •  0.0 
RETURN 

4  continue 
xspct  »  o.o 

YSPCT  •  0.0 
l  -  l 
RETURN 
£  NO 

function  as  xnou i 
ASUC-ST.295?B*ASIN(X) 

RETURN 

ENG 

FUNCT ICK  RONCGI A) 

CIPCRSICN  XI 51 t  « VI 5  1 1 

CAFA  IC-7777T7777777C00CCCC0BJ 

CAT  A  tY *0.0, 0.1, 0.2, 0.3. 0-4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0,1. 1 ,1.2, 1.3,1. 
14,  US.  1.6,1  -  7. 1  .6,  1.9.  2.0,  2. 1, 2. 2, 2. 3, 2. 4. 2. 5, 2. 6, 2, 7, 2. Br2. 9, 3.0, 

2  3, 1,3. 2,3. 3. 3. 4.3. 5, 3. 6,3*7, 3 .8*2. 9,4. 0,4. I, 4. 2,4. 3 ,4, 4,4. 5, 4. 6, 

3  4,7,4.814.9,5.0  ) 

CAT*|X-.C,.5  39  8,.5793f.6We*,65  34*.69l5,.7  25  8,.7  580,.78Bl,,8159, 
l.eMl,  .B64)..8349,  .9032,  .9192,  .9332,  .94  52,.  9554,  .964 1 ,  .  97 1  3  ,  •  9  7  73  , 
2.9621,  .9661  .  *9693,  .9916,  .99  3tJ,  ,99SS,»9465*  .9974»»99fll  ,.9907  ,  •  9  9  90  , 
3.999J».999S»«9997,  .9998*. 9998,  ,9999, .9949,  1 .  ,  L . . I  * , t . , U • l . , l  .  • 1 . , 
41..1..1..I8  I 
I  •  0 

e  -  RCNCti.oi 

lFIC.GT.k5)  GO  TO  1 
I  •  1 

e  -  i.o  -  a 

l  CONTINUE 

CALL  CLINET  5i,S,T,Y,AY) 

IFU.EC.il  AT  -  -AY 
0  •  A  Y  •  A 

RCNCG  •  0  .AND.  C 

RETURN 

ENC 
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SUBROUTINE  CUNEIN,*,»,B.Y  » 

on eos ion  *111.0111 
IFIM  100.100,  101 

101  00  102  *  •  2,N 

I  ■  J 

!F|X-*IJ)I  103.10J.102 

102  CON  T 1  HUE 

10’  It’  f11'1'  *  •  <81  I  1-81  l-U  l/UI  I  l-*(  1-11  > 

too  r  •  o.o 

RETURN 

END 

function  ronoiii 

CUPFCNC  T*J10  I 
EQUIVALENCE ICI 2000), T) 

C*l*  IC.7117T777T7770OOOC000BI 
0*  T 1  I*  .  S77777777777BI 
CA71IK1. 177777777*0081 
If IT.Gf.O.ISCTO* 

2  If  1 1X17.3.0 

7  I  X  -  -  I  X 

8  If { IX. £7.494)6070) 

1 1. 1 1.2. 1 

core* 

J  ix.lx/io 

CCTC2 

*  CONTINUE 

J*l«*2621*7 

IXoJ.ANC.K 

IT.  lx. INC. Nl  I 

8»f L04TI 121/3. N3S97J9E10 

8.  e  .* 

RCNC«8.1NC.f 

RETURN 

End 

subroutine  C2 

CCffCn  ClUIO) 

DlFtNllCMl.no) 

EUvllf>>CEIC«2000l,Tl 
c • *w t no  Inc  costs  focuce 

C  •*  I NPUT  CRT* 

ECU  I V1L  f NCt (Cl 5  4 I.PNINCI 
ECUIVICESCE 1059], SRI  ) 

ECUIviiCSCEtcISn.VNTCV) 

ECUIVlltNCCIMSS),X2K  7 
ECulVll INCf IC.X  1 

fCuIvn.ENCClCIS51.5N  I 
ECUIVllENCE  1C(  50I.0PTNM  ) 

ECUIVllENCE  1CI  MI.BPSIW  ) 

C.FCul PUT  ClTl 

ECUIV*CEnCE  ICI  5  2 ) , VRTE  ) 

ECUIVllENCE  ICI  S3). RHr  1 
ECUIVllENCE  ICI  38C1.R1NCO  ) 

ECulviCENCE  ICI  13)  7  I , R»Il  1 
ECU  I V*c ENCE  ICI  1001, VHXE  ) 

ECUIVllENCE  ICI  101I.VV7E  3 

ECUIVllENCE  ICI  102I.VN2E  7 

C«»1NPUTS  FRC»  OTHER  F0CUIE5 

ECUIVllENCE  ICI  2091, RH  ) 

C*oC»CCUUTS  R I  NO  VEICCITY  COhpONENTS 
IFITiEC.O.IORS  I rv- BPS  1m 
IPI I.EC.O. |VNTEV>yklE 
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IF t ftH.GT.RH*  IGOTOIO 
1ft  CftAKCC  .1.1.  RAH  )  GO  TO  10 
If Ck2k.EC.1969A32)C0TG5 
R2*»1969432 
R  A  KG  T  fi ■ PASCO 
S  CCKflKLE 

I F I RASGTR  *RW  INC.GT.R ANGG1G0T07 
RAKCIR^KANCT ft«RW INC 
C  AL  IGiuSS  I  S‘*,CPS  l  t.  a  PS  UVI 
CALiGAbSstswi.vwre.rfwrEvi 
T  vh*I«-vm£v«CU$CUjfsiwvj 

V*YE»-Vi.f£/*SlN01BP$lWY) 

VWZE  •  0. 

RETURN 

10  VWXl  •  G. 

VrfYE  •  0. 

VwZE  •  0. * 

RETURN 

£KC 

C«*AIR  CATA  HCCULE  G3 
SuBRGuT IK6  Ci 
CCPRCN  C 1  4 3  L 0 ) 

C  •  •  I  N  PL  T  CAT  A 

EQUIVALENCE  lt<0?08J .RHfftO  » 

C«*1NPU1S  fKLP  CTH£R  mOCULES 

ECuIWUtNCE  I  C I  0 100  )  .VWXE  1 
EQUIVALENCE  fCCOlOl 1 »VUYE  1 
EQUIVALENCE  JC(0102).VWf£  1 

CCi'lWAifNCE  (Cl i601 I -VXE  1 
ECUlwtl6SC£  IC(!607)#VY£  J 
EQUIVALENCE  rCMOIlUVIE  » 
cWc  1  »*uSCc  iCI  16ni«Aic  1 
C**1NPVJTS  FRO  MAIN  PROGRAM 
C**StAT£  VARIABLE  OUTPUTS 
C**NCNE 

C**CtHER  CUTPuTS 

EQUIVALENCE  I  C.  (  C200  >  *V.y.WX€  1 
EQUIVALENCE  ( C ( C  20  I  )  .VHWY  f  1 
EQUIVALENCE  IC102C2  l.VMV/E  1 
equivalence  ic i0203i .PCYfiMC* 

EQUIVALENCE  tt(C-204  UVNACH  1 
EQUIVALENCE  CC(C205 » ,CRHO  1 
ECu  I V  A  L  E  f«  C  6  (C<0206».VSCUN01 
E  G  U I VALENCE  !C<C20T WvaIRSPI 
EQUIVALENCE  ICC0209  I  » 

C»«CALCULAT6  present  aliituce 

ft  H  *  -H l £ flQ 

C«»CALCULATE  MISSILE  VELCCITY  WRT  AIR  MASS  IN  EARTH  AXES 
v“wxE  ■  vx6-vwx6 

VHwYf  >  VYE-VWYE 
VKV./E  «  VZE-vwZE 

V  A l R  S  P  ■  SORT f VMWXE«VH*Re*VMWyf«VMWYE*VMWZE*VNV2E) 

C ■" A{ R  CEKSITY,  SPEED  Oft  SOUND  #  CYNAMJ  C  PRESSURE,  AND  MACH 
CRhC»l . 076975 )/( l.*.3325E-04*RHtRH*RH#RH*.023l^E-l2l 
vSCcNC  •  -.  00392««***l\U.3 
POtKNC  *  l0R^0*VAlRSP*VAIRSPl/6*.344 
VWACh  •  VA 1 RSP/VSGLNC 
flETuRK 
EKC 

SuQROUT I NE  CA 

EKO-CP-RUN  CALCULATIONS  SUBROUTINE  GA 

c*»  this  is  a  subroutine#  not  a  hooule. 
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C**  IT  IS  CALLED  £Y  STAGE  3  TO  COMPUTE  MISS  OISTANCE  AND  STOP  IH£ 
C«»  PftCGBir  IF  PANCE  IS  2ERO. 

. . . . . 


COkkCN  C I  A) 101 
C  *  *  I  N  PU  I  CAT  A 
CvtNCNE 

C  *  A  1 N  PU  T  S  FRCP  GTi-fR  POClLES 

CCUlVAiesr.L  IC  1 0  3S  ?1 ,00AMH 

equivalence  fci ose ) »ocamv 
EQUIVALENCE  ICIQ3T1I,  PANCE 
EQUIVALENCE  lt<l&0?l*VKE 
EC‘U  I  VALENCE  LC  1  ibO?  )  .VY£ 
ECLIVA  EnCE  (C816UWVZE 
E  C  U  I V  A  l  E  N  C  E  t  C I  ltbC, I.VTXE 
ECU  I  V  A  )t‘«CE  ICIU6M.VITE 
EQUIVALENCE  IC( «6b?).VTZE 

ecuiva  ence  icmtis 1 »rxe 

ECUIVALENCE  ICIIA19WRVE 
EQUIVALENCE  I  c  <  1A23URIE 
EQUIVALENCE  <Ct  1631 )tRTX£ 
EQUIVALENCE  €  C l  1 6S 5) ►RITE 
ECUIVALENCE  fCIU59l»RlZE 
C** INPUTS  FACP  PAIN  PROGRAM 
ECUIVALENCE  I C I 2C0 0  »  * T 
C««ST*TE  VARIA6LE  CUIPUTS 


C«*NCNE 

C«»OThEP  OUTPUTS 

ECUIVALENCE  f  C  I  202O.LCGNV  ) 

Miss  CISTANCE  PARAPEIERS  ARE  OUTPUT  0 1 REC  TL  V  AND  ARE  NOT  IN  COMMON 

C  •  TEST  fCR  INCREASING  RANGE  AnO  SOLVE  f  ")R  TIME  AT  MHltM  RANGE  IS  *£*0 

5  FORMAT  (1»0,  1  bh  MISS  CI3IANCE"  1PE 1 7 . 8/ 1H0 # LJH  Tl m£  f INAL  * 

C  i?Ll7.6> 

6  PC  P  M  AT  ( 1  NO  »  10X  •  10FXM  6ARlN-lPei7*e*3X,lOHYM  EARTH*iP£ 1 7 . 9 , )X  * 

C  1 0  N  i  M  EAATh-lFUT.fl  I 

T  FCR-AT  I  1N-J.A0*  .10FV  F  L  T  P  A  T  H  •  1  PE  l  7 . 8  •  3  A  T  OHZ  F  L  1  PA  TH-  1  P£  J.7  .  8  I 

C«*  IEST  FOR  INCREASING  range  AND  SOLVE  fjft  lint  AT  WHICH  RANGE  IS  <M 0 

IF  CR  A>.CE  «  G  i  *500 • I  GO  TO  20 
IF  <RANCS-uRAN*>E.LT.O.  1  CO  TO  10 

TCEL  *  -CJJ  <!•  *  UXT-UXE  )  ♦UVYTMUYT-UYE  »  ♦  UV  L  T  •  IU2  T-  Kill  I  *  U  VXE  •  I  UX£- 
C  utT) *IJVYE • l uvf-u YT i ♦  UVZEMUZE-UZTI ) / I UVX T«UVXT-2. *UVXT»UVXC 

C  #UV<  t *UVX£ *UVYT  #UVYT'2.*UVYT«UVYEtUVYE»UVYE*UVZI*UVZI-i.# 

C  UV2  I •uVZE*UVZ€*UVZ£l 

UXTC  •  U'VxT  •TCcL«U*T 
UYTO  «  VVrT •TOEL*UYT 
uzrc  ■  uvzr  »tccl*u£t 
UXfO  ■  UVXE •TC5L»UXE 


U  y  /  C  •  LVre •TDEtRUYE 
U2M0  »  uVl E *TD5L ♦UEE 
PCX  -  UX^C-UXTQ 
RCV  •  U>MC*UYTO 
RC/  •  uiMQ-UZTO 

R»MSS  ■  SORT  |ROx*RCX*ROYtRCYF.ROZ*RDZ) 
TiERC  •  TC£L*UT 
USTHT  »  SlNCieCAMY) 

UCTHT  -  COSCICOAMv) 

USPSl  •  S INOI 8GAMH ) 

UCPSI  -  CCSC  <  0CAMM 
UCll  •  UCTHT.UCPSI 
L/C  1  2*  -UWmI»uCP$I 
UC13*  USPSl 
UC2l«  USTnf 
UC2  2*  UCfHT 
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CCH*  -UCTMT  *US  PS  I 
UC32-  USTmT«uSPSI 
UC33-  UCPS1 

P«FP  -  vC l \ ■  RC.x  «uC  l?«ROYtUC13*ROZ 
Pvpp  -  uC2  t  *ftOX  «UC  ?2»RCY 
AJFP  -  lC  31  «AOx  «L‘C  32*RCT*uC33*ftOI 
IT£  (  6.6  )  PMISS.  UCRO 
N  P  l  T  £  (6*6)  RCX,  PCY,  ROZ 
N*[TE  (6,7)  RYfP.ftZfP 
LCC.W-2 
REToRh 
10  Lit  »  1 

VJ  K  E  ■  R*E 
tve  •  RYE 

uzt  •  m ze 
u*  t  •  r r x e 
ox  r  ■  a rvc 
uzr  •  ptze 

UYXf  •  VXE 
LVfE  ■  VY£ 

LVZ£  ■  VZE 

u v <t  ■  vrx; 
uvyt  •  ^ t r 6 
uvzr  -  vwe 
20  L»*Kce  •  range 

If  <«Z6  .GT.  100.)  CCONV  *  2 

RETURN 

ENC 

C«*CCC»C:naI 6  CONVERSION  H0DU16 
su'PctriKE  C6 
CCPNtN  C  *.«.  3  10  ) 

c 

CMlNPUTS  FRCP  C  f  H  E  R  PCCUltS 

EtuIVALEKC**  ’  *(0200  UVHMXg  I 
EGUIVALENC-  :i.<  020  |  UVN»Y€  ) 
EQUIVALENCE  IC(C202»tVHW£E  ) 
EQUIVALENCE  ICIC20  7)rVAI«SP> 
EQUIVALENCE  ICt  160  3  J.VXE  ) 

EQUIVALENCE  IC116Q7I.VYE  ) 

EQUIVALENCE  ICI 16  1 1 ) . VZ  E  ) 

ECwIVALENC:  (  C  (  l  6  1  6  1  t  R  X  F  ) 

EQUIVALENCE  (Cl  1  *  I  S  >*.Ry£  J 

EQUIVALENCE  IC<  I-4  2  2  ),*ZE  ) 

EQUIVALENCE  (C  f 1 6  3  6  ) »  RC  lX  ) 
EQUIVALENCE  ICI 16351, RDELY  ) 
ECjIvAlC.sCE  (Cl  16)7  ),RC£LZ  1 

equivalence  ici  1663). vcxa  ) 

EwulvAiENuf  (Cl  1 6  6  <■  )  i  V  C  Y  ft  ) 
EQUIVALENCE  ICIl666)*VGZ8  ) 
EQUIVALENCE  ICI  16.6Et.RX0  I 

EQUIVALENCE  1C  I  1666  )  #RyQ  ) 

£CUI VALENCE  (Cl  1 6  TO  I  . R/0  1 

E  C  l  I  V  A  L  t  S  C  £  ICIIT03UCFAH  ) 
EQUIVALENCE  1C  i  t  70  7  WCFA12  ) 
EQUIVALENCE  (C  I71IUCFAU  1 
EQUIVALENCE  (Cl  ITli  I.Cf A2I  ) 
EQUIVALENCE  ( C I  i 7  l  9  ) , C F A  2 2  t 
EQUIVALENCE  ICI 17231, CFA23  ) 
EQUIVALENCE  1C (  1 72 7  I.CFA31  ) 
EQUIVALENCE  KCITM  l*Cf A32  t 
EQUIVALENCE  ICI 1736 UCFA33  ) 
ECuIvAiENCE  (Cl  1  7  39  J ,WP  I 
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€CUJVHENCE 
ECU! V*l EVE 
ECO  ivnescE 
ECU  I v  *1  ESC  £ 
ECU1 VAifsCt 
ECUlV>LEV£ 


iciuui.ng  } 

I C I  1 74  7  I  »WR  1 
lC< 1751  t»C*AO  I 

(CI200C1,T  J 

lC<266*),0ER  1 

(Cl  350* UOPFN*  I 


c 

c«»QfHE»  curruTs 
ECU!  VAiEVE 
ICUIVAI EVE 
ECJiVAl  FVE 
ECU! W ?s:f 

equivalence 

EQUIVALENCE 
ECulv&i CSCE 
EQUIVALENCE 
ECU! VALENCE 
EQUIVALENCE 

equivalence 

EOUlvAl  EVE 
EQUIVALENCE 
ecu! vai ENCE 
equivalence 
EQUIVALENCE 
EQUIVALENCE 

Etu I VAI escc 

ECuIvai  f.SCE 
f Sul VAL  ESCE 
ECU! VAI (NCE 
ECU  I  VALENCE 
EQUIVALENCE 

squivalq-ncc 

EQUIVALENCE 
EQUIVALENCE 
EQUIVALENCE 
ECul VALENCE 
EQUIVALENCE 
ECul VALENCE 


IClOBSOl.BTH  I 
ICI035  l  1 ,0PS  I  I 
1 C  (  0  35  ? 1 I  | 
1C  I  0  352  J  »C ThTO  ) 
10  054 1, OPS  10  I 
1C  10355  I.e^lD  1 
IOC35  61.VT0TC  » 

1C  1C  35  7  J .BCan*  ) 
(C I  0 15fl  I  ►ftCAHV  ) 
1C 103*0  I , V«wU  » 
1C  I  06  I  I  »V^MV  I 
1C  10)6?  I  *.V*4W  I 
10063 1.B1W  I 
1C  i  0  36*  Vl!pSLV  I 
<C(0)65),6LAMy  > 

1C  I  C  366  )  >01  A*H  ) 
ICI0367  U9ALPHA I 
(CIO  368  I, 0ALPHY  I 
IC(O369  :t0ALPHP  » 

IC(C370  W0^MP  I 
fC ( C3 J  I  J ,p anCE  ) 
100372), M0A  ) 
10073  )fRTOA  1 
JCI0374  URfCA  ) 
100)75), OLCS&P  ) 
(00)76), BLOSBY  I 

(C(C377l,CALPO  J 
100)78), fiALYO  ) 
ICI0379I.0AIPPC I 
Id  380I.R»NG0  ) 


C»*C«UUIAUCN  OP  HEaCInC.  PITCH.  ROLL  EULER  ANGLES  IN  0£G*E£S 
BPhJ  *  ATAuCICf A?),  CPA  3)) 

BThT  •  A(  ANCI-CFAl  JtSURT  ICFAlMCf  All4CPAU«CEAUM 
BPS  I  •  At  AUCtCf  A12.CEAU) 


c 

iPicasci ®tht i.cc.o.oi  co  to  s 
bpsic  •  uu«siNCie?HiiMR*coso(0PHm/cosoteTHU 
5  CCAMNL'E 

e pm i c  •  wp*opsic*s  jnoibtht  i 
BTHTO  -  WCVCSOICPM  l-MR*S  indibpmii 


c 

C»*C  ALCUL  ATI  C.N  CP  TOTAL  VELOCITY 

VTCTE  •  SCRT1vxE*vxE*VYE«VYE*.V2E»Y?E| 


c 

RASGO  •  5  CRT  1  I ftx  E -RXC) •• 2  ♦  (ftrf-RY0)**2  *  IRXE-R J0»4*2l 


c 

C*»TRAASFCRN  piss  he  lcs  mo*  EARTH  TO  BODY  axes 

PtXBA  •  PrFLX«CFAll  •  P0ELY-CPAI2  r  ftoCL4*CFAL3 
RT  8  A  •  RCELX*CPA?1  ♦  RCELYVFA22  ♦  ROEL2»CFA?3 
R?  0  A  -  RCELXVEAH  ♦  R  CEL Y  *CEA  32  ♦  ROEU*CFA33 


C 

C««*lSSll£-TGf  U  IN  BOEV  AXES 

6LCS0P  *  AT  *  NO! -RZ  e  A^CRTI  RXBA»RX8A*r,YeA*RYBAn 
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BLC38Y  •MANCr  AYBA.mce*  ) 

c 

OVP  1  •  VX£«BOflK*VYE*«Cfcr 
OVP2  ■  SC6LX*«»CEl.X  ♦RCELT»RCElY 
UVP3  •  VIE«RCfrl2 
Uv P A  -  SCRTIUVP2) 

RANGE  •  SCanuvP2*»CK2»«2J 

{♦•VERTICAL  A.SC  MCRI/CNTPL  LINE  OF  SICHT  ANGLES  DEARTH  AXES) 

c 

8L«^H  •  AT ANC(-RDfcL Y.ROfLK 1 
ei«*/  •  AT  ASC( *RGE  L  Z *UVP4 ) 

c 

C*»VERTICAL  a.\c  horizontal  prupurtional  navigation  angles 

IF  I T-LT  .CEP  I  GOT  C30 

IF  I  RANGE. EO.C.O 1  CC  TO  2 

VXP^IGVPl*JVP3)>RAKCE 

2  CCMINl-E 

|F IUVP4.EC. O.O)  CO  TO  l 
VY?  *  I  VYf  «RCELk-VXE«RCELY  I/UVP4 
V/P  •  IV/E *UVP2-ftCEL2*UVPI  1/IRANC6*UVPA I 
i  COST  I. SUE 

6IHLV  *  AT ASOl VZP , VXP) 

8  P  S  L  V  -  Af AKO( VrP.VXPI 

c 

CCAHV  •  ATANCC -V2€ » SCRTI VX£«V*E *VYE»VYE ) I 
BGAPM  •  it  ASCI  V Yc i VXfc ) 

c 

C*«VElOC I TY  kRF  A | R  |A  6  CCY  AXES 

WFWU  •  CFA1  1  «VHtaHE *CFAI?»VHWyE ♦CFAl 3*v*NZE 
VNfcV  ■  Cf  A2  1«VMhx6  •f.FA?2*VNWYE  *CFA2  J«VNW26 
V^H  •  Cf  A3  1  •VNMKg  •CPA32»VM*n  E  «CFA  3  S’VMHZE 

c 

C**VCMlC«l  ASC  hOKUCNTal  angles  of  ATTACK 
8ALPHA  •  AT ASC( VM^k. VNrfU) 

8  A  C  P  H  Y  •  AT  A»,0(  VHUV.VHWU) 

c 

use  Y  VP«'J*"2 

|FluSw.EC.0*O.Af.C.VPMW.E0,0«0)  CO  TO  } 

8  A  L  PC  •  |VM»i;«Vti3  -  VM»W#  VlJXP  )  /  (  USQ  ♦  VHNW#  *2  )  »C*AO 

3  CONTINUE 

IF IoSC,£C.0#0»anC. VPwy  « tC . C  .0 )  CO  TO  A 
EAlYC  ■  lv«MU*VOYO  -  vn«V*V0*3 >/l US0*VNWV*«2UCRA0 
A  C  C  s  T  I  N  i;  C 

0 A l f  PC  *  0. 

If  TCAlPhP.OT.O.  )  EALPPC-  <  0AL  PmA*  8ALP0  A  6ALPHY*8ALY0I /3ALPMP 
C 

C  •  •  A  l  PIJ  PRINf  A.NC  PN1  PAINE  (WIND  TLNNEL  AXES) 

If  (  I  P  *  i.P  h  z-  &  al  PhY  )  . E  C .0  *  I  CO  TO  JO 
0Ph  I P  »  ATASOI VHWV, V^WH  I 
30  KALr^P-SC3T(BALPHA*«2*fALPbY»»2) 

IF  I  APS  I CAUPFP) .or* 20. ldALPH»«20. 

fit  Turn 

ESC 

C«  *  SEEKER  ANO  PLATFQRN  I  Nil  MCCUlE 
SuDRCuMSf  S 1 1 
C  CP  v  CN  CU310I 
CIPENj  ICS  I  HI  100) 

CvUlVALESCE  *C  <  *tS)»UI  ) 
ecu! VAt 6SC6  ICf  4SIUC«CE  1 
ECUlVALfACE  !C<  <.60. GUIDE  ) 

EGuIvAttNLE  <  C  (  A6l)»SA>*P  I 
EQUIVALENCE  ICI25M1#N  ) 
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€CU1VU€HCE  1C  <  2562  ) »  1  PL  i 
SZ-O.# 

SY-O. 

UT  •  0. 

CUICE-i. 

CAGE  *0 • 

SAPp  •  0* 

Cl  452)  •  o, 

!pllN  1  •  4?* 

IPLCNM  J  •  428 

h«N  *2 

RETURN 

INC 

MC£fl  PLATEOAN  AND  TSAC*£A  NODULE 


3ua«GU(IK(  51 
CONPCN  CI4310I 

c 

C»«lNPur  data 


c 

c» 


EQUIVALENCE  K  (  054?  )  .StGB  I  S  I 
EQUIVALENCE  ICt  443),QPTkA  J 
ECulvALE'.CE  ICI  444),QPrBKLI 
EQUIVALENCE  ACC  4451, UT  I 
EQUIVALENCE  (Cl  4461. COT  | 

EQUIVALENCE  ICI  4471, CKCB  I 
EQUIVALENCE  (Cl  44  6  ) , CF  CV7  1 
EQUIVALENCE  ICI  44S),CF0VY  I 
EQUIVALENCE  ICI  16771, INSLI 
EQUIVALENCE  ICI  4501, CSM  » 

equivalence  ici  4iii, case  j 

equivalence  1C (045V  I ,oo«EAK 1 
Equivalence  1C i 045 7 I, R BKLOK 3 
EQUIVALENCE  1C  I C454  1  ,OECC£  ) 

EQUIVALENCE ICI 405|, AVI 
EQUIVALEMCEICI  404  I, Ri  I 
EQUIVALENCE  ICI0455 ) , v  E  A  SHX I 

ECulVAl  ENCEICI  4701  .SI.-ICL  ) 

ECU!  VAL  C.NCE  I C I  4  55  I  ,  7L  ASA  3 »  I  Cl  4671,  YLASft  I 
ECUI valence ICI  4  63  I, CffoCSI 
EQUIVALENCEICI  4.591  .COTAAOI 
EQUIVALENCE  I C  I  C4  55 ),CK  SNA  | 

EQUIVALENCE IC(47)),VLA7AP| 
EQUIVACENCE(CI477I,CVSMI 
EQUIVALENCE  1C  I  045  I ),CRCS7P I 
EQUIVALENCE  1C  I C45e > , CRO SP I  I 
EQUIVALENCE  ICI  4601, CU ICE  ) 

EQUIVALENCE  ICI  4611, SANP  1 
EQUIVALENCE  ICI  464  I .CCAHVS I 
EQUIVALENCE  ICI  4  6  5  ' , C  C  AMM  5  I 
EQUIVALENCE  ICI Oe55  ) , CJ  t 
EQUIVALENCE  (CIC056I,GT  1 
EQUIVALENCE  ICI 5504), 0PTN4  I 


•INPUTS  FRO  CTNER  NCCUlES 

EQUIVALENCE  I C I  0  3  7  1 ) , RANGE  I 
EQUIVALENCE  1C  I  0  JTP  l,RX8A  1 
EQUIVALENCE  IC (C J 7 3  I ,RYBA  I 
EQUIVALENCE  IC(0374),AiBA  I 
EQUIVALENCE  ICI 16151, AXE  1 
EQUIVALENCE  ICI 1449), ATE  | 
EQUIVALENCE  ICI 16771, RAC  J 
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ECUtVAlfNCE  ICC  1714UW*  ) 

ECUlVAtE'.CC  ICll?AO»»WOO  ) 

equivalence  ic(i?43»^wg  > 

EQUIVALENCE  ICIl?AA),WftC  » 

EQUIVALENCE  ICI17A7  )rUR  ) 

EQUIVALENCE  I  C  <  2000  )  ►  T  I 

c 

0*$  TATS  VAfi  I  iPie  CuTPuTS 

EQUIVALENCE  1C  I  0*24  WBThTGO  I 
EQUIVALENCE  IC(  C*2  Hrd  TMG  ) 

EQUIVALENCE  1C  <  0*28  I,  BPS  I  iC  I 
EQUIVALENCE  1C  t  QN  3  n  ,  BPS  10  1 

c 

C«»OTMER  CUI0UTS 

EQUIVALENCE  t  C  (  0  SO  3  > .  c  Z  ) 

EQUIVALENCE  I C I GSO  7 ) #  E  Y  ) 

EC-IVAl-ENCc  !C(Cs32»,RXG  I 

€QaJ  IvallnCE  ICI0sj3»tRTG  I 

EQUIVALENCE  IC!0s3s>*RZG  ) 

ECU! VALcSCf  1C! os 35 J ♦ BEPSZ  l 

fCuI/AitNCE  1C I0A1* ireEPSY  ) 

EQUIVALENCE  !Ct0O7),WZ  \ 

EQUIVALENCE  IC<OS39>fWY  \ 

EQUIVALENCE  !ClOs5t*),6GCEFL  I 

EQUIVALENCE  ICt  *b2>rUl  ) 

EQUIVALENCE  I C (  N63I.UY  » 

C 

C»*0! REC  T I  ON  COSINES  FOR  BOOT  10  PLATFORM  TRANSFORMATION 
|ftl,GT«0.)CCT030 
$L  T  $»0  . 

$2-0* 

Sv*0. 

>0  CCNTINUE 

ue 3 1  •  siNOieiMTG) 

UB33  •  CCSCICTHTG) 

U012  •  S  I  NO  J  l»P  S  IG) 

uez?  •  ccsoiopsigi 

UBli  •  LB22-U933 
Uen  •  -U03MU92Z 
UB21  •  *U0  3  3  *U  0  12 
UB23  *  UBJ1«U012 
UB  3  2  •  0. 

C 

C*A  CALCULATE  total  DEFLECTION  QF  GIMBALS 
BGCEft-SCRT  <  BTHTG**2FePSIC««2) 

C 

CA»TRANSFCRM  LCS  FROM  fiOCY  TO  GIMBAL  AXES 
RAG  •  U211 •RxQA*UB12-RyOA*uB13«R2BA 
RYG  *  ue2L'RXuA.u6;2*RYaA*Lfi23«ft2aA 
R2G  *  UC Jl»RxCA»ue23*R28A 
PCl-ATANC(-ZHSL.RANCE) 

c 

C»»CHECN  FCR  MISSIVE  AT  SEEKER  BREAK-LOCK  RANGE 
IF  (RANCE.GT.RCNLGKl  GO  TQ  AO 

If  IORSAk.lE.O*)  CC  TO  35 

IF  ICIChCK.GT.O.  \  CO  TO  4<J 

C 

C«»L!f»E  CF  SIGFT  RATES  AFTER  9R  E AK-LCCK 
LEP5  -  I  T-UT  IM£  )*w(PSMX 
IF  luEPS.GT . BE 0C6)  UFPS  B  EOGE 
BE  PS  2  •  U6P52*UEPS 
BEPSy  -  U  C  P S Y • u E PS 
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OIChek  *  1. 
cc  to  50 
c 

t,,^‘ ”JLU»tlC«  op  BREAK  10CK  yMUHU 

UEPS?  »  66P3Z 
l/EP$v  -  OtPSV 
CQPfAK  *  1* 

4,3  i6»2D0)  T  «  R AfcGf 

c  100  FC*"'  '1OH0BREAKLCCK  PAS  OCCURRED.  TIRE.,Fe.9.»H.  RANGE. .PW.A, 

C**tOS  ERRORS  IN  PlAtfORN  COORDINATES 
*0  BEPSE  -  AIANCI -AEC.fUCl 
eEPJy  .  ATANOI  RYC.RKGI 
CO  10  50 
AO  C1CPE*  •  -I. 

50  CCMIMJC 

ip  itnn  .gt*  o. )  go  ro  «o 

SE  .  CKSKR*bfPS2 
it  -  CKSKDtOfPSY 
CAGE  .  1. 

GUICf  .  I. 

CO  la  92 

•0  IP  ir  .LE.  Uf I  CO  to  92 
or  »  at  »  cct 
SEfiO. 

A  A  A  «0  • 

C»l(CAUSSISEf , AAA, AAV) 

IflCIlOJI.GTkVLAERPIGOTOfl 

s$r»o, 

SjE’ECAj^'jINuIFCI)  /RANGE 
SSY*VLASft/BANCE 
*0  PEG>REC>SSE.RAnCE 
p<C«RkG 

BEPSE •ATASOI-REC.PKCI 
PYC-arC«5sr*«ANCf 
8E  PSY.A 1 A NO I PVC, PXC I 

c*l.A^e“*cocD’',C’C,f0C‘,**'#B’CC,0D'eE',w,sei>ST> 

|F<Cf.ll.l.e-5)CF*lO' 

I CC -2 e ) *CKSKR/CP 
3Z*  tCOAA  l«CK$KR/Of 
CAce-i. 
ccrc?2 
91  Si-O. 

SY*0. 

92  Ci  •  Si 
UY  •  $Y 

c  IP«CF.Gr.O. )CICHECR>!« 

CPItCH  PPIJCRA-AINC  ANC  SEEKER  CAIN  SOItCNlNC 
If  <  GU I Ct.Ct .01001020 
ip  isipp  .or.  o.i  oo  to  19 
If  ICAOe  ,L£.  0.1  CO  10  21 
Of E  •  VI 
SAPP  .  1. 

Utl  ’£S*  £  ICN 1 1.  a  UEZ » I  Co  to  21 

20  0AFCC$»5?.6».00ft67/C0r/C0T 
ll*  UZ/CK$fcft«OAf ocs 
Cl'»UY/CK$K««CAfOCS 
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Ot-82 
ur-ir 
CO  TO  21 

21  £2  •  CCAPVS/Gl 
EY  *  CGAMmS/GV 
02  t  CS44U2 

OY  a  CSk"UY 

22  CONTINUE 

c 

CA4NCN-UNEAP  »0  I  c  T  [  COUPLING  OF  GIMBALS 

02«*SICMO»CSP2,8PSICD) 

UymS  ICMCRCStP.ETMCO  ) 

c 

c*»Hissue  eccr  rates  in  gimbal  axes 
NY  «  02  2  1  *WP  .UB2?4l,(;  ,u823aIiR 
Ni  i  ofi31aNp*ud32»NO*OQ33«Nft 
C 

0ilPC  •  (02  ♦  S2G8 IS  ♦  U2K 1 

setAc  •  (oy  a  sYcais  -  lyx)/ub22 

c 

C*4GIMBiL  ANGIE  DERIVATIVES 

If  ICAG8  .0 1.  0.1  CO  TO  89 
SINTCC  •  I8ALPQ  -  NYI/OB22 

ePS ICC  •  I 8ETA0  -  Nil 

C  PIG  OPTION 

c  et(-rGc--CKSK34arwTG 

C  e“SlCC--CK5X34BPSIG 
RE  TORN 

btmico  4  o. 

epsiGc  ■  o. 

RETORN 
f  Nr. 

C«*  TIGER  AOTORILOT  IN  IT  I  Al  HAT  ION  NODULE 
CT'f'T'ICa  FREQUENCY  PCTEL444... 

SuePColiNE  cn 

CCE-CN  C (43 10 ) 

DIpENSICN  IPKIOOI 
EQUIVALENCE  TO  835  I  ,E2S  ) 
equivalence  ici  idi.Evs  i 

EQUIVALENCE  JC!  88  3  1  . E25S  I 
EQUIVALENCE  ICI  8871, EVSS  1 
EQUIVALENCE  (Cl  »64),CCANY$) 
EQUIVALENCE  ICI  4S5I.CCAMMSI 
ECl.IYAlt.NCE  ICI 35041, 0PTN4  I 
ECutYALENCE  ICI  256D.N  ) 

ECOIVAIENCE  ICI2562I.1PL  I 

C 

NN SOP  •  H 
IPlIN)  4  800 
IPLINaII  .  820 
I  FI  I  N. 2  I  •  8  24 
IfllNAll  4  828 
I  PL ( NAN  I  •  832 
I  F v  I N  a 5  I  4  834 
IFLINA6I  4  $40 
IFLINaTI  4  880 

IPlINAj)  4  884 
N  •  Na9 

E2S  •  0. 

EYS  4  0. 

GC  TO  22 
21  EJS  •  CGAPVS 
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€YS  •  CCA*nS 
22  E 2 S $  •  E/S 
CTSS  ■  EYS 
Cl  803)  •  0. 

Cl  8231  *  0. 

Cl  827)  *  0. 

Cl  831)  •  0. 

Cl  839)  »  0. 

RETURN 

ENO 

C»*  TIGER  AUTCPllOf  POCULE 
C»*«***LCw  fRECUEhCY  POCtL**«**« 

SUBROUTINE  Cl 
CCH^Ch  C (43 10) 

oiPeNSicN  eccircM ),vari  ion 
c 

C«*lNPUT  CAM. 

equivalence  ic(O8S0),hli«o  ) 

EQUIVALENCE  ICIO05 U ,*L IMS  ) 

EQUIVALENCE  ICI  052)tO6IAS  1 
ECulVALENCE  ICI  85}),ftB|A$  I 
equivalence  (C[08S3),02  ) 

EQUIVALENCE  I C l 085  C ) , GY  ) 

C  C(  85  7  )  T  )-Au  Cl  85  Q  |  AR£  uS60  OY  ECNTHl ( l ) 

EQUIVALENCE  TCI0843), TAU2  ) 

EQUIVALENCE  (ClC064  ),MUY  ) 

EQUIVALENCE  ICI  065  )rTCVl  ) 

ECUIVAI.ENCE  ICl  86  8  )  *  TC  Y  2  > 

ECuIvALENCE  ICl  8  7  7  ),  TAUL  ) 

c 

C#* I NPUT  S  ERCP  Cl>ER  POCULES 

EQUIVALENCE  t  C I  0  IS  7 ) , rpn |  | 

EQUIVALENCE  ICl 0)55 ), 0Hh 10  l 
EQUIVALENCE  (CIQ40)),E/  ) 

EQUIVALENCE  IC104Q7),EY  I 
EQUIVALENCE  IC(0848).N)-SUM  ) 

EQUIVALENCE  1C l  1739  )»WP  l 

equivalence  icii740),wcc  ) 

EQUIVALENCE  (CIIT43UWC  I 
EQUIVALENCE  ICllTAAl.WRO  I 
EQUIVALENCE  ICl  174  7  1  ,  WR  ) 

EQUIVALENCE  ICl  1751)  *.t  A  A  0  ) 

c 

C«*lNPUTS  ER  L  H  PAIN  PROGRAM 

EQUIVALENCE  ICl  2000  1  r  T  ) 

EQUIVALENCE  1C l 2965 ) » VAR  ) 

ECUIVALENCE  1C  I  2664  I  *.GEA  I 

C 

C»*  STATE  VARIABLE  OUTPUTS 

equivalence  ici  80C)^bpmi5C) 

EQUIVALENCE  ICI  8C  3  )  »>G?r  I S  \ 

EQUIVALENCE  ICI  8  2C I #E  SUHQC I 
EQUIVALENCE  ICl  8  ?  3  I »  E  SUMQ  I 
EQUIVALENCE  ICl  824)vESUMEC) 

EQUIVALENCE  ICI  827), ESUME  I 
EQUIVALENCE  ICI  628  )  »E7  SCO  I 
EQUIVALENCE  ICI  83IWE/SP  ) 

EQUIVALENCE  ICI  832), E/SO  I 
EQUIVALENCE  ICI  835), f/$  I 
EQUIVALENCE  ICl  836UEYSCO  I 
ECUIVALENCE  ICI  «39>,EYSP  ) 

EQUIVALENCE  ICI  840), EYED  I 
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ECutVALENtE  TCI  043). EVS  I 
FCUl VALENCE  1C!  08O.E7SSO  I 
EQUIVALENCE  ICI  6831, EZSS  ) 

EQUIVALENCE  1C!  834). EY&SO  ) 

EQUIVALENCE  IC(  68?),EvSS  ) 

C 

c*40ufputs 

EQUIVALENCE  (C I  85?),BCELTCI 
C 

C«*OThE A  CUIPUTS 

€CUv£LENCE!C«Ce7t)*fZAft) 

Ca»platfcrh  p a r l s  in  inertial  space 

SCUIVALENCE  It  IOfle>@  J.EVRR  ) 

EQUIVALENCE  ICI  e6<?),uGC  ) 

EQu  W&i C\Cs  ICC  0  70  ), wRQ  ) 

€CuI/ALt?.Ce  I C I  C  S  7  3  )  .  EC  CCR  ) 

EQUIVALENCE  IC10e7i.»#6vNCA  I 
EQUIVALENCE  (ClOG7S)tbCELPC) 

c 

c 

C  »*Gu I C akCE  SIGNAL  SHAPING 
f  z  sc  •  ezsp 

E V SO  »  EVSP 

EZSCC  -  TAUZ»ITAUZ«<GZ-€Z  -  SZ$l  •  2.«  F  Z  $0 ) 

E  V  SCC  »  TAUYaI TAUY •IGY«£Y  •  £YS)  -  2.«£yS0I 
EZSSC  •  TAJZMEISU/fAUL  «f  EZS  -  EZSS) 

EYSSO  *  T  Aijy.(  EYSC/TAUl  •  EfS  -  EYSS) 

C 

C**GRAVITY  AhC  PATE  BIAS 
WQC  ■  EZSS  ♦  081AS 
H8C  ■  6YSS  A  R0 I  AS 

c 

C*»acP*  rate  shaping  and  gyro  dynamics 

*CS  ■  HQ 
MRS  •  MR 

If  (AflStwCSl  .GT,  30.)  MGS  *  SIGN! 30. ,  MQ$| 

I  f  U3SINRS1  .GT.  20.1  MRS  v  S1GNC3G.,  HRS) 

C 

C»*SUHHATICN  Cf  Rate  camping  ano  cuicance  signals  ano  their  DERIVATIVES 
EZ3R*  HCS-WCC 
EyR*  -  MRS  -  MAC 

c 

UKR  •  .85 

IFIT.LT.TCV2)UKA«4.25 
IFIT.LT.  TCY  I  HJK R » 0  • 

EStHCC  a  UKR«<EZRR  -  E Y Aft ) 

FSuMEC  •  UK R  •  I E Z RR  *  EYftR) 

c 

C«*TCTAL  CUICANCE  SICNAL  SHAPING  ANO  tlMtTlNG 

EQCCk  »  <FSuM0u/6.  ♦  ESUH-O) 

e v ncr  •  <  esuKEOye.  *  esukf) 

IF  rAeSlFCCCR)  .GT.  HUMC)  FOOCR  •  SfGNIHllMO,  EOOCft) 

IF  UeSIEVNCR)  .CT.  HLIME)  EvnCA  •  SIGnIHLImE.  CVNCR) 

c 

C •  • RCL L  SIGNAL  SHAPING 
URP  %  .33 

IF|T.LT.TCy2)UKP»1.65 
IHT.U.TCrl  I-jKPaO. 

UPHIS  -  UNPM8PHIC/12.  a  8PHII 
BPl  ISC  ■  t6.MUPH!S  -  OPMS) 

IF  I  DBSI8PMS  I.GT.2.  )  BPh  I  S  •  S  IGN I  2.  *BPMJ  S  I 
BOELPC  •  -6PHIS 


c 

C'MuTOPIlCr  OUTPUT  CURRENTS  TO  EACH  ACTUATOR  {FROM  SUMMATION  AMS  I 
BOElTCI l)  •  EOCCR  -  BOUPC 

eoftrcm  ■  cvncr  -  bceupc 
Rctircm  •  egocr  *  eoeipc 

BUUCI4)  -  EVNCR  «  DOEIPC 

RETURN 

ENC 

C«*  TICER  SIMPLIFIED  ACTUATOR  NOOEL. 

C«****UCW  FREQUENCY  MQCEl  •••••• 

c 

SUBROUTINE  C A 

c 

COMMON  C 1*3  10) 

C  OIMENSICN  BUEUTOI  * )  ♦ODEt  M A)*BOElTCIAJ * VARI 101) 

DIMENSION  BCELTIAl.nCCuTClA) 

c 

C*A I NPUT  CAT* 

EQUIVALENCE  TCI  1121 ) ,BDMAX) 

EQUIVALENCE  ICIllAO  WCEUPBI 
EQUIVALENCE  1CI  11A1 )»OELTOei 
EQUIVALENCE  I  C  t  1 1*2  J ♦  OELTR0  I 

c 

C*%INPUTS  FRCP  C  T  E'ER  MCCUCES 

EQUIVALENCE  iCl  IIUWNFOEITI 
EQUIVALENCE  I C I  8S  ?  >  *  BOELTC ) 

EQUIVALENCE  ICI  IU7).BSURF|) 

EQUIVALENCE  ICT 1110 ) rBSURFil 
EQUIVALENCE  TCI 1 1 19 ) »  8SURF  3 ) 

EQUIVALENCE  4CU120  >.B$URF*) 

C 

C  •  •  F  *.  *  p  PFPLECTION  RfAS 

eceiTcii)  •  eceiTcm  -  ceitpb  s  oeutqb  -  oeltrb 

BCELTCI?)  •  BCtLTC I  2  I  -  DEL  TPB  a  DEUTQB  a  OELTRB 

eo€irci3»  •  eOEiTcm  ♦  oeltpb  a  oeutqb  ~  oeltrb 

BGElTCUl  •  BOEUCIAI  ♦  celtpb  a  oelvqb  *  OELTRB 

c 

COACTUATCR  CVNAMICS 
OC  30  1*1,4 

BOEUTI I  |  «  BCEuTCI I » 

C 

COSURFACE  POSITION  LIMITER 

IF lABSIBCELT II) ) .LT «B0MAX)GOTO30 
BCE  IT  I  I)* SIGN!  BCM.AX  *  BD6LT  (  III 
30  CONTINUE 

C 

8SURFL  *  B06LTI1I 
0SUPF2  *  B  D  E  l  T  <  2  I 
8SUPF3  •  6DELTI31 
0SUPFR  *  BDfcLTIA) 

c 

Cl  110  3  3  *  eCELTIl) 

CI11C7)  «  8CELTI2) 

emu i  •  eoELTiji 

C(lllS)  •  8CELTUI 

RETURN 

ENC 

SUBROUTINE  A 1 
C 

COMMON  C 14310) 

c 

C**T ABLE  ICOKUP  FOR  BGCY  FORCE  COBFFICIENT* 
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COPHCN  /KCKO  /NCXO  /Cx  0  A  RG /  CIOA  /CXOFUN/  CXOF 
'  /NCXCP  /  NCX  /CXAAC  /CXA  /CXfUX  /Q  XF 

*  /NC2  /NCN  /C 2ARG  /CNA  /C2FON  /CNF 

*  /.sncr/NcrcN  /ccarg/cocna  /rc/f un/cocnf 

^  "  /NCV2  /NCY2  /CV2AAG/  Ct2A  /CV2FUN/  C»2F 

OMABLE  CCOlLP  f  C«  CLOY  HOP  ENT  CU£Ff|C!CN1S 


1  /NCL2/MC12  /CE2AAG/  CL2» 

*  /  NC  L  3/  SCL3  /CL  313G/  C  L  31 

*  /NCH/  NCR  /CHAAG  /  CHA 

B  /NCCH/  NCCCN  /CCH1HC.  C0CH1 

*  /  NQ  N  2/  NCSA  /CN2AAG/  CN21 

•/NCLP/  NCLP  /ClPTPG/  CLP1  /CL 

•/NC*G/  ncrq  /chcipg/  Choi  /Ch 


‘  /CL2FUN/  CL2F 
31  /CL  3 P LN/  CL3F 
/CMFLN  /  CHF 
chi  /OCHFLN/  COChF 

2*  /C'NTiV  CN2F 

/CLPFOM/  CL  PF 
/CHOFLN/  CHOP 


C*priSL6  LCCXUP  FOR  SURFACE  COEFFICIENTS 


/HC2C/  F.C/C 
‘/NCIC/  LUC 
/NCHC/  NCHC 


/C2CA1G/  C201 
/ClFAaG/  CLOl 
/CHC4SC/  CHOI 


/  C  2  0  F  L.N  /  C2DF 
/CLOFI.N/  CLDF 
/CHOFLN/  ChDF 


c*«  input  nn 

EQUIVALENCE  tC I  1 232  I  ,X INTER ) 
EQUIVALENCE  I  C I  12/3  I  >  CXERR  ] 
EQUIVALENCE  1CI  1261 1  ,C2  ERR  I 
E  QU 1 V 1 L  t NL  t  I C I  126/ I , CV ERR  I 

EQUIVALENCE  ICl  12631, CLERR  I 
ECulVILENCE  (Cl  1261 1 . Ch£RR  I 
EQUIVALENCE  IC1 12651.CNEHR  I 

C 

C**1NPUT5  FRCP  CTeER  HOCLLES 

EQUIVALENCE  IC( C201 1 ,VH1CH  1 

eclIvacencl  ic 1036 ti , eicPHii 

EQUIVALENCE  (LI  20  2  0  t.LCOHVt 
EQUIVALENCE  (C I C 3NP  I  , C AlPMV 1 
EQUIVALENCE  ICIC3i1  I.OALPHP  I 
EQUIVALENCE  I C  I  0  3/0  l .  CPt-  IP  l 
ECuI VAt  EncE  IC I  l 1 1  1  I .CSORP 1  I 
ECU  I  VALENCE  I C (  11  l  ?  I  #  CSURF  2  I 
ECuIVAIEncE  IC I  1 1  10  I , SCURF  3 ) 

^  equivalence  io  i  120  i.bsuafii 

CPPINPUTS  PRC*  HUN  PROCRAH 

EQUIVALENCE  I C  I  2000  I  •  T  1 

ECUIVAIENCE  10  2661], OER  I 

CiiCC'TPuTF  -  CCEEFIC1ENTS  FOR  BQCY  FORCES 
ECUIVAIENCE  <0  12031, CX  | 
ECUIVAIENCE  (CI1212UCX0  1 

EQUIVALENCE  Id  12131, CXC  ) 

ecuivalcnce  io  12111, cnpt  1 
ECuIVAIEncE  ICI  12151, CV2  I 
equivalence  10  12  30  i  ,  eoefi  1 
EQUIVALENCE  <Cll’35),tCCN  1 
ECUIVALENC?  ICI 12111, CNPU  ) 
EQUIVALENCE  (012151, CYPU  I 

c 

C««CUTPLTS  -  CCEPFICIEMS  EOP  BOCY  HOHENTS 
EQUIVALENCE  ICI  1206). CLP  I 


EQUIVALENCE  10  12021, CHO  I 
EQUIVALENCE  1012061,  CNR  | 
EQUIVALENCE  201200). cl  I 
EQUIVALENCE  1012101, ch  ) 
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nnrf'O^nof' 


ecuivalence  iciijiu.cn  » 

ECulVALENCE  ICIWWI.CNO  1 

ECulvALfNCE  ICIIJIII.CNJ  » 

ecuivalence  iciwjw.cccn  i 

e  CU  I  V  AL  CNC  6  ICl 12CO I.CL2  I 

ecuivalenci;  ict  nm.cu  i 

ECulV*l.6l.CS  ICl  lJwll.CNP  I 

ecuivalence  iciWabi.Cnp  l 

ECUIVALENCE  ICU2<iS>.CLR  1 

C..CUIPUI1  -  COEFFICIENTS  FOB  SUB  FACE  EFFECTS.  ANO  IQTAL  EFFECtS 
FCulvAlENCE  ICIWOAl.CY  1 
ECUIVALENCE  IClWOSl.Ct  I 
ECulVALENCE  ICl  120S 1 .CL  > 

ECulVALENCE  ICl  WIOI.CH  1 
ECUlWAlENCt  ICIWUI.CN  ) 

ECulVALENCE  1C  I  W  IN  I ,  CEO  ) 

ECulVALENCE  ICIW2CI.CER  1 
ECUIVALENCE  ICl  W2W.C-0CP  I 
ECUIVALENCE  IC(W2?I.C“B  1 
ECulvALfi.ee  ICl  1223*.  CVA  ) 

ECulVALENCE  ICIW3«1.C<6  > 

ECjIvAlEnCC  ICl  W2EI.CL CAP  > 

ECulVALENCE  let  W2  A  I.CNQ  ) 

ECulVALENCE  ICIW2M.CLC  » 

ECUlVUEl.ee  ICl  1228  1  .CLNP  1 
ECUlvALfcNCE  ICl  W2NI.CLNP  I 
SCUIVAIENCE  ICIW32..BCL  1 
ECulVALENCE  ICIW3JI.BCH  I 
ECulVALENCE  ICU2JAI.8CN  I 
ECulVALENCE  ICl 1250  I .CEE  I 
ECulVALENCE  ICU2SW.CVP  » 


Input  vabiable  VINTER  IS  T6E  interpolation  control 
LESS  IRAN  1 6  BO  -  STRAIGl-r  LINE  INTERPOLATION 
POSITIVE  -  PARABOLIC  INTERPOLATION  .  WITH  END  INTERVAL 
INTERPOLATION  10.  TO  I.J 
0.0  -  STRAICHT  LINE 
l.o  -  TULL  PARAEOLIC 


IE  IT.lE.CER)  UT1N6-  0. 

If  It-uTlPE  .LE .  0.1  RETURN 
UTIPE*  T 

C  NULTIPL6  ANCLE  FORMULAE  ANC  ABSOLUTE  VALUES  OF  ANOLG  Of  ATTACA 
US  PE  I  •  S  INCI8PMB  I 
UCBPI  *  CCSCI3PE IPt 
US  2  PN*  I  v  SINC  I  2-  •  SPUR  1 
US2PH2  •  U>20HJ  ••  2 

us^fhi  «  $i?»c  (  *.  •  6P»-iP  j 

IF  l  ieSlPM-P^l  . Qt • 20*  IdAlPFY-SIGNt  2G*«0ALPmY1 
.GT  •  '"‘•JOALPFA^Sl  GN I  20  •  »  6  AG  P  HA  I 
UAL  PHA  -  A0S<0A».Ph 
UALFHY  •  A3SI6ALPHVI 


c 

c 


CALCUL A  T  |C*  OF  pOCY  FCBCG  COEFFICIENTS 

BGEFL  •  1  A  Q$  1  SSU*M  1  ♦  APS  I  B  SUftF  2  1  ^ABS  l  B$U&  F  3  >  »A8  S  16SG&F*))  M* 

C*0"CCC  IH2*  vpACh*  CxCAiCKOf  »hCXO»X|NTE<l»  3«C  *0  > 

y  *  gL  2  ISAtOKP*  VH  AC  h  » Cx  A  i  Q  XF  t  nC  X  ?  x  l  N  1 E  P »  3hCxC  #C  XC  I 
If | LCCNV.E0.2)CALLCUPT2 
If  lLCCNV.Ea.2IC  120001 -Cl  20011 


IFUCC*V.EQ.2)ft£TuRM 

CALI.  TAGL2  iaAtPFP,VMACH,CNAtCNf #NCN,X|NTER»AhCNPT»CNPT3 
CALL  1  A  PL  2  (  8AlPhP,V*ACH,C0CNA*CCCNF  ,  *:C  DC  N  »  X I  N  TE  ft  ,  AhCOCN  ,  COCNI 
CAwL  TAEL2  IBALPKP,VMAr.H,CV2A,CY2f  ,NCY2.XINTER,3HCY2#CV2) 

C 

Cx-Cxo  *CxC 

CKPl»CM’’T  *CLCN*US2PH2 
CYPU  •  CY2-USAPM 

c 

C  CALCULATICN  CP  BODY  *C*»rNT  COEFFICIENT* 

CALL  T  A  F  L  2  I  8  A L PK P . V* aCk, C L 2A , Cl  2 F . NC 1 2 . X J * T£ ft . 3HC 1 2 , C L2 ) 

Call  tael2  i 8 al p**p# v*»a>. Cl >a,Cl  )f.nCl  j, x inter  .  3hCl) #Cl3 I 

CALL  TA3L2  JBALPkP,  V-AC^.  C^A.C^PfNCf*,  *  l  N Tfc R  ,  3HC *0 , C H 0 > 

CALL  1  A  Bl  2  I  d  al  P‘  P.V^ACH  .CUCPA.COCMp  .f.COC".  xiKTEnI«.HCOCNlCOCF,> 
CALL  TAFL2  MALrr-P.VHAC»^CS2A,CN2P.NC*:2.XlMeR.  3MC.S2*CM2» 
call  TACL2  I  BALPHp  .  VHACH.  CLPA  ,ClPP  .f.CLP,  x  IN  IE®  ,  ?MC  L?  #CLP) 

CALL  TACl2  luALPf  A,  vmaCh,  Cm  ;a,C^OP  .NCH  j,  x  in  I6«  »  3»*CkC,C>*2I 
Call  TA0l2  IUALMHY,  VMICH.  CM<JA,CM3f  ,NCMw,  X  INTER  ,  JHf.Nft  ,CNR) 

c 

Cl®  »  Cl2«USAPm|  ♦  CL3«b$PM 
Cpp  •  C*0*CCCs»cS2FF2 

CKP  *  C*»2»US<PHl 

c 

C*aCALCUI.AT1CK  CP  SURFACE  COEFFICIENTS 

8CL  •  I -eS‘J*M-fiSUPF2»8$L.RF  I'BSCftMk'A. 

BC**  •  I  CS=JRF  L »  dSURF2*dSURF  3-aSLRf  4)/4. 

BCA  «•  <-8SURfl*0SURF2-8SURF3»6SLRF4)/A. 

BCNP  •  -ecu 

CALL  TAELi  IBALPhP.eOM,  VMACH,C20A,C20F,NC/0»XlNr£R,3MC20.CIQ) 
CALL  T  A  8  L  3  I  8  AL  PH  P  ,  BChP .  V*  ACH  .CZO  \  .  C  l  OF  *NC  TO .  X  I  NTE  ft  ,  3HC2R.CIRI 
CALL  TAil?  <?alp^P,8C>.vmach,ChCA*CMCF  ,*,C  HO .  XI  N  TE  *  ,  ShC.  HOCP  .  CHOOPI 
CALL  T  A  :-L  3  M A L PHP . StSF . VH ACH ,C“ DA . CHfF  *NCM> .  X  I  NT  Eft  ,  )mCHR  *ChR  I 
CAcL  I  A  _i  L  )  IBALPhP,3D'i.  v- ACM.C2DA  .C20F  .NC /0 .  *  l  U  T E  ft  ,  JmCYR.CXRJ 
CALL  I  ALL)  |0ALPFP,J'.i;>*,VMACHiCfi>A,C2DF.NC/O.l|Kf£ft,3HCYC.CYQ> 
CALL  T  A«!L  3  1  6  At  P  h.P  ,  »•  CN »  YH  A  C  H  •  Ch  C  A  #  C  wp.P  »-NC  HQ  *  X I  N  T  t  R  *  ShC  LOR  P  #  C  L  DRP I 
CALL  fAEL)  I  0  AL  Pwp,  00.“.  VMACh,  CHOA  ,CMOf  ,NCMO.XlNr£ft,)HCKfl,CNQJ 
CALL  TALL)  IfiALPFP.aCL.VMACM.CLOA.CLOf ,  NU  0  »  X 1  N  T  E  ft  ,  3HCLO ,  C  L  0 ) 

c 

C*C«-AES(CNC» 

CHCCP--AesiCFOCP) 
cl  c  * as s { clo i 

CVR  *AQS  fCYA I 
CYC»ASS  C  C  VO  3 
C2R*AeS(C/RI 
c/c*Aescczoj 
CHR^-AeSlCMR) 
cl  c  rp« • Ags ( clorp  > 

CZP^CC^FU  •C2C"L,CPh!#9CH*CZR  •LISPm  l  • SON ) 

CyF  •  ICYF u»c Y«»uCPM I •ec*»CTC«U$PM’*OOM) 

CL  -  (CL ft'ClC«80l i 

C.  Y  f  «  c HP#  I  C^CCP  *UCP3  1  •  eCM- CMR* LSPX I  •  BO* ) 

CLNF»C^*I  ClCRP«UCPH  t»dON«CHO*USPHl*BO«l 

c 

C**AfcRC  COEFFICIENT  ERROR 3 
CX  •  C<  *  C x ERR 
CZP  •  C  2  P  ♦  CZ  ERR 
CYP  •  CYP  ♦  CY ERR 
CL  •  CL  ♦  CLERR 

Cl. HP  •  C L M P  ♦  CHCftR 

CLNP  -  CLNP  ♦  CNERR 
C 

CM  TIUNSFCRPAf  !0H  FROM  W|HD  tO  0QOY  AXIS 
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Cv  *  CYMUC*H!-C!P*u1PH 
C t  *  -C2P«UCPMl-CYP«uSPhi 
Cn  «  CLf*P*UCI‘Hl  •ClhP»USPHl 
CN  *  CLNP»UCPH|-CLPP*IJSPHI 
MIUON 

en  c 

t«*AE«G  FC«CE  ANO  HOMENT  KOOULE  OOOV  AXIS 
Sue*CuTlM  A2 
COPFCh  Cl  A) 10) 

c 

C««!NPUT  CAT* 

EQUIVALENCE  ICU106),AF  AREA) 
EQUIVALENCE  I C  I  non,RFLGTh  ) 
EQUIVALENCE  JCU30I  l.ftCElCG) 
EQUIVALENCE  1C  I  l  3  1  3  )  ,  RMCG  ) 
EQUIVALENCE  fCmi4l»ft’VCG  l 
ECUlVALEN-CE  ICtmSl»AT|CG  1 
EQUIVALENCE  ) 

EQUIVALENCE  ICinm.RML  ) 
EQUIVALENCE  4C(  1627  UAGRAV  ) 
EQUIVALENCE  IC(  350A  )  rOPINA  ) 

c 

C«*INPUTS  PICK  CThCA  NOOUL6S 

EQUIVALENCE  <C(0?OH,PCTNMC| 
EQUIVALENCE  <C(C?OM,VAIRSP  1 
EQUIVALENCE  <C<  3*01, BUT  ) 
EQUIVALENCE  ICI  3B0).AAnC0  ) 
EQUIVALENCE  IC||20!UCX  ) 
EQUIVALENCE  ICCI20O.CV  ) 
EQUIVALENCE  ICC  L20S1.CZ  ) 
EQUIVALENCE  IC( 1206  IfCLP  I 
cwUUaiUCc  iC<  UuMrCM^l  l 
EQUIVALENCE  (C<  L2GII.CNR  ) 
EQUIVALENCE  !C  (  l  204  I » CL  I 
EQUIVALENCE  tC  f  1 210  I  »CH  I 
EQUIVALENCE  CCU2UI.CN  ) 
EQUIVALENCE  ICI  1236UCNI  J 
EQUIVALENCE  ICI1237I.CN2  ) 
EQUIVALENCE  ICIU38).Ch3  I 
EQUIVALENCE  ICimtWCt-4  I 
EQUIVALENCE  ICIIM  |  l,F  UX  ) 
CCi;IVAlENCE  IC<  I«*12  1,M»>Y  ) 

EQUIVALENCE  CC«:4l)l»FTK2  ) 
EQUIVALENCE  <CI142?I'UCG  ) 
(CulVAiESCC  ICI  1628  I.DnaSS  t 
EQUIVALENCE  ICIUi*  I.CEA2J  I 
EQUIVALENCE  ICC  IMS  I.CPA33  I 
EQUIVALENCE  I C  I  1  7  39  )  P  ) 
EQUIVALENCE  ICI  l  7Q  J  )  »vQ  ) 
EQUIVALENCE  IC117W),*«*  1 

EQUIVALENCE  ICt  17497, F *JY  I 
EQUIVALENCE  ICI  1750  1. EMU  ) 

c 

C««otme«  cut  PUTS 

EQUIVALENCE  ICIUOOI.FXBA  ) 
EQUIVALENCE  ICMIOII.FTM  I 
EQUIVALENCE  4 C I  l)O?t,f20A  l 
EQUIVALENCE  ICm03l>FHXBA  I 
EQUIVALENCE  ICU30A  l.fNYBA  | 
EQUIVALENCE  I  C  t  1 305  I  •  Ft*  2  BA  ) 
EQUIVALENCE  ICI  13041, EMU  I 
EQUIVALENCE  ICC  1310), MH2  I 
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(CUIVALtKCe 
CCUIVM.6NCE 
ECU  IV ALfNCe 
f Cut  VALENCE 
ECUtVAltNCC 
ECuIVALENCE 
ECu  t  V At  f KCE 
EQUIVALENCE 


KmiUrf^J  l 
icinui»F«M  i 
tcimoi.fNKTN  i 
icti)2i i.fnyth  i 

let  1)2.  > 

ICI  l>2n.EHALUC  I 
ICE  I32A1  9FmvLUG  » 
tctnj^i.FHuucj 


c 

C**FC*CE  V6CKR  CCHPUNEnT* 


UCS  ■  PtrKMC »PF AREA 

ucsc  •  ucs*«FtcrM 


F*eA*uCS*l-CH)*FThX 
Eve  A*uvs*cr  4MHV 
F28A»uCS*Ct «FThZ 
IF  I  v  M  n 3 p  , v.  E-0 .0  )  CO  TO  72 

c 

C««AE»0  PQPEM* 

Fwiba  •  ICLAlClP/VAl«SP)«AfLGTM*P  )«UOSL 

*  •  ICHilC^C/VAtflSPI-RFLGTh^tfO^UOSLAF  rBA#RO€LCC 
tn  BA  •  (CM(CNS/VAIRSP)*RFtCTh#hR  UOQSl-F  YQA«AOUCG 
C 

C^HO^EMS  CAUSED  b Y  THRUST  HI  S  At  t  GNmEhTS 
fKKFH  •  *FTFY»AF/CC  i  F1kZ»AFYCG 
FPYTH  •  fTHX«RFKC  A  FlHZ«AFxCG 
FPZIm  •  -fthX*RfYCC  -  FfHY«AFXCG 
C 

(••POPE  MS  ANC  FORCES  CuE  TO  LUGS 
IF<  RASGC.LE  .fUUtALUGKOTOTO 
Ft luG  •  0% 

•  •»  0. 
fpxlug  •  o. 

fHYLUG  •  0. 

F>*ZlUG  *  0. 

CO  TO  7A 

70  IF  IPAACC  • L E*  RAILS  CO  TO  72 

FTLUG  *  -(FYEA  ♦  ChASS*AGRAV«CFA23  ♦  IFPZOA  ♦  FH J TH I • 
•  *L€G*0»ASSfFHUI/(  1.  ♦  0MAS3*RLCG*RLCG/FHtn 

FZCUG  ■  -(F26A  ♦  DPA$$*IGRAV«CFA)3  ♦  (  FMVBA  ♦  FNVTHlt 
p  RLCC-OHASS/FMIYI/U*  •  OmASS*RiCG«RlCG/Fm|  Y| 

FHXLUG  ■  -IFHX3A  «  FHXTHI 
FHYluC  a  FZ LUG»RL  CC 
FHZlUG  i  FYLUG»ftlCC 
GC  TO  7A 
72  CCKTIKUE 

FtLUC  •  -(F Y0A  ♦  CPA$S«4GAAV«CFA23) 

FZtoG  ■  - <  F  2  8 A  ♦  OP  AS  $ • ACR  AW«C  F  A  J  > 1 
FPXIUG  •  -tFPXSA  ♦  FHXTh) 

FPYIUG  ■  -tFPYQA  »  PMYTM 
F  ►  Z  LUG  ■  * ( F  HZ  0 A  t  FHZTp» 

7A  CONTINUE 

c 

C««TOTAL  force  A*0  MOMENTS 
FY0A  •  FyCA  *  FYLuC 
FZflA  ■  FZ0A  ♦  FZLUC 
F  P  L  8  A  ■  fmiA  ♦  FpxTm  ♦  FHXLUG 
FPYgA  *  fP  VBA  ♦  FHYTH  ♦  FHYLUG 
FPZeA  -  FPICA  ♦  fPZTH  ♦  FM2LUG 
C 

CAICUCATE  HlNCC  HCHf NTS 
F  PH  1  •  CH*UQM 
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Fpn2  m  0*2*UCSL 
Fph3  •  CFnuCSL 
FPN*  •  CF»»uOSl 
RETURN 
END 

CM  INI  T1  EC  1  2  AT  ICN  FOR  ENGINE  ROOULE 
SUeocuriNE  EJ| 

CD  ppCn  cnMOl 
C I  PENS  ( CN  IHI  ion 
ECUIVECENCE  ICI2S»1I,N  ) 

ECUIVECENCE  ICI  254  2  1 1  1RL  | 

Cl  1*99 1  •  0. 

IFllN  1  .  |«ft 
N  •  N*  | 

RETURN 

ENC 

CpiENCIne  PCCulE 

SU6ROUT1NE  I) 

coppcn  cmioi 

c 

Cmucok  uf  r»eiE  for  thrust 

COEFCN  /NTF/NTN  /TheAG/TFA  /THPUN/THP 

£«•  INPUT  c»r* 

ECUIVECENCE  ICI  1*01  ).  8ECPMT  l 
ECUIVECENCE  ICtl*02),»PHT  1 

ecuivecence  ici  hoi  i.cnecgni 

ECUIVECENCE  <CI  1*0«>,PCFTM  | 

ecuivecence  icmosi.oejRN  i 

ICOIVELENCE  ICIItHUtllE  » 

ecuiveleuce  ICI  nisi, chi  » 
eouivei t nee  tcinn>,o*p  i 
ECU  I  v  Ec  ENC  E  ICIIMTl.RCCOC  I 
ECU  I  VEtENC  E  IC(  Hie  I  ,PCCGF  I 

ecoivucnce  icinj9),FH|<o  i 

EQUIVALENCE  ICI 1*20  I.Fmiyo  I 
ECUIVUENCE  ICI  l*21).ALCGO  I 
CCUlVELENCE  ICI 162 7  I ,  ECR A  V  I 

CM  INPUTS  f«c«  OFFER  POCULES 

ECUtVElENCE  ICI 12S21.XINTER> 
e  EQUIVALENCE  ICI2C00I.T  I 

CM  OUTPUTS 

EQUIVALENCE  ICI  llOSI.AOELCCl 
ecuivelence  ici  Homuc.p  1 
CCulvEtENCE  ICI  1*10), FTi-RSTl 

ecuiveceiice  ici  mi  i, fine  > 
ecu i v ec ence  icmm.im  i 
ecuivecence  ici  ui ? i.fth  i 
ecu i vei ence  ici  nm.Ricc  1 
ECUIVUENCE  ICIl62i),OP*SS  | 

ECUIVECENCE  lCll7N8l,F«i*  ) 

ECUIVECENCE  ici  1748 1 ,FF|T  | 

ECUIVECENCE  ICtlT5Cl,FPIl  I 

CM$r«7E  VARIABLES  ENC  TPEIA  CERIVETIVES 
EQUIVALENCE  1C  I  1*96  I , u IPPO  ) 
t  ECU  I  VAC  ENCE  ICI 1*991, UI"P  I 

IF  ICEURN.CTkO. I  RC7URN 

F?NR5l-CC0IP2iT,  7FEf7MF,HTN.X  INFER,  EHF7HS  HI 
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fTHRST  *  FTHftST'll.  R  PCFTHI 
C 

io  u$ina*$inci enPHti 

f  Tm*  FT  NR  ST  #COSCC  QALPhT  ) 

ftnn«-fthhst  *U  S I  N A  *  $ INO I bph  I  t 1 

FTh*»fTNRST  «OS |NA*C0SCI8PMlrl 
CC  TG  JO 
20  F  IHS^FTNRSr 
FTMY^O. 

FTb2*0. 

30  CCMINUE 

c 

UII*P0  •  fthrst 
UCWf  •  UINP/C1SP 

c 

CPASS  •  l C w  T  -  UCwPI/AGRAV 
ROELCG  4.  RCCGC  -  (RCCCQ  -  ROCGF )«U0*P/D*P 
C 

ff»U  •  FNIXCMCWT  -  U0*P)/CwT 
F*IY  ■  F*lYO«tOuT  -  UOWPl/CbT 
F*U  »  FWIY 
RL  CC  ■  »LCGC  *  ROELCG 
IF  1F1KRST  •  C-T  •  0.  J  RETURN 
C 

WRITE  16. 1001  T 

100  fCftPM  (//MM  ailRNCUT  T  iHE«*F8.4,5M  SEC.) 

ceuftN* i .o 

FTHRSI.O. 

F  T  H  X*0  a 
FThY«0. 

F  Th  A«0  • 

RE  f  u*N 

two 

c«*  T R A N$i AT l CN At  DYNAMOS  INITIALIZATION  NODULE  FOR  01 

SUBROUTINE  Cll 
CGRRCN  CI4J10J 

EQUIVALENCE  1C12561UH  » 

EQUIVALENCE  ICC  2562UIPL  I 
CHEKSICN  I  PL  (100) 

C 

C«»  INPUT  CATA 

ECU IVAlENCE  (Cl  100). W*E  ) 
equivalence  1C «  10 1  >• vwr6  I 

EQUIVALENCE  1C  <  102). VUE  1 
ECU (VALENCE  (C  <  20 A  )  , VM ACM  1 
ECulVALENCE  ICI  20ei.Rh£S0  ) 

EQUIVALENCE  I C 1  36.?  )  .OALPMA  t 
EQUIVALENCE i Cl 3?  1». RANGE) 

ecuivAicNceici  1752  i.bpw  io) 

ECulVALENCE  2  C  C  ISO!  i,EMARG) 

ECulVALENCE  (C (  36  E  ) . fiALPHY ) 

ECulVALENCE  I C 1  4271.BTFTG  1 
ECulVALENCE  ICI  421),BPSIG  1 
ECulVALENCE  I C 1 1403 )#CCU«N  ) 

ECulVALENCE  IC<  163<?1.0PTARGI 
ECulVALENCE  ICM666  ) ,  BC  IVE  ) 

ECulVALENCE  1C1 1667 ) .RSLANT ) 

ECulVALENCE  ICUT39)»WP  ) 

ECulVALENCE  ICI 1743), «Q  1 
ECulVALENCE  (Cl  17  4  7  1. MR  1 
ECulVALENCE  ICI >502),0PTN2  1 
ECUIVAIENCE  ICI 3304 l»G*TN*  I 
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no 


c 

C»  CUTPul  1C  PCOUlES 

ECUIVALENCE  (C(  S70I*BPP1P 
ECU  IVALENCE  ICII615).R«E 
ECUlVALENCE  ICt!61<»l  aRVE 
EQUIVALENCE  ICU62il,RIE 
ECUlVALENCE  ICII605I.VRE 
EQUIVALENCE  ICIItCtl.VTf 
ECUIVALENCE  ICI16U»aVIE 
ECUlVALENCE  ICE I62VI.AEHOSI 
ECUIVAlENCE  ICI  1655  I.RCEl* 
ECUlVALENCE  ICI  1616). RCEL7 
ECUIVALENCE  ICI  16111. RCEU 
ECUlVALENCE I  Cl  UNNI.ATARGt 
ECUIVAlEN  .  ICI  1657  I.VTAAG) 
ECUIVALENCE  ICI  16511. RT»E 
ECUlVALENCE  ICI  1655  I. RITE 
ECUIVALENCE  1C  I  1656  I  . RI2E 
ECU IVALENl  '  ICI I66E IfrAXD 
ECU  I V  ALENL  ICI1666WRY0 
EQUIVALENCE  ICI 16701. R20 
ECUIVALENCE  ICI 17551 .BTHTO 
ECUIVALENCE  IC<1T5*>»BPSI0 


C<16N7HC<16<I6  I 
l P  L INI  •  1600 
I  PL  I II*  1  l  •  l  60N 
I  PL  I N  *  2  I  •  1600 
IPIINMI  •  1612 
IPI  (NANI  •  1616 
I  PL  I  ♦  5  A  •  1620 
IPI 1 N*  6  I  -  I 6N0 
IPLIN‘1)  -  li** 

IPLiNNBl  •  16»B 
IPLINA9I  a  1652 
I  PL  I  Pin  10  I  •  1656 
1PLINMIH16T2 
N«  N  n  1 2 

IFISP'TAPGIIO.V.IO 
10  ATAOG-AINRST/ENTARC 
CCTOll 

9  AIIPG-0. 

it  continue 
BPH  IC-O. 

C*aCalCulaT£  PI  55  I LE  PARAKETER  INITIAL  CONDITIONS 
BH-IC'BOIVEaSAIPMA 
BPSIC--EALPNY 
R«€n-RSIANIaC05HBCIVE» 

RAhCE • R5L  AN  F 

RVEnRSLANTaSI'ICIIIALPHT) 

R2  E  a  NSLAM  aSINCIBCIVEI 
20  Rh  a  RH 2R0  •  R2E 
CIN27l"-eALFHA 
CINIIIaEALPPT 
CIOll-'BALPPT 

c 

USTNT  •  SINCISTMOI 
UCTHT  A  COSCUTHTO) 

UCPS  1  •  COSO  I  EPS  10  I 
USPSI  •  S 1NCIBPS 10 ) 


«*8«  .  -UOII.UCThT.RXE  .  UJTMT.Rte 

USPSMRXE 

«ZB»  •  -UOSI.uSTHr.RxE  -  UCTHT.RIE 


-  - 

VHhTE  .  V  NACN . VSOUN  C 
VH.XY  .  VHWTE.CCS0(8AIPHA  ■ 

vxe.vwxE.vmvxy.coscibalrhyi 

VYE.V.XE-VHk.Y.SINCIBALPHYl 

c  VTE.VvIE-vmvIe.sinCIBOIVE) 

JO  RCEIX  .  RTXf.RXE 
*CEIY  »  RTvE-RYE 
ROEtt  . 

RXO  •  RXE 
*YC  .  RYE 
R2C  •  Hit 
RE  Turn 
eno 

C..TRANSLAIICNAI  CTNAHICS  H0CUL6 

sueRcun#.E  ci 

CCCHCH  Cl  A3 101 
C.RISPUT  CATA 

ECulVALBNCE  ICI  14211, RAIL  ! 
EQUIVALENCE  ICI 1 62  7  I ,AGRAV  | 
ECU  VALENCE  ICI  1628  I , CRASS  1 
EQUIVALENCE  ICI 1629 l,ATHRST 1 
EQUIVALENCE  1C  I  1633  I , AIURHT 1 
ecu  VALENCE  1C  1 1 6  3 1  I , BC AMf  1 
ECUIVAL64CE  ICI  1480  I.IIPINO  1 

,..u , ehce  ;c;  lioi  i,ac ive  i 
ECUIvaiehCE  ICI 175 l l.CRAO  I 
£  ECUIVALENCE  1C  I  3504  I »0RTN4  1 

C.RINRUI  5  ERO  CIEER  NCCCLES 

EQUIVALENCE  ICI  3  7 | 1 , RANGE  1 
equivalence  ici  oaoi, banco  i 
EQUIVALENCE  ICI 13001, FXBA  ) 
ECUIVAiEnCE  ICI 1301 |,FYe*  l 
EQUIVALENCE  1C  I  I  30  2 ) , FEB*  1 
ECU (VALENCE  ICI1667I.RE  J 
ECUIVALENCE  ICI170D.CFA11  1 
ECUIVALENCE  ICI  1707  I.CFA12  ) 
ECUIVALENCE  ICI 17111, CFA13  1 
ECUIVALENCE  ICI  17151,0*21  | 
ECUIVALENCE  10  1715  1,0*22  I 
ECUIVALENCE  ICI  17231,0*21  I 
ECUIVALENCE  ICI 17271, CFA31  1 
ECU  IVALENCE  1C  I  1 7 3 1 l.CFA »2  I 
ECU (VALENCE  ICI  1 7 35  1, CFA33  1 
ECUIVALENCE  1C  I  2000  I , T  | 

C*»S(ATE  variable  OUIRUTS 

ECUIVALENCE  1C  I  1 600  1 , vxf  D  1 
ECUIVALENCE  ICI1603I,VXE  I 
ECUIVALENCE  ICI  160*1,  WED  | 
ECUIVALENCE  ICI 1*071, VY£  | 
EOUIVALENCE  ICI 160  8  I , VE EO  I 
ECUIVALENCE  1016111, V2E  I 
ECUIVALENCE  ICI 16  12  I ,RX£0  1 

ECUIVALENCE  ICI 1615  ),RXE  I 


6THT01 
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ECUlYftCENCE  1C  I  1 6 1 6  I » AYEO  I 

ecuivucesce  jci  > 

ECUlYftCENCE  1CI  l®*0t.«*i£0  I 
ECUlYAlEhCE  1C  I  l*2!l.»tE  * 

ECU  I  VAC  ESCB (Cl 16*1 1>V4 »«GI 
ECC.IVAEESCEICI16T2  ),S>  .101 

CCUlv*liNCE(CU6T5).0f>S!T  * 
ECUlVAltSCE  ICDSOll.ENT**C» 
ECUlVACESCE  ICI  l6ftft)^t«E0  I 
tCUlvftlBsCE  (Cl  It'S l  I  .ftTXE  > 
ECUlvftlENCE  ICI  1652  I.RtTED  l 
ECUtVftlESCE  ICI  16561. «TYE  » 

E  CU I Y  AC  ESC  E  ICI 1656 l.ftl «E0  » 
ECul YAlENCE  •'( 16591, »rtE  > 

C 

C»ft07i>e«  cltputs 

ECulYAlfhCE  ICI  ll!M.I<l‘  > 
f CMYAiCNCE  ICI  1625I.AYBA  > 
ECUlYftCENCE  ICI  16261  . »/B»  » 

ECUlYACENCE  ICI 16)21. VCEEX  » 
tCLIvftttACE  ICI  16)2 l.VCECY  ) 
ECUtYAlENCE  ICI 16)6). VCEll  I 
ECUIYACESCE  ICI  16)61. ECECX  1 
ECUIYAlElCE  ICI  I6)M.RCEH  * 
ECulYAcESCE  (Cl  16)11. BCEC2  ) 
ECC I  V  Ac  ESCE  1C  I 1 6)0 1 • VC15NCT 
EC'J  l  Yft  l  ESCE  ICI  16601. Vl*E  1 
ECuivfttEuct  icii66ii.vr»e  l 
EQUIVALENCE  ICI 1662). vite  * 
EQUIVALENCE  ICI  166)  1 1  VMS  • 
tCu.v»l.6hCE  ICI  U6*I.Yil<«  A 
ECU  I  YM  EsCE  ICI  1666  I.YOI0  » 
ECulVACEsCE  ICI  16)61.  ASCII  ) 
ECU  I Y AC  EsCE  iC I  U17 ) .ASCY  l 
ECUIVAlENCE  ICI 16T8I.AN6J  ) 


CAOO^^o'^-  EOIICES  10  GET  TOTAC  .CCECOUT.OS  IS  OOOT  »«* 
AlEA  .  EXEA/CMASS 
AYEA  »  EyEA/CSASS 
ACOA  »  FZEt/CSAES 

Cftft«E$CCVE  EMCK  8CCY  TO  E*RTN  AXE* 

>»«E  •  CF All Y*xa»*Cf a21**y»a«CFA)Iya26» 
aye  •  Cf A12<A*BA.CE A22*»r8A«Cf A)2«A26A 
A*E  •  CEAl)YftA0A*GEA23*AY0Af.CFA))*AE8* 

C.*l NT  EQUATE  ACCEIED ft T IONS 

VAEC  •  »»E 

YYEC  •  »1E 

VZ^C  •  «  AGRiV 

Cftft  CftLCUCftTE  TOT  AC  "ISSUE  *CCEcE**TIOS  IN  BOOT  »AES 

Y0A8  •  Cf  A 1 1  ft’/AEO  *  CF»12*YTE0  ft  Cf*lS*v2ED 

YOYt  •  Cf *2 1 YVft  E C  »  CFAJJyYTEO  »  CfAJWlEO 

YCE»  ■  CfAllftYAlC  •  CfftJJftYTEB  ft  CFAJ)«y2E0 

If  (  AGMV.LElO  *  )  GO  TO  10 
ASGX  »  VC*e/AG«*V 

JhCY  •  VC*e/AG«*V 

ANGt  •  VC1»>*G*AV 
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C««!hT6CR#T6  vf IOC  I T  I  ES  TO  EARTH  AXES  P03IT10N 
AxEC*vxp 
ftY£C  •  WE 
*ZE0  •  VZE 
!F(EMTARGUO«4»10 

10  AT  AflC-ATMST/EHTAAC 
cqicii 

9  ATARC'O. 

11  CCNTINUf 
EPS  ITC*  0. 

IF  ( VT  ARG«GT »0* I  BPJITC*  A  TURN  ? »AGR  A V*CR AQ / V  TARG 

c 

RTXEC  •  VTX  E 
RTYPO  •  vTYE 
RTZEC  ■  Vr2E 
C 

vceix  •  vtxe-vxe 

VCEtY  •  vTY€-vr6 

vcei /  ■  vtze-vee 

c 

RCEl*  *  rtxe-axe 

RCEL**  a  ftlYE-AYE 
fl  C  E  l  Z  ■  AWE-R2E 

VCtSNG  •  <RCtLX«vO6LX*R06Lr*YO£tY»AOELZ»VOeL2)/RANG6 

RETURN 

EnC 

C**R0TA7ICNAl  CYNAHICS  IN  I  T I al IZ AT  ION  MODULE  02IEUL 
SU9RC0T  1st  02| 

CCppCN  C 4  *  3  1 0  I 
ClUNStCN  1PC  U00> 
c ••NChE 

£•• Inputs  frc*  pain  program 

EQUIVALENCE  TCI2561WN  1 

E  CU I VALEnCE  ICIZ562I.1PL  1 

C«»S1ATE  vJRMCLE  OUTPUTS 

EQUIVALENCE  <CU701)«Cf  All  1 
EQUIVALENCE  ICI  1707  I.CFA12  ) 

EQUIVALENCE  f  C I  1711  UCFA13  » 

EQUIVALENCE  ICI  1 7 1  5  I  • C  F A?1  I 
EQUIVALENCE  ICI  1?19)»CFA?2  ) 

EQUIVALENCE  JCI1723 »rCFA?3  I 
EQUIVALENCE  ICI  17?n.CFA)l  1 
EQUIVALENCE  ICII7J1 J.CFA32  1 
EQUIVALENCE  ICI  17)S),CFA33  ) 

C**CTHER  CuTRuIS 
C«*NCHE 

C**lMT!ii  CALCULATION  Of  EULCR  ANGLE  MATRIX  OF  DIRECTION  COSINES  4CMI 
USPhl  •  SINUI6PHI0) 
uc ph i  •  casciBPFiO) 

USTbf  •  SlNC«OTMQl 
UCThf  •  CCSCI6TFT0) 

uspsi  •  si-icibpsioi 

UCPSI  •  CCSCIBPSIOI 
CF Al  |  ■  UCPS  |  «UC  TFT 
CFAI2  •  U  5>  P  $  I  •  U  C  T  h  T 
CF  A  1 3  ■  -uSIbT 

CFA21  •  -uSPS!«OCPfIAUCPSI«U$ThT»USPHI 
Cfk?7  •  u«'  PS  !  auCPb  I  *USP$  t»USTWl*USPHI 
CFAJ3  *  UCTM  -uSP*  I 

CF  A)  l  -  L'CPS  I  ■uSThT«UCPhI*,USPSI»U$PHI 
C  F  A  32  •  uSPS  I  aUSTMaUCPhl-LCPSI*USPHI 
CF  A  3  3  -  UCT  ►  T  »UCPh 1 
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1 


I 


\ 


t»»tnter.B»tec  u*t  MBL»  eo*  wPD.xoo.^o.iHO  Cfw 

t*iOD  <•  1700 
ftklh*  )  «  170' 

1 PL 1 h  •  i  1  *  1  70S 
IPUM))  •  17W 
I  PL l**'  l  •  17V6 
I  PL  IN* 5)  m  1720 
| PL  I N* 6  1  «  172' 

I  PL  IN*7  1  ■  l7?® 

J  PL (N»d  1  •  1732 

IPLtN*"  ■  l  7  J' 

J  PL  I  10  *  *  17'0 
t  Pt 1 N ♦ 1 1  7  *  77" 

<«,  BESS*  »KCUl»»  «A»l  0W»V»TIVM  to  ie*o. 

C (  1700 1  *  o. 

0. 

0. 

0. 

0* 

•  0. 

.  Q, 

•  0. 

»  o. 

»  0, 

*  0. 

*  o. 


CIITO" 

Cl  170V 1 
Cl  1 712  1 
CU7161 
CH720  ) 

C 1 172'  1 
Cl  1728  1 
C  1 1 7  3  2 ' 

C  1 1716  1 
C (  l  7 AO  7 
CUT"  I 
PE  TURN 

SCTaTICNAL  CtUAHtCS  PODVJLt 
SuARflwT  C7 

CCP'Ch  C I  S3  10 1 


C«#0P7*  I KPU  T  3 

ecu IVAlENCf 
£CU!V*L  ENCE 
t  cu 1 V*L  est€ 
eCUlXACEKCC 
ecu lv*L€NCE 

ecu ivale*ce 

equivalence 


•  C t  l'2l>»^*77.  1 

IC<  IT'BU  PHI*  7 

«C 1  17"I. P*Mt  ) 

ICIlT^OJ.^U  1 
icu7^n.c«A0  i 
ICt  J80JI.ORTN31 

let 250')r0PtN'  1 


C##|HPUTS  fBCK  CTvEft  *°ci16*  rrt  . 

E QU l Vi L E'CE  ICt  30OI»R"kO  7 
E GUI V *L E^C6  10  1  1305  1  .^*8*  » 
CCUlVAlENCE  4C 1 130' 1 » FH70A  ) 
EQUIVALENCE  1C l l 335 1 ♦ F*l BA  l 

EQUIVALENCE  tC «  13081  *  acei.ee> 

c  •INPVllS  MC#»  rMN  PPOGAAN 

r «*S 1A 7 E  VA«M8lE  CutPUTS 

equivalence  ici  itooucpauci 
equivalence  tci  no; wce au  i 
EQUIVALENCE  ICC I70'l.cf A12C1 
€CulV*LENCE  cr.l  17071*07  A\2  * 
EQUIVALENCE  ICI  17081. CfAlJCj 

EQUIVALENCE  ICI  ITllJ.CMll  » 
cquivale^ce  ici  i712)*cea2ici 

EQUIVALENCE  CCl  171^  1#CPA?1  * 

eCutV*LkNCe  IC'- 1716  )  .CP422C  I 

EQUIVALENCE  ICI  m*UCf>22  1 
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€CU!VAieNCE  TCI  17?0),CFA2301 
EQUIVALENCE  ICC 172?  I.CFA23  J 

equivalence  ici  i  u*  )*cfahd  1 

EQUIVALENCE  f  C  4 (72?),CfA31  ) 

EQUIVALENCE  IC  I  17  2  8  >,CFA32C ) 

equivalence  ici  irn  >*cfa32  ) 

EQUIVALENCE  t  C  (  1?)?)»CF  AJ30I 
KUViLtNCP  ICI  1715  »,CFA3)  I 

EQUIVALENCE  ICI1736J.WP0  I 

EQUIVALENCE  ICI1T3S1.WP  ) 

ECU l V  At  f  Kf-E  1C  I  l  MO  1  *WCC  ) 

EQUIVALENCE  ic (  LT4 J),wQ  ) 

EQUIVALENCE  ICl  1  7*4  I,wRO  ) 

EQUIVALENCE  «C(  17*91, Mil  l 

c 

C**!NTEGftlTE  eCCY  ANGULAR  RATfcS 
WPC  •  CftAOf  *MQA/EMX 

55  tfOC  -  tCftAO«EHveA*tFH  f  I -EH  I  X  I ■ wP*W*/CRAO I /FH I V 
65  HftC  »  CCRAO«FMZCA»  (  FMtX-FhlT  l»>»P«WQ/CRAO»/FMW 
C 

C**  I  NtECR  AT  E  A  f  T  l  T  C'OE  OlRECrtON  COSINES 
*9  CPillC«(CFA2l*«:t-CFA3l*MOl/CRAC 
Cf A12C«CCF a22*w«-CFA32**GJ/CRAC 
CPA  I 3C« (CF A?3»WR-CTA13*WC)/CRAC 

CF  A  2 1 0  *  ICf a>1«wP-CFAII«wRI/CAAO 

CFA22C  •  <CFAi2«wP-CFAI2*wR)/CRA0 
CFA23C  •  IGF  A)  J*v,p-Cf  A13»wR  l/CRAO 
CfAJlC  •  I C  F Al l -wg-CFA21»WP l/CHAO 
CFA320  *  (CF  AW«KC-Cf  A22«WP  l/CRAO 
CFA33C  *  C  Cf a!3«kQ-CFA23*wP)/CRAO 
RETURN 
*NO 

C  BASIC  INPUT  SUBROUTINE  OlNPTl 

SUBROUTINE  CINPTl 
COPPCN  Ct*3lO» 

EGUlVAlENCECCI  2BOO  1.6) 

EGulVAlENCElKC  4310  ) .  KIC  > 

EQUIVALENCE  ICC  3218J.CNAHE1W  <  C I  3268)  ,CNAPf  21  ,  IC  I  331 81  ,0NA#*E1  > 

C  Id  J32B  I,Onake*U  ICI  2  36U  .NCMCO  1,  IC  I  2  362)  ."CONG  I 

C  ICI  34*0  WNORNQNI,  I C I  34* l ) , ANOpNO)  ,  IC  (  3  1  6  7 1  .  NOJUt  I 

C  ICI 3168), OUTNQ  ),  ICI 2*61 ) ,NCSU6  ),  <C I  2 462 ) , SU3N0  ) 

C  4Cl  2.355  ),|R  I,  ICC  23  57 1  e  VR  ),  IC  I  3 1 39  J  .NOS  T  A?) 

C  JC(  333B),lCSlAT  I,  ICC  3340  ,S!ATNO»  ,  (C  (  3C6M  ,NUL  I  ST  I 

c  iCl  3067  1  ,  USTNCU  ICC  31 1?) .vallE  ),  !C  I  2  GO*  J  ,?vOTNC) 

C  ICC 2009I.NCP10T ),  ICI2325) .ylahlE; ,  <C,M 

EQUIVALENCE  lC(  190A  l.NPLOT  ) 

EQUIVALENCE  ICC  1965 ) ►OUTPUT  I 
CCl/PLE  PRECISION  ALPHA 
CIpEnSICN  Ona«C 31 IO),ONAHE*l 101 
OPENS  ICN  LIST  »IC !  5C  ) ,  VALUEC5QI 
01  PENS  I  CN  SUCVJI99  >•  1RI2).VR(2> 

OpEnSICn  RNC^NC I 5C  I 
ClPENS ICN  91  ICO  ) 

OpEnSICn  ALPKAI  3),CNANEH50UONAH62I50»,OUTNOI50I  ,*00nQ( 991 
OlPENS  ICN  KIOLOI 
OpEnS  ICN  STATNC  I  100) 

OlPfNSICN  VLA9LEI2.15) 

0  I  PENS l CN  GUTPLT I  151 
REAL  PCCNO 
INTEGER  OUT  NO 
INTEGER  ANOPNO 
Integer  cutplt 
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IMSCCa  STATNO 
lF(*fO101M6»M«96 
35  KK«0 

3*  consul 

JA»  »  0 
**17616.31) 

31  f LP^A M 11H1  INPUT  tATA/> 

1  ftoits.n  i  r  i \) ,  at  ( 1 1  .alpha  i  2 1 ,  alpha  1  3 ) ,  1  r  <  2  1  ,v«  1  n  ,  vri  d 

Nfi I T e  t  6  ,30  )  I  PI  l I , ALPHA!  1)  .ALPHA  12) .ALPHA! 3  I » 1  A  I  2 ) * V* III , VR  l 2 ) 
30  FCPkaT! 12.3A6.I5, 5X,1P261S.?) 

2  FC«nat  I  12,  JA6,  IS,  ^X, 2615, 9) 

7  m  [Kill  .N6.  I  )  CO  10  3 

NCSUO  •  NCSlf 6  A  1 
SteKC<sCSU8»  •  IM2) 

K  1  •  NIC  «  3510 

C !  M  )-  tM2> 

MChl  )"24*l+NQ3UQ 
KR • K  K ♦ 2 
CC  TC  l 

3  (M  IR  I  U  .N|.  2  )  CO  TO  4 
NCR'CC  -  NONCC  *  1 
>'CC».CIACHCO)  •  UU  2  I 
Ki«*KOSKO 

MKl*l  )  *  236  1  ♦NQMOC 
KK  ■  KK  t  2 
CC  TC  l 

4  If  I  mi)  .Nil  3)  CO  TO  <5 
L  •  IR12) 

cm  •  vriii 

K1 *KA ♦ 3310 

Ccr.mvAs;) 

KIN  l* 1 UIRI 2) 

K  K ■ KK  a  2 

IF  IVRI2)  .10.  rt. )  CO  TO  I 

NCI  1ST  •  NOL  IS .  *  l 

LI  STS Cl  NCI  1ST)  «  l 
ViUf  ISCLIST  I  -  VRI  1  I 
*  t -r.ic ♦  >5 10 
r. :  •  l  i  » l 

.<«“  1  I  •  3G66*NCl!ST 
C 1 K  L  ♦  2  JaVRI  1  ) 

K(Kl»3)aJll6»N0LtST 

R  K  a  KN  *  4 
CC  TC  i 

5  IF  I  I«U  l  I  .N!  .  4  )  CO  TO  6 
NCCUT  •  NCCV.T  ♦  1 

IF  ISCCUT.CT.SO )  CC  TO  * 

CNAPt  l  INC  CU  T  )*ALPnAI?) 

GNA^IKCCUT)  •  ALPHAI3J 
OUTNCl.NCCUT  I  N  |R  I  2  I 
M 3510 
CIM  )  •  ALPi-Al  2) 

K  I  *  l  ♦  l  )  a  32 l 7  ANCOUT 
C  1  A  l  ♦ 2  1  N ALPHA  I  3) 

MAI*  J  1*326  7  +  -SQCUT 
Kt*l  *4  )  ■ I R I  *  ) 

K!Kl*5)NJ167 ♦NQOUT 
MR ■ KK  *6 
CC  TO  1 

6  IF  IIRI1)  .hi.  5)  CO  TO  16 
IF  I VR 111  .60.  0.)  CQ  TO  IT 
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LCSTaT  «  10STAT  •  | 

17  hCSTAT  *  NCSTAT  •  I 
ST aTnC I KCST Af )  •  I  ®  f  2  1 
CKA^f  3  m~ST  il|i  AtPF  A(  21 
C  N  A  P  £  a  l  hCST  If  >•  AlPMt  31 
CC  TC  L 

14  IF  MM  1 1  .Ng.Tl  CO  TO  19 
SPLC »• A P^Ct «l 
IP  UPlCT.GT.m  CC  TO  1 
CC  20  Ul»2 

20  VliHE  l  l  *NKOT  )•  ALPHA((«I) 

CUT  Pi T (NPCCl ) • I R (  2  I 
CC  TO  l 
19  CONTINUE 

K  1  4  30 ^  )  •KCSC'B 
M43Q4I«KCI*CC 
HU  307  J-NCl  1ST 
K(  4  3C6  )  «NCO(JT 

h ( a  ) *ics i tr 
M  4  3C*  >  •SCSI  AT 
K14303  UNPl.CT 
K143e2>ANC*NCM 

IF  I !»< 2 ) .£0.0) a€TURN 
h  -  IRI2J 
CC  12  I  -  1*  N 

«£AC  14.13)  J.Y»HANO«m(ER»SIGNO»BETA 
*»ire  16.13)  j.T,j-ANo.Hi€«,sicso^eeu 
13  FCPpbT  I  15,  f  5.3*2<  3Jf.U2).  2615.41 

•  NCAKcH  •  \ 

ft.NCPNC(l)  ■  J 
C  {  J>  •  3 
KlJfl)  A  HAM*. 

M  J  •?)  A  H  I  i* 

cuo:  a  sicno 

M-*<*  )6lO 
MM  )*f»AKC 
»m«i*i  )»j*i 
i«  i  k  l  *2  )  »n  iER 
K(R1*3)«J*2 
Cm*4HSIGNQ 
MMtMMO 
CIMUUY 

ki*i«h«j 

c » K  1  ,4  I *J 
Kl R If 9l«39A0tI 
KK  «  KK  *  10 
12  CCMlhut 
RETURN 
€HC 

C  CuTFtr  KNIT  MU2IT  ION  SUBROUTINE  0UPT2 

SuCRCi  riNc  CU? T 2 
C  C  P  •*  C  N  C  I  43  10  )»  GRAPH 

tCUlVAlCNfe  fC< ?0l?»,0TCNT  1»  1C  1 316T)  iHOOUT  1. 

C  iC(20141,  IfCNT  ).  IC(  2CQU.PCNT  ). 

C  1C  I 201HI.IAPE  )*  ICI 20191. UPEND). 

C  ( C  (  2C0C I  r>f  )»  1C  ( 2C2U  .KCO^V  ). 

C  ICI 20041, PLCTNC1.  1  C  <  20091  .N0PU3T)  . 

C  tC I  2004  ),PPHT  ).  1 C  1  2023) .OPOINt 1 

0I*»F*SICN  GRAPH  12,21.7 1 m E 1 2  WOuTNOnO) 

INTEGER  PCO.T  ,  CTCNT  ,  OUlNO  ,  OPOINI 

fCulvAlCNCt  <  C  t  l945)fCUTPLTJ 
INTEGER  CuTPlT 


IC  12016) .PGCNT  I 
(C  (2015) »CPP  ) 

1 C  1  201  3  3  ,  o#;c  1 

1 C  1  202  6 1  ,  T  1  HE  1 
1C  1316BI .OJTNQ  1 
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Cl ►tKSICK  OufPlTI 15) 

KC  C-NV-0 

ITCM  •  ncc  *  1.0 

PC-M  •  T-O.CCQOOl 
PP.n  I* PCM 
PCCM  *  l 

DTCNT  »  IKOCur  ♦  A  )/5 

?f  !  ircM  .ce.  rt  co  rc*  2 

nMf{l6.4HI»';m.cnUI.UH2».CMOJ»CUM»,CII*5).CIJ*6M>1.35 
•  IC«  7  » 

6  fC#*»ATI  l  ►*  X  X  I  15,  IP  7  t  H.  n  ) 

2  1 1 pf <  |  )«T 

CfCSHT  *1 

cc  to  jm.ncpiot 

K  «  C  V.  T  P  L  T  I  J  I 

10  C4ArHU«Jla.CIKI 

Of  I  C«N 
ESC 

CuiPur  SUeaCl'TJKF  CUPM 
soe«cuTiNb  cum  3 
C <;p»o»  CJ4  3  LO )  tCaAPM 

C G'j  I  VU&'cCE  ICI  SlhSl.OuTNU  1,  ICI 3210) ,ONANC 1) .  K  I  J25b  I  •  C  2  )  • 

c  ic  i  ?e  1 7i,  orcraT  u  ici  tutkncckt  i,  ic  i*:i6J  ,p:'.nt  », 

C  f  C I  ?01«  U  I  TCNT  I*.  KC  2CC)I  ,KCH1  1,  (C(«?1*)tC  t?  1, 

C  <C.(?GCCI»f  W  KI2*6M,06R  ».  I  C  I  2Cl  *  )  •  TxPE  >. 

C  ICC  2^1^  *.  riPFNO)^  ICC  20C8)  rPLCtNC)  .  I C  <  2  CO*  I  *  V_  PlC  T  J  , 

C  f  C  C  200  5  » »  PP P  )f  i  C  l  ZOO*  I ,  PPN  T  i,  (C  .  UH£  >» 

C  ICI 202  2 l ►CPC  IN  T  I 

CCUVUEKCE  ICC  1065),  CtTPLTl 

C  I  .«  F  n  >  I  CN  P-  1  50  ‘  ,C  J  f  vil  50  I  *  ONAME  l  (  50  )  >ONAH£  2150) 
rtvP'.'.rh  r,  a  a  o  f  i  ;  •  p  i .  T  I  a  P  l  2  ) 

Cl  *  E  *iS  ICS  CLIPLTI  15) 
iSlECcR  CTCM.PCCNl  .CulNO 
i  f*  r  f  o »  o  cooler 
iMfOfc*  CUTPtr 
IF  f ITCaT.  CT,  t )  CO  TO  7 
I  I  C  N  T  1  I  f  CM  ♦  1 

*R  IT?  16,6  M  UCM  1,  Cl  IM).CI  1<2>,C<  I  ♦  3 )  #C<  I  tAl.C  C  I  *5), Cl  I  '61 .1*1  M5 

tiC , 7 > 

6  I  15#  1P7E  15.7)1 
PCCM  -  l 

7  If  ICES,  EG.  Ctfil)  CO  TO  8 
Cf«M  -  C€« 

Ml Tfl5»20)TtCER 

20  fCSMMiH  ,5)T  I«E«f  U.7*2X,  10HSTEP  S  1  U  -  lPt  i  9. 7 ) 

6  If  (T  ,U,  PCWTJCCIC15 
9  PCM  -  KM  •  f-PP 

I  F  C  »CCM  .  M.  1  >  CO  TO  3 

1  fc?;TE(5,2J  ICM*E  li  I  ).Gl»AM{2l  I  ),  I-l.NOOUT) 

2  F  rp  •#  *  r  |  1M  .3<.4M  IM,  5X,5l  7X#.2A6  1  /  /(  20X,  2A6 » 7X , 246 , 7 X ,266 , 7 X , 

12X6 ,  7X,?*6) /  ) 

FCCM  •  2*CtCM  •  4 

3  IF  C  FCCM  .Cf.  66)  CC  TO  1 
CC  A  |  «  I.KCOUT 

J  -  CO  T  *»0 1  l  l 

4  Bill  •  Cl J) 

h  a  I  T  |  U,5)  T.(9(l).  I  •  l.NOOCT) 

3  FCPt-xr  (/// 1P5E19.7/4  IA**1P5C19.7)  ) 

PCC a T  •  PCCKT  ♦  CTCM  ♦  A 
IS  IF  I  T.CT.PPNT*CR.hOPl.Or.£C.C  )R£TCRM 


PPNUPPhTAPPP 

WPCJKt  •CPOINT  *\ 

tr  I K  PC  I  NT  - SOO I  UtlJ.ll 

13  WRITE  U.U) 

14  FCA*AT  I//71P  ••••  WARNING* PLOTTING  ARRAY  MLLEO-ONLV  HAST  103  P 
CCtATS  PLCUEC  •«••,//) 

16  CPOIKT -KF'Jt  NT 
1t*£  ICPCINTUT 
CC  10  J-l.NCPLOT 
R  *CUT  Pi»T  1  j  > 

10  CBAWHICPOIN?  .JUCIMI 
18  AEToftN 
EnC 

SuSPCuTlNf  ZERO 
CO*CK  C  1 4  3 1 0 1 

Ui/lVAitsCC  4C1W8AJ  .UPtOT  I 
EQUIVALENCE  IC(  2021 1»OPO!NT  I 
iLUlVAUNtl  C  C  <  2  36  I  )  »  NCmOO  I 
EQUIVALENCE  (C42A61 )(NG$U9  ) 

ECU  I  VALENCE  IC1  3066 J .NCI  1ST } 

EwUlv  AiENCE  ICI3U7WNCCCT  l 
EQUIVALENCE  IC133?A>.LCSTAT  1 
EQUIVALENCE  Id  J3IS  J.NCSTAT  1 
ECulVAlCNCC  Id  3*.<.1  )  ,NC*sO*  ) 

EQUIVALENCE  Kl2C0b)fPcGTNQl 
I M t OCR  PtCTNQ 
INTEGER  CPU  I  NT 
LOSTAI  *  o 
NCSIAf  %  o 
NCsua  •  o 
KCPCC  •  0 
SCC1.?  ’  0 
NC*  n  O'  *  0 
KClLf  -  o 
CPC INT*0 

NPl  U  T • 0 

RETURN 

CNC 

iuescutiNe  subli 
CCP"C6  Cl  43 io » 

EQUIVALENCE  ICI2A61  UHCSUfi  1.  <C  T2A62  1 ,  SU8N0  ) 

CP’FNSiCft  lUBNCiRR) 

DC  1  I  -  t.  NQSUB 
J  •  SudNCI I > 

CC  1C  (  1.  2 »  ),  4,  3,  A*  ?»  9,  9  I.  J 

2  CALL  INPT1 
CC  TO  l 

3  CAll  OUPTI 
CC  TO  X 

A  CALC  STGI1 
CC  TO  1 
3  CAUL  C*TRl 

CC  TO  \ 

6  CALL  RNCMl 
CC  TO  l 

7  CALL  AUXA1 
CC  TC  1 

0  CALL  AU< 8 l 
CO  TO  l 
9  CALL  AUKCl 
1  (CM  !MJ( 

return 
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END 

soe»curihE  sueci 

CO»PCN  CM)10> 

€C'.jlv*LffiC€  IC(  2A61I* NQiUI  I,  4C<2442)*SU9N0  ) 
OI^eNSICN  SuBNOI?*! 

cc  i  i*i,  fiosua 

J  •  SUOKCI  1  I 

cc  rc  i  i,  a  3«  4*  4*  im  a*  9  }»  j 

2  CUL  1SPT2 
CC  TC  1 

3  CALI  CUPT2 
CO  fg  l 

♦  CALL  SfCI2 
CC  TG  l 
9  CALL  CM#! 

cc  ra  i 

4  CALL  RSC#2 
CC  TO  1 

7  CALL  AUXA2 
CO  TO  l 
fi  CALL  AUX02 
CC  TO  l 
9  CALL  AUXC2 

1  CCMIhUE 
RETUPn 

EnO 

SuesCuTIhC  SU8L1 
COP^CN  CU310) 

fC'jIVALtKCE  1C  I  7A61 UNOSUB  U  (C(2462) >SUBNO  I 
CIPPNSICN  SuBNCI 99 l 
CC  1  I  •  1 »  NOSUB 
j  -  ;  i 

CC  TO  I  1,  2.  V«  4,  3,  6 ,  7*  8,  9  U  J 

2  CALL  I n P r B 

CC  to  l 

3  CALL  OuPTJ 
CG  TO  l 

4  CALL  STCI3 
CC  TO  l 

3  CALL  C-STRB 
CC  TO  l 

6  CALL  RNC#3 
CC  TC  1 

7  CALL  AUXAJ 
CC  TC  1 

8  CALL  AUXBJ 

cc  ra  i 

9  CALL  AVJXC3 
1  CONTINUE 

RETl«N 
E  NO 

SUCCOUMNE  3TCE2 
CO*-PCN  CUJ101 

ECjIVAlPncE  4C( 201l)»KSIE#  ip  I  Cl  2020  I *LCONV  4, 

KCCNV  •  0 

LCCNV  •  0 

ASTEP  •  1 

RETURN 

CnC 

SucacuTlNf  MC.fi 
COPPCN  CUHO).  graph 

CGUlVJltftCE  AC  4  200CUT  if.  (Cl  2001 )  #  TF  >, 


1C  4  2021  1  ,kCONV  ) 


t C I  2003) »PCNT  1 
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EQUIVALENCE  J  C  «  2010  l,»T€P  )»  <C 1 201 i  )  »KSTEP  1,  1C  I  2020 )  ,LtCNV  } 

EQUIVALENCE  JCI  2021  I.KCCNW  It  «t(256il»N  I,  I  C  I  266?  1  ,H«l  N  I 

ICUIVAL5NCE  /Cl  266M,*-kAX  U  <02*64), OER  It  ICI2765),EL  I 

EQUIVALENCE  JC(20661,6u  1»  (C I  2965 > , VAR  1 

equivalence  ici  nr? i,ka*e  i»  icii9tm,nj  i,  10119751, npt  i 

CI^ENSU'N  C6RM01I  t  V  A  R  (  1011  t  (U  100) 

CIpEMICS  EUUQQ  ),  C*  APM  2,  2) 

AUXfcUB 

CALL  04 

IF  « ae  S  t  T-TF)  .16.  0.01  )  CO  IQ  20 

t  F  <  <  T  F-T )  til  .  0.1  CO  TO  10 

IF  (iCCNV  .EC.  2)  CC  fQ  20 

IF  IlCCNV  .CC.  1)  CC  fO  10 

IF  UC*  (IJ.UtO.  1CERIII—  C6«U  1*0.5 

RETURN 

10  IF  (C6iU  l  J.Cr.O.  )Ct"(  IM-DERIIMO** 

KCC.SV  ■  HCCNV  ♦  I 
tfUCCKV  *G£  •  10)  CO  TO  20 
RE lUPN 

20  PCM  -  1.0 

C  CU**  •  RECRRILOCIX  n.LOUXAII 

C  WRJT6  <6,30>  X3»  X4 

C  30  FCXRAf  (38-0  RESTART  INITIALIZERS,  Xj  AND  X4,  ARE  2FJ1.0) 

IF  I STEP  WCO.il  .  1  CO  T □  AO 
PRfcCt*  n  CfcRI 1  I 
OEP(l)  a  0. 

NJ*N-1 

NPT*0 

CHI  AHRMAUXSUB1 
otRm  PR6ce« 

40  CALL  LUPT3 

xst:?  *  2 

RETURN 

EM) 

S'JSPCor  INC  APR  A  !  AUX  SUB  ) 

C  •<  I  NG  L  E  PRECISION  V  6  R  S  I  CN  •  INDEPENDENT  VARIABLE  IN  OOU0LE  PRECISION 
CO'-vCN  CI4310) ,C«aph 
C  CCune  PRECISION  CHt.TNE 

C  CCuClE  PRECISION  NfwCt 200), NEwPI 200), OL0I200) 

CtuPlE  PRfClSIOh  Cfi.T,TPE 

CCueie  PRCCI6ION  NE*C< 200),NEwP( 200 ) ,0L0(2001 

01  pens icn  r 1 1000 ) 

OIpCnSICn  0(1011  *EUIC0»  tCUflOO) 

OIPESSICN  VI 101  ),GRAPh(2,2) 

C  OIpEnSICNRCmPI  200 ) *  CC  01  200  )  ,N£ WC I  2 00 > 

t^UlVBLtNCE  <C (  2662 ) , HN IN  I,  I  C  <  2  A  0  )  )  ,  H*  A  X  1,  (0(2664), 0  I 

ECjWAcENCC  f  C  I  2 76 5  )  *.€L  l»  (C120651,Eu  I,  (C(2465),V  I 

equivalence  i c 1 1 o / 1 1 , r 1 1 e  u  uund 

EQUIVALENCE  C  C  <  19?  1  l,KA$E  )#  <0(19741, NJ  It  4  C  I  19751  ,NPT  > 

Cat#  rCukt/c/ 

CATA  Pl,P2fP3,  P 4/2. 29 16 6667, 2. 4583  3  3  33* 1,541666? *0. 375/ 
data  c 2. C3.C4/0.79 1666 7, 0,2063 31 3,0.041666?/ 

IF iKASt -CT.0 ICO  IQ  20 

M  »NJ 

J2  •  J 1 *N  l 

J  3  •  J2  ♦  Nl 

J4« J3#Nl 

J5- J4»M 

J6 •  J*  •  N  I 

J?-J6«M 
J8*J7*M 
J9  •  Jfl • Nl 
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K*St»KASe  ♦! 

C*NPT«E5.0  ACAM$ -POOL  Tos  INTEGRATION  node 
C»NPT.EC.l  PU*C€-*UTTA  INTEGRATION  *GOE 

C»NPT.EC.2  BEGINNING  A C A P S-PCUL TON  k  I  TM  AUNGE-KUTTA  START 
20  IF  lf.PT  .EG.l  I  GO  TO  AO 
IFUPT.EC.2 ) GO  TO  30 
|F{  SNGLIOELI  ).NE<(O.S*Cm  )  ICO  TO  30 
IFUCUM.LT. 3)G0  TC  40 
GG  TO  2C0 
30  KGUNT-0 
NR  T  «0 

(•START  RUKGE-KUtTA  INTEGRATION 
(•CQPPuIE  Kl 

40  OC  SO  Hl.Nl 
KEhFI  (  )  V4i  HI) 

50  CCNT1.NUE 
TH£«VU  I 
KCuNT»KCuNT*i 
DC  60  I  '  1  iNl 

60  cum  »ciiiioi  nn 
C  CCPPUTE  *2 

061  T»O.S«C( 1  I 
T»E«  T  n  t  »C£u  T 
VI  l  JkTPI 

CC  70  MI.U 
If  ULOKT*KE.2)G0  TC  65 
HI- J9*  I 
TH1  HNE*F  I  I  1 
65  TIII*Clhll 

NEAP  III  v’.EwPlI  H.O,S«OLCI  I J 
70  VI  I  *  l  HMWMI  I J 
cal:  AL  XSL'3 
CC  CO  K  •  l  •  N  1 
80  N 6 w C I  1  )  •CIU*01  1*1) 

c  ccppufe  Ki 

C  C  9  3  I  •  t , N  1 

Nek  Ml  HNEkPtl  HO  «5»1  NENC1  I  J-OLOII1  i 
90  VII  *1  I  “SE**P  I  I ) 

C*LL  A'jtlOB 

cc  i:o  ui.ni 

*  2  •  J  7  *  » 

100  IU?HCMI«CI1*U 
C  COPPLTE  PA 

TPt-TPE «CELT 
V I  l )  ■  T  p  E 
CC  110  Hl.Nl 
K2  «  J  7 ♦  | 

I  )*,-ErfPI  1  )  *T(K2  l-0«5*NEtfCC  1 1 
110  V 1 1  ♦  l  I  •NEuPM) 

CALL  A'-'TSUa 
CC  120  HUM 
K3-  J8  ♦  I 

120  TlK)l»Cin«CU*U 
C  CQPPUU  value  of  function 
CC  1)0  IM.M 
R2-J7*  I 
R)“ J8  ♦  I 

HEMP!  I  I  •NEWpfl  I -T I *2  1*0. 166 6666  7* 

XI  CL  Cl  1  I  AN  EM  Cl  I  MNEkC  (  1 )  at  1X2  HUM  )  MIX  31) 

110  VI  ni  HKEnPI  I) 

140  CALL  AUXSU8 

00  ISO  IMtNl 


I 


t 

! 

| 


RS*J1*  I 
K0«J2*I 

*2- J4*  I 
*J» J5*  I 

*<.•  J6  ♦  I 

TlKU«nK)l 
T I K3  I  *  T  (K2I 
TU2)«TUll 
T<kD«T(KOI 
T(KOI*TUM 

t  i  rs  t • r  i  i i 

ISO  T( I  1 • c  f 1*1) 

RETURN 

C  AOAP»S-hCuI  TCS  INTEGRATION 
200  KCu.M-kOAMI 
0EuT»O.5*Cin 
CC  210  IM.NI 
RL*J2«1 
K2«  J  3 ♦  1 
K4- Jl ♦  I 

C  CO'PuTfc  Y-PMfClCf EQ 
CLC  I  I  )  1 1 

NenPl  nxCLLM  I ltCll)«IPl«1l  I  T-P  2«T<R4 1^3*1  tKU-M*Tllc21 » 

210  VI  I *1 1 »NE*P( n 
T*E*Th£<CIL1 

vm«rce 

CALI  AUJlSoa 

rs-o 

cc  2so  im,m 

R2- J2»  I 

K4- Jit  I 

C  cuppuit  t-cc*<«tcrEO 

NE»C I l  )  *CLO I  n*C<  |l*tP4«CI  HlMC2«rm-C3»TU4UC4tf(K2»  ) 

IPlt-NlN.tu.HHAR  JGO  TO  2 SO 

TtPP-AeStSNCClNtvCI  II-NEWPUJM 

IE1TEhP.LT. EDI n )CC  TO  24 o 

IF (q ITE.iE.O.OIGO  TO  230 

SOT  IHc  »SvGL  1  THE  ) 

W«!7E16,220I I.SPTlPE,TfHP 

220  FC®haT11>-  .UHSTATE  VAR  l  AQL E *! 3 *26H  6XCCE0E0  TOLERANCE  EAftO*k 
X7h  T  in- »F  1*.7,9H  TEHP»  .1PEIT.B) 

2 30  IF  l  ABS  (Sn.UCEI  T  n  .GE.HHIN1C0  TO  270 
240  1 F  f TtHP.LT.EU  I  )  IKS«KS» 1 
2S0  CCM  INUC 

If <K$. IT. M ICO  TO  290 
in  A8S  «  C 1  *.C|  1  )  1  .CT.HHAX1G0  TO  290 

C  SET-UP  FCR  CCuflLlNG  STEP  SUE 
r ( RCUVT.L6.61GC  TC  290 
CC  260  l«l,Nl 
K1 • J 1 M 
K2-J2M 
KJ*J3M 
RS- J5*  I 
t 1 1 1-iiKD 
tiki  1  •  T  I R  3 1 
260  TU21-1USI 

ci  n-c(  i)  tom 

KCUM-4 

CElT*0.S*C1 11 
GC  TO  290 

c  SET-UP  FOR  HALVING  STEP  SIZE 
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JIO  if  mCuM.LE.%1CO  TO  110 
Uf-TKE-CU  l 
VI  1 l-TfE 
0(1  )*C6lT 

o t l  mo  ,s*di  1 1 

CC  280  1*1. M 

Kl»Jl •  I 
A2- J2  •  I 

m- ji*  i 

Af.el 1  I  »CLO 1  II 

V(  l«t  >*CLC(  11  .... 

n,>|.T(,!|.0.5,ITml»-l(Ull 
T(«!hTI«ll 

180  I  (  «  1  l«t(K  1 1  «0.5*t  T  I  1  l-UM  >  1 
KCC»f •* 

c  |MEC«*nCN°*S.  f  INl&hEC.  SET  UO  DERIVATIVES  RhO  E*IT. 
2^0  CC  100  i»UM 

NgwPi  I  I  «h(aiC <  l  l 
100  vihl  )  t  * 1 

C  RETURN  !c  i>C  PRECE CIHC  R01M  RESTART  RK 

JIO  OC  320  1*1. •*» 

M'JVM 
KE»M  1  >*T|«  l» 

110  VI  l  *1  •  "TUI  1 

TRE»T  RE-A.O *CI 1 1 

viimm 

oi  i  mccii 

CALL  AUXSUE 

CC  To  30 

c  a.  r* 


S^PPCUTINE  *U«t 

ECumiuJa'KUTMUNOROO  ».  4CI2J41I.XKOONOI.  ICIH4U.N  > 

OmehSICH  XP0CNCIV9) 
h  •  l 

CC  1  I  •UNO^OO 

1  ,2*115.14. 17,  18.H.  10.  31,  31, 13. 1*.  35.  >4' 1 71 'C 

2  T  ACL  All 
CC  TO  1 

3  CALL  *11 
CC  10  l 


CO  10  1 
5  CALL  *M 
CC  T  C  l 
4  CALC  »M 
CC  TO  1 

7  CALL  CU 
CC  TO  1 

8  CALC  Cl  1 
CC  TO  1 

V  CALC  CM 
CC  TO  l 


10  CALC  C*l 
CC  TO  1 

11  CALC  C51 
CO  TO  1 

12  CALC  C4i 
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9 


CC  TO  i 

11  cut  CM 
cc  ?o  i 

14  cut  cal 

CC  TO  l 

15  CUl  C91 
GC  TC  I 

16  CUl  CIO! 

CC  TC  l 

1?  CUl  Cll 
CC  TC  l 
la  CUL  021 
CC  IG  1 

19  CUL  C  3  I 
CC  IC  1 

20  CUl  C 4  I 
CC  TO  1 

21  CUl  CM 
CC  TO  l 

22  C *t l  Cll 
CC  TO  1 

23  CUL  G2I 
CC  TC  l 

24  CALL  G 3 1 
CC  TO  l 

25  CAU  G4t 

C  TO  1 

26  CAU  G5 ! 

CC  TO  1 

27  CAU  G&I 
CC  TC  1 

Cam  si  I 

CC  TC  1 
29  Cau  S 2  I 
CC  TO  l 
JO  CAU  Sll 

CC  TC  I 

31  CAU  S 4  I 
CC  TC  1 

32  CAU  S5I 
CC  TO  l 

33  CAU  S6I 
CC  TC  1 

34  CAU  S 7  I 
CC  TO  l 

35  CAU  sei 
CC  TO  l 

36  CAU  $9  l 
CC  to  l 

37  CAU  SIOI 
i  co?iMje 

RETURN 

£N0 

subroutine  Auxsua 

CG*Kn  C  14310) 

EQUIVALENCE  ICI2000I.T  I.  IC I  2)61 ) , NOHOD  ), 

ECtlvAlcNCf  fCC  2361  It  ■!  |.  (CtZ*62),|Pl  ). 

EQUIVALENCE  IC(296)I,VA«  > 

EQUIVALENCE  ICI202CI.LC0NV) 

OIPEKSICN  CERIJOU  ,  ¥6*ll0l> 

ci'Cascccc  xMOCNomt 

CC  50  I  *  2.  N 


IC 12)621 .XNCOHOI 
IC I  2666  I »0Eft  I 

IPII100I 
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J  •  IPUI-ll 
50  C  c  J  ♦  J  >  «  v*ft  f 1 1 
T  «  VASILI 
CC  1  ]-l.hO*CO 

IFllCCNW.eO.imTufiN 
i  •jpccsc( 1 1 

CC  1C  l 
3  CALI  M2 
CC  10  l 
*  call  M J 

cc  ro  i 

5  OLL  AA 

cc  ro  i 

6  CALL  A5 
CC  TO  1 

T  CAU  Cl 
CC  TO  l 
a  CALL  C 2 
CC  TO  1 
9  CAU  Cl 
CC  TO  l 
10  CALL  C* 

CC  TO  l 
U  CALL  C  5 
CC  TO  l 

12  CALL  C6 
CC  TO  l 

13  CALL  Cl 


vC  t  U  l 
U  CALL  ca 
CC  TC  l 
15  CALL  C9 
CC  TO  l 
lt>  CAU  CIO 
CO  TC  x 
l?  CALL  Cl 
CC  10  1 
l«  CALL  C2 
CC  TC  l 

19  CALL  CJ 
CC  TO  l 

20  Call  os 
CC  TO  l 

21  CALL  C5 
cc  ro  i 

22  CALL  Gl 
CC  re  i 

23  CALL  C 2 
CC  TO  1 

2*  CALL  G  3 
CC  TO  l 

25  CALL  C A 
CC  TO  l 

26  CALL  G5 
GC  TO  l 

27  CALL  G6 
GO  TQ  l 

2«  CALL  si 
CO  TO  1 
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29  cm  $2 
CO  TO  1 
JO  CALL  S3 
CC  TO  1 

n  c  *  1 1  s* 

CO  TO  1 
J2  CALL  S5 
CC  TC  1 

33  CALL  S6 
GO  TO  l 

34  CALL  S 7 
CC  TO  l 

35  CALL  S 8 
GC  TO  l 

36  CALL  5 9 
CC  TC  l 

37  CALL  S I 0  . 

I  CC**  T  1  SUE 

00  60  l  *  2»  N 
J  -  IPUI-IJ 
60  OER  III  4  Cl  J) 

RE  TORN 
ENC 

SueacGTinc  RCSEf 
Cl>*C.\  Cl  63  10) 

ECOTvAL€KCElCUJiO)*KH),<CtM 

ECUlVAtENCE  CC  I  3O66UN0L1STU  ICO067)*tlSlNft|,  (CMUTI. VALUE  > 
CCUlVAlCNCC  I  Cl  2900  1,NORN0H  J  *  <  C I 2441 1 iRNOHNO) 

FCUlVAlfNCE  1C  I  2  800  I  »  8  ) 

C !  ►  £  N  S  l  C  N  *14310) 

CIKFSWCN  USTNCI50)  *  VA*  «'c  1 50) 

r.i  bhK<,  ten  a <  ino» 

OIVCNSICN  RKCHKCI50) 

CC33I-1.3509 

JJ  cm*c. 

JC»  3*09^* 

OL«JI«J5IO»JC»2 
43  MMIMU-KU) 

K|  2461 )  ■  *  |  4309 1 
M  236 l  )•* (6308  ) 

M  3066  l*KI430n 
M  3  16  7  I  »M  4  306 ) 

Kf  3  3  38  )  *K I  A  305 ) 

Kt  33  )9)»K{4204) 

IT  [  1984  1-K14303) 

K  l  J4sO  )»M430 2  ) 

K| 20 W  l«C42010) 

K I  1975  )«2 

RETURN 

END 

C  OU^Y  SUBRCUI |  ME 

SueRCUflNC  CUKNV 
C  ENTRY  A  1 

ENTRY  All 
C  ENTRY  A2 

ENTRY  A  2  f 

C  ENTRY  A3 

C  ENTRY  A)| 

ENTRY  14 
ENTRY  A4( 

ENTPf  A3 
ENTRY  A5I 
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c  ffiTRV  Cl 

c  entry  cu 

ENTRY  C2 

ENTRY  CM 
ENTRY  C  3 
ENTRY  C  3  t 
C  ENTRY  C* 

ENTRY  CM 
ENTRY  C$ 
ENTRY  C  5  I 
ENTRY  C6 
ENTRY  C 6 l 
ENTRY  CT 
ENTRY  CM 
ENTRY  C8 
ENTRY  C8l 
ENTRY  C9 
ENTRY  CM 
ENTRY  CIO 
ENTRY  CIOI 
C  ENTRY  Cl 

C  ENTRY  Cl! 

C  ENTRY  C2 

C  ENTRY  C 2  I 

ENTRY  Cl 

entry  cm 
entry  C* 
entry  cm 
ENTRY  cs 
ENTRY  CSI 
ENTRY  Cl 
f  N  1  BY  C  1  I 
C  ENTRY  C2 

ENTRY  C 2  I 
C  ENTRY  Cl 

ENTRY  CM 
C  ENTRY  G5 

ENTRY  CM 
C  ENTRY  G5 

ENTRY  CM 
ENTRY  C6 

Entry  cm 
c  ENTRY  si 

C  ENTRY  S  1 1 

ENTRY  S? 
ENTRY  S2I 
ENTRY  S3 
Entry  $31 
ENTRY  $R 
ENTRY  SM 
ENTRY  S5 
ENTRY  $51 
ENTRY  S6 
ENTRY  $6t 
ENTRY  S  7 

entry  sti 

ENTRY  S8 
ENTRY  S8I 
ENTRY  S9 
ENTRY  S91 
ENTRY  SlO 
ENTRY  SIOI 
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fh?>V  AU*Al 
E^T^Y  AviA2 
EhMY  1^11) 
fMoy  ,11.x  3  l 

*v,xe2 

(M«1  A'.*  21 
Fs'XY  K'Cl 
f.’^Y  Ai_*C2 

tNNY  4  ■„  A  (  3 
f»Mr  C*>T«l 
t  N  T  ft  Y  :Sf42 
FMrty  C.sMl 
tMdr  is  *  \ 
f  t  ^  Y  I  •> 

£  S  T  «  Y  l  V  ; 

f  M  nT  i  Ui»  T  I 

C  fsr»Y  Ct>T2 

C  E.sr^Y  KPtl 

cS'^Y  PHCCtS 

c  f* r «y  *e st t 

ts  f  - Y  «SCm1 

t  M  *€  Y  P  S  C  •*  2 

£st^y  asr*j 

C  s r 4  T  «'Cfl 

C  ts  r  A  Y  SICs? 

C  E‘.ny  s  l  C  t  J 

C  £N?hY  SLfLl 

C  fst-v  SVC -2 

£Mq»  ».|«ur 
£n*=t  CCcsf* 

«t  lefts 
«  *.  ■• 

C*4  K*  USE  blTH  C00l«2.  KNJ,  KN) 

C 

5L£-CLM».£  TERROR  <ll*0Et> 

( C-*Cs  » ( *310) 

-IVAUKCE  KI  202C1.1CCNV) 

*  I  T  £  is,  13  I  xiaBEL 

10  r  i)r  ••  A !  <  s)HO  fcQ  Af*0  JOINTS  S'KIMCD  *0*  arc  ,  9*4 

C  7nla?i€  IA6  ) 

cm  €  *  i  r 

C 

EyIPt  aE««0* 

*»[TC  XLABfV 

20  FC-*AT  &  *3hQ  0V»T  OB  A«R0  TA6L6  ARGLHENT  ARRAY  .  Bit 

C  ?hIaCi€  .J6  ) 

C:  AO  2*l2?JtlZ>l.T 

40  bi.lK6.iJ)  Ct  I  ),C<IUUCf  («2l»CMOt*C(|#4»tCn*9l.CU»6> 

10  K  *•  «*  at  I  lb  , ?t|*. 7 1 

b4|I(  u.!0»  Ct2C00)*U)67)»Ct34«ltCI2C4)9CI)44).Cn70)»C«lllTI» 
C  C«  IU4I. Cl  HHI.CI  11201 
ICE  XV*. » 

At  T v«H 
ISL 

n v;?k*  jikc  id 
si‘-c»  bin  «  i/)T.2*>) r« i 
**  !  i.s*k 

esc 

l-.M'flfls  coic  ID 

CCS*-  iX/57. 295731 

81  I  L  AS 

tsC 
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r.  rsrsnr\nr>rsnr%r\nr* 


fukctic*'  »>**c  tt.xi 

lf(l.fC.O..ANC«Y'EE*0*EVCICl 

XtiNG*  M.2XB  7fl*M4N2  I1.*l 
CChM  NUt 
•E<Of X 

Su5»CI_riN«  TXOI .E  II.XI.Tt,NX,X«,XC»8£L.YI 
Cl'ENSUN  JLi-’.tl  I?)  ....  „ 

y  .  CCC l x  2  (  A.xl.Y  1  .  NX.XX.XUBU ) 

BET  OSN 

l^,c .  ,»ol2<«.t.«*i.i u,><x.x,'*te*.xn9et..«) 

XLietL  » 

fit  »U*N 

luCXCLtlSS  ilJtX.Y.2.»»IUXI.Nxr2.X!lXlt»,Xl.*tt(..Xil 

OIXENSU*  **  t  («l.MYfl*].fcll2>.*>-^Bet«2* 

«o,. 

^XTZ  12  1  l 

P€  TURN 
tsC 

S'jCPCUtlfcE  tipevickti 

ft  £  T  •*  f*  h 

Su^BCbflht  i»PIT6l 

M tuft* 

So?Ct-LT  I n E  PLOT  4CG**«>>‘»0<>OIHT,VLA»L€wT  l«6  ,*.?U0T*I 

*(  1  L  «H 

CNC 

Sue«CtT!K«  PtOT2t*PLOT2) 

»t  TU«n 

6NO 

5ue«»CLT|KE  PlCTMUCPLOT  ) 

INC  . 

,  . . . . . .  *  * 

SuB»CutlM  C CO  1*2 

™To“.t  *  SM  ct  rci-.u  wn»  »  COMMAS  f  UNCI  ICS  ih»i 

4l.0l.XUS  *  Mt'.Cl-  C v«vl  lY»f  c,|1,vE  MT> 


tOO,.2  l  X  .  XI  ,  »!  .  N  .  f  •  ' 


y  x  COOI»1  l  x  «  xl  ,  xl 


f  ,  XEABEl  » 


nt  *.  C  *  l  P  T  l  CN  C*  P  *R  A**£  TC  *  5 

.  AOfvj^tsr  -  iNiePt^CN?  VHUCU 

.  ABflAY  Of  IM-tPtNtfXT  VUIWU  .  * 

t  *  r  f  £  <  j  ( *<  *.  t  N  f  V  t  A  ’ll  €  •  ^  w  _ 

i^SfB  Of  POINTS  •tPPtSf.UO  ftt  H  *N0  Tl  A1AAYS 

«  (hT(PPCL  A 1  lljN  CC’ilROL 

1  ex'.  1H-X  /tao  -  S;»*IONX  ClNt  IhUHill'.x" 

»cim*i  -  enc  nii'i'i  iNit*pot*T icx 

0.0  ST»*IGhT  line 

1.0  Evict  EX»*#OllC 
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C  X  l  AS  £  l  HOLLERITH  FIELD  OF  UP  TO  6  CHARACTERS 

C 

C  REMARKS 

C  EXTRAPOLATION  IS  CONE  by  PASSING  A  STRAIGHT  line  thru  the 

C  TtaC  POINTS  AT  THE  ESC  INTERVAL* 

C  THE  ARRAY  CF  THE  KOEPEkCENT  VARU8LE  »  XI  »  MAY  BE  IN 

C  EITHER  increasing  CR  DECREASING  OROCR. 

c 

C  •  . . . 


FUNCTION  CO  C  I M  2  I  X  »  XI  »  Y1  ,  N  ,  f  t  XLAQEL  I 


OIPENSICN  X 1  (  N  )  ,  YIINI  *  P  C  2  »  »  E  1 2 1  »  ISU.21  *  XlACEH2) 
logical  cur 

CATi  IS  /  -W  U.  -2*  -If  0.  W  -I.  0  / 
sqo  out  ■  .fav.se. 

Nl  •  N 
XX  *  X 

J  •  l 

IF  T  M  -  2  »  190  .  1 2CG  •  300 

ISO  CALI  TERROR  IXLABEL I 


200 

CCCl* 

2  • 

Y  1 1  J  } 

RETURN 

300 

KPL 

•  1 

RFU 

•  3 

If  1 

Aim 

-  x  1 12  » 

)  ACO  ,  150  •  600 

400 

CC 

500 

J  •  l  « 

NT 

If  I 

XX  - 

X  II J  >  » 

9C0  4  2C0  ►  500 

500 

CCM  tsi:f 

GC  TC 

POO 

400 

CC 

700 

J  -  1  ♦ 

Nl 

If  < 

*  I  (  j  ) 

-  XX  ) 

900  ,  2C0  •.  700 

700 

CCMlNUE 

600 

J  * 

M 

CUi  AERRCR  IXLABEL  I 
GC  TO  1300 


900  CUT  -  f  tlT.  0.0 

IF  I  J  -  2  1  1200  ,  10CO  •  X ICO 

1000  XPL  -  2 

GC  TO  1300 

1100  IF  f  j  -  h|  I  1500  ft  1 ACO  »  1300 

1200  J  «  2 

1300  CL T  •  .TRUE. 

1400  *PU  •  l 

1500  At  •  I  XX  -  XIIJ-11  >  /  I  XIIJ1  -  X  •  I  >-l  I  1 

CCCIM2  •  AL  •  YIIJJ  ♦  (  1.0  -  AL  J  •  YltJ-l  ) 
IF  1  CUT  »  RETURN 

CC  lflOQ  HP  ••  R PL  •  KPU 

PIRPI  ■»  0.0 

CO  1600  *  •  l  ♦  J 

JO  •  J  ♦  XF  ♦  X  -  4 
If  I JOI  1900 •  1550, 1900 
1850  XC-O. 

YC-O. 

GC  T  C 1950 
1900  K-iUJOJ 


no 


r 


TC*T1 ( JO  » 

050  Jl'JHSm.HPJ 

J2  *  J  •  ISt**l**P» 

1*00  firm  *  piKp)  ♦  vii  •  (  u  *  nun  )  /  <  io  -  nun  i 
•  •  I  XX  -  XIIJZ)  1  /  I  no  -  X | 4  J  2 1  ) 

if  I  KPl  .NE*  KPU  »  CO  10  1700 

J1  »  3  -  *PL 

Fiji)  »  ccci*2  ♦  f  •  i  Pt<F)  -  coomz  > 

£1411  *  *ei  1  FIJI)  -  COJlnZ  1 

ITOO  El*P}  •  A?S  I  FtRFl  -  C0CIM2  1 

1100  CCKTlMjC 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

in 

no 

1)2 

c 

c 

m 

m 

140 


IF  I  (111  ♦  cm  .EC.  0.0  1  RETURN 

CCCIHJ  •  l  I  €11)  •  H  )  •  Pm  ♦  I  £12)  *  (  1.0  -  At  i  ) 

•  •Pin  I  /  I  l  fill  «  H  )*16W)»11.0-Al.))l 

MTbftN 

€N0 

2-C  IhENSICNAL  INTER  PCL  AT 1CN  SUQPftOGAAH . . . . FCC0N2 


C At l I NG  &CUUENCE  - 

l  ■  f  CCC».2U.y,jh,yuz  i,Nxc^Nr,sx^xx*xtAacu 


x  ■ 

Y  • 

xt  ■ 

Y  l  ■ 

zi  ■ 

Mf  ■ 

XT  • 
AX  • 
IK  • 

meet.  • 


ARGUMENT  -  1ST  WAFUOlC 
nCl.HCNT  -  2NJ  vaKlABU 

A  •  A  Y  Of  1ST  VAR  IAPLC 
ARRAY  Of  2*sC  YARlAfite 
ARRAY  Of  CfPfUfM  VARIABLE 

ci-cxsicNtr.  si/e  of  xi  array 

C  f  VftuCl  IK  ARRAY  Y I 
CP  VALUES  is  ARRAY  Xl 

£Sc  INTERVAL  1NT£RP0lAMUN  CONTROL  CONSTANT 
►•CLtCMfH  flElC  CP  LP  lo  6  CHAfiACItftC 


THIS  AC'JriNC  ClfffRS  fRO*  FCCO*Z  Is  ImaT  TmE  21  ARRAY  DOES  SCT 
►  AV  6  TO  BE  PAOcC  -  1  .£ .  •  M  CUC  S  NCI  HAVE  To  OCCUPY  CON- 
SlCUlvC  tCCAMCSS  1ft  CORE,  A  NO  IN  T  >«4  T  EITHER  OR  OOTH  Th£ 
XI  ANC  Y |  ARRAYS  nay  BE  IN  ASCENDING  Oft  OESCENOING  orcch. 


FONCTICN  f  CCGNZ  I  X*  Y*X1  ,T|  P.1  UNXO.NY.NX,  XK  ,  Xl ABEL  , 


OlfENSICN  Xllllt  Ylin,  II(NX0»1U  114*  »  XIA&EK2I 


IF  (NY.Cr.4)  GO  TO  120 
*0-1 

M  IS  TH&  INDEX  rru*9C«  CP  ThE  FIRST  Y  CLRVc  TO  BE  USED 

h*»XY 

NV  !  9  THfi  COUNT  OF  THE  NLFI8EA  OF  Y  CC*V£J  TO  9£  CSEO 
CO  TO  200 


NA  ■  A 

IP  *V1<U-VII2>  I  1)0»  ISC*  133 
CO  112  RxUNH 
If  i  v-r  lull  15Cil50»n2 

ccm  isue 

co  re  no 


CC  134  X«UNY 
IF  lYlUlot)  190*  HO#  l>4 
ccnt inue 
M*AM 


111 


GO  10  200 


!50  IP  ln-31  15SU55tUO 
115  M-l 

CC  to  200 

160  IF  IK-NTJ  U9.M0.140 
165  N3-X-2 

700  L*N3 

CO  300  M1.N4 

T  I  I  MCC01N2IX.XU2  II  Ut  UNXtXK.XLAftEU 
300  L*LU 

FCLCN2*CCCIP2(  V  »Y  J  I  N  3  )»T,NA*XK,XLA5£L  I 

RETURN 

CNC 

3-CMtNSI  CNAL  INTERPOLATION  SuBP AOCR  AM  ... «  F  COON  3 
CALLING  SEQUENCE  - 

«  •  f  C  CC  N  3  I  X i V . 2  *  X  I • Y  t  •  £  1 I «  N2  »NY  t  NX  »  XK . K  L  4  BE l ) 

X  •  ARGUMENT  -  1ST  VARIABLE 
V  -  APGu-Cnt  •  Jn 9  VARIABLE 
7  •  APCu'fNT  -  3«C  VARIABLE 
M  •  aRRAy  OB  1ST  VARIABLE 
tl  •  array  of  2k C  VARIAOLE 
21  -  ARRAY  OP  3ft  C  VARUBLS 

m  amray  of  ceoeNCcfT  variable 
fci  *  NL-'tM  Of  PC|NT$  IN  21  ARRAY 

K  V  •  KifRE*  of  points  in  YI  array 

kx  -  number  of  pcints  in  xi  rrav 

x*  ■  ENC  interval  INTER^ULAI  .on  CONTROL  CONSTANT  10.0  10  U0> 

XLABEL  •  HOLLERITH  field  of  up  TO  6  characters 

f  CUCN  3  CIEFERS  FROM  YCUUM  IN  THAT  THE  Wl  ARRAY  OCES  NC  T  NEEO 
TO  EE  RACKEC,  . .  E  •  •  wl  NEE  C  NOT  OCCtPr  COnSECuTIwE  LOCATIONS 
IN  CORE*  A.NC  ANY  OR  ALL  ARRAYS  NAY  0£  IN  EITHER  ASCENDING  OR 
DESCENDING  ORCER. 

FUNCTICN  FCCCN3IX,Y,2,X|*YU2I  *Nl  ,N2*NYfNX, XK.XLA5ELI 

oipessics  xiii>»  y i ( ii»  mu,  hmi,i*iu  iui  ,  kla»eli2> 

IF  (Ni.ct.4l  GO  TG  120 

M*l 

B5-S2 

GC  TO  200 

120  MU 

IP  mill-2  1  12  )  )  130*150*133 
130  CC  132  lUl.Ni 

If  I  2-21 (K I  I  150,  1 50* 132 
132  CONTINUE 

GC  1C  MO 

113  CO  134  *•!.*£ 

IP  (2IIKI-2)  150,190,134 
134  CONTINUE 

|40  M*N|-J 

CC  TO  200 


L2 


ISO  It  ( K  •  C  T  •  3 )  CO  TO  MO 
M«l 

CC  TO  200 
C 

1*0  It  ( *.  C6.NZ  1  CO  10  MO 
*4*A-2 
C 

200  14*4 

CC  300  MifNS 

•  UM  FC0CK2  U,T,XUTl,XM  U  1  »i.l»H**NYpNX,Xt,XLA0EL» 
100  l*L«L 

C 

*CCCh3-CC01«2l  ZvZ  I  \  N4  ) ,  T  rNZ  *XK  tUAQEU 
f>E  ru«N 
MO 


COPT  2,3 

3 

GIGS  2,3 

A 

G2-T 

23 

G3-H 

2* 

G5-H 

26 

Sl-T 

2  B 

Cl-T 

7 

C4-T 

10 

Al-T 

2 

A3-T 

A 

A2-H 

3 

Ot-H 

IT 

C-10  H 

16 

CZ-H 

18 

TF 

2001 

7  .0 

tec 

\  3 

/  . 

HP  I  H 

2662 

»CG  75 

HKAX 

2*63 

.COM 

Cf 

26ft  A 

.C02S 

c  p  r  r  a 

Aft  3 

1. 

c*cr 

A  ft  ft 

.05 

GS* 

450 

1. 

R0M.C* 

453 

UO. 

t*S*R 

A  5  6 

2. 

V/t 

16U 

t  10. 

V«fc 

1603 

44  0  . 

solve 

1666 

-20. 

RSI thl 

166  7 

9CC0. 

204 

.82765 

TH£ |*-C 

42  7 

.  1 

QAIPHA 

367 

2.3 

CB1  AS 

852 

1. 

HI  MC 

850 

15. 

HLMC 

851 

15. 

GV 

856 

6.0 

G  / 

853 

3.0 

TAUT 

B  6  4 

6 .0 

TAU2 

86  3 

6.0 

TAUt 

677 

•  5 

BC*  A  X 

1121 

20. 

At A*€* 

1  306 

•  267 

AttCTH 

1307 

.564 

Af  xCG 

1313 

2.M 

At  ZCG 

1315 

-.0063: 

ACiCC 

1421 

2.54 

AA4Y0 

1420 

15.4 

113 


3  RQCC# 

1418 

3  ROCGC 

1417 

3  0*P 

1416 

3  OUT 

1415 

3  CISP 

1414 

3  VAIRSrO 

207 

3  C*SK2 

668 

3  CSAlGK 

1403 

3  ft  A 1  L 

1317 

3  ftlUG 

1316 

3  CEFCC3 

466 

3  CETfliO 

469 

3  CKACT 

1163 

3  T0Y2 

866 

3  TCLV2 

866 

3  0*5*1 

472 

3  VI A  2  ft  ft 

473 

3  C*S*2 

668 

3  C*fti 

2664 

3  f"l*C 

1419 

3  G 1 

855 

3  ACftAV 

1627 

3  C  ft  AC 

1751 

3  SUP 

2010 

3  ec«6 

OOCO  1 

3  fiAOtuS 

or, 3 0  3 

3  vanO 

OCC02 

3  A  If  AC 

11 

3  MIL 

7 

3  SPLf*C 

6 

3  CPP 

2015 

A k 

a 

1739 

4 

PCU  FLAP 

1232 

4 

phi  fismi 

352 

4 

VX  £  Aft  TH 

1603 

4 

X  rsi 

1615 

4 

c 

1743 

4 

PITCH  FLAP 

1233 

4 

T*£T»  MSSIE  350 

4 

vv  EARTH 

160  1 

4 

v  FSL 

1619 

4 

ft 

1747 

4 

Vl»  FLAP 

1234 

4 

P3I  •’ISSUE 

351 

4 

vl  EARTH 

1611 

4 

l  HSl 

1623 

4 

ALPHA  OCG 

367 

4 

SETA  CCG 

368 

4 

ALPHA  PftIPC 

369 

4 

PhI  PftIRE 

370 

4 

A2  6  A 

1626 

4 

P  ACH 

204 

4 

AIR  SPEEO 

207 

4 

POVhPC 

203 

4 

MHRST 

1410 

4 

THETA  LCi-V 

363 

4 

CPS  l 

435 

4 

TMI*  CIKtll  *27 

4 

C  CCH 

869 

4 

P 3 1  LCS-V 

364 

4 

IPS  V 

436 

4 

ft  3  I  G  I  *  8  AC 

431 

-.1165 

•  •3165 
15.6 
133, 

153.8 

880. 

6. 

1. 

7.5 

2.37 

3. 

7.3 
i  . 

1. 

5.0 
1.0 
l  .CO 

•  6 

•  CC75 

•  216 
6.0 

32.176 

57.295778  -0. 

2.  -0. 

0. IQ  1*0 

0.250  1.0 

0.20  1.0 
0.0 
1.0 

1.0  1.0 

l.c 


114 


»  CC**  »Tfl 

CELT  8  COH  B?S 

CAHAAK  357 

G*l«n*V  158 

<18*  1625 

H  AC  J 

81  AO  7 

*T  AO  5 

M  AO* 

-0 

ASLANT  1*67 

C88  2015 

SOIVf  1*66 

-0 

Cl  8  55 

C»  856 

a 


-c. 

scoo.o 
.  1 

-10. 

-0. 

5.0 

5.0 

0. 


-0. 


-0. 
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APPENDIX  I 

COORDINATE  TRANSFORMATION  FROM  BODY  TO  GIMBAL  AXIS  SYSTEM 

I 


121 


2 


Figure  1-1.  Angies  Between  Gimbal  and  Body  Axes 


Transformation  for  Gimbal  Pitch  Angle  (#5) 
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w 
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Y 
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-Sin^j 

Cos^f  0 

Y  * 

B 
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Z  ' 

B 
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Axes 

rx  1 

G 

Cos^ 

Sin  Iffy  0 

Cos  Qy 

yg 

-Sin^ 
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0 

z 

G 

l 

0 

0  1 
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1 
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0 
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ion  for  both  Gimbal  Pitch  anqle  and  Gimbal 
or  Range  Determination  is  given  as: 


Sin-^  -Sin#Cosy? 

RXBA 

Cos  "\f/  -Sin^Sin$ 

RYBA 

0  Cos  6 

RZBA 
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appendix  II 

HIGH  FREQUENCY  ACTUATOR  PROGRAM  LISTING 


C**  |M(T|iai/4TtOM  *UOULE  FOt  T  I  &ft  (IMLIMCO  ACTUATOR 
Frequency  NOOCL****** 

C 

i'liitCUf  INC  C A  t 
CU**C*  CI4110I 
DlMt.\SIC*  I  Pll  100) 

0 1 n*  NS  I CN  V>  C  f  L  T  I  4  ) 

FCulVALFNCt  It!  I  1 1 6 ) • NMQEL f ) 

ELcH  VALcNCS  ICt  1160  )  .OUTPO  1 ,  (C(  1161  )  .deltqui 
?  JUlVALENCE ICl 116? I.OElTRu) 
tVul\Mt  if.LL  ICUSMI.H  ) 

LvUl  VALENCE  ICI25621.1PL  ) 

c 

NMyf l T*H 

d,  tl T I  l  I  •  Cl  110)1 
HOUM2I  -  C  (  1 1  0  7  ) 

30.-11  (si  ■  cunt) 

HliflMAI  •  Cl  1  LIS  > 

tfr/ilTIU  ■  ftOElT(l)  -  DEL  TPB  ♦  OfLTQB  -  OCLTRO 

QUUH2I  ■  8l)EL  T  I  2  )  -  OELTP&  ♦  0ELTG3  ♦  06 1  TAB 

JPHTMl  «  ftOCLril)  «  OELTPO  ♦  0ELTG3  -  OELT«B 

tuncrm)  «  oottr (4)  ♦  utLTPu  ♦  dfiTjh  ♦  uCltab 

ci  i ;-)3»  *  qcelti  i > 

C(ll(<n  -  6 CELT  I  2) 

Cl  11  m  *  ikUUH  H 
Cl  l  US)  -  b  DEL  T 1 4  I 
I  PL  I  N)  •  1100 
ll'UN»ll  •  1104 
|PL1N*2)  •  1103 
IPL  I  hi)  »  1112 
lPLlN*y)  -  1124 
I Pt I N*  S )  •  1126 
JPLIN*6)  •  1  1  >2 
!Pl(N*7)  -  1136 
l Pi  IN*  u«1140 
1Ti.(N»))«1144 
l?L  1  -it  10). I  148 
IPL  INM  l  hi  152 
N«\«12 
Cali  MksCT 
BclLR.S 
EKC 

C  ••  TICER  S  I mpl  l  P  l  E  D  ACTUATCR  MODEL 
C*M.««hlGK  FRECUESCY  MCCEL»«*»»« 

C 

Sl.dRCliTlNE  C4 
C 

CCPPCNCI4310I 

ct  MCNSICN  nc EL  TCI  4 ) ,00ELT14 ),OCELTCl 4), vARllCl) 

OI.-‘NS!Cf«  Q£ELTPt6)*£0L?CC(«) 

CIP'fNSlCN'lCfLTCIA) 

C I  PENS ICNPMPI4) , BCLTC31 4) 

C 

C • • 1 S  PU  T  CAT* 

EQUIVALENCE  €  C I  1163  I, Cl  ACT ) 

ecu  I  VALENCE  CCI  1309  ». Fm-H. ICC  1210) , f«H2 ) . 1 C 1 1  111 ) ,Fmm3 I » 1C  1 1112 1 fc 

«PP»-4  ) 

ECUVALENCE  ICI  112  l  ),  BTMAX  ) 

EQUIVALENCE  ICl 1 122  I , CACTO  ) 

ECuIVALCNCEICI UAOl.CELTPO )f I  Cl  1161 )f DEL T08) 

EQUIVALENCE  1CI  l  I62)»CELTPJ)  ) 

C  BCE  L  T  ARE  MEASCREC  IN  RAOjANS 
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c 

Cl'tKPufS  H*LP  CTM*  PQCLIES 

€CoW4irsCC  iCIllUl.SKUT) 

ECctvHffcCf  ICI  0)7  1  *01  EL  TC  ► 

CCljkVAifKCr  ICI  l  l  l  1  I'QSlHF  I  » 

ccuiv*itNCL  mmei.nuHfj) 

€uu1vtUfiC€  ICI  U!9)«0^U»FSI 
EtUlVAltNCC  ICI 1120)  ♦  BSuft Ml 

c 

C«*SUT(  VA*l*ecE  CUIPUTS 
eourpi  it  «  cm2Ti 
BCClfPI?)  -  Cl  im  I 

fcctimji  *  cimsi 

9CEITPIA)  *  Cl  11391 
BDklTl l l«C( 110)1 
BCCU  I2)»CI  no?) 

ecu  t 1 1 ) •  c i  nu ) 

QCClf <A  I  *  C  I  I  115) 

PT»6, 21)16/ 360 . 

FPF-  l  1 »  »f  *K| 

FfF  1  ))  9f  Ph  ) 

FPH I A )  •  F  NH  A 

c 

*i-i./ib. )/ieo./ieo. 

44*1* 

*)•  .h/180.*l./u.3 
A2-1./U0./  180.4.6/16.3/180- 
ei'Ctucr 
ei  *  tifA<«PT 

C** INPUTS  F*rp  FAIN  PROGRAM 

FCLlVMtKCe  K.I29<»5»»var  J 

Cl  1  ICO)*Ct  l  127) 

C  I  I  1 1 2  )  -  C  I  1  1  )9 1 
Cl  1106  I ’Cl  11)5) 

Cl  l  10 A l “ C I  11)1) 

C(ll?Vi-CIII55l 

cim?i*cui5ii 

Cl  1  12?)  *C 11147) 

Cl  l  1  2 A ) *C I  1143) 

C 

C • *F l fl  P  CfflcCT  ICN  BIAS 

BGElTCtl)  3  BOEltCII)  -  C6LTP0  ♦  DEL  TQB  -  DElTRO 

0CCITCI?)  *  CCCITCI2)  -  CtLTPB  ♦  OFiTC-8  ♦  OEITAB 

BCslTCl*)  =*  CCCLTCm  ♦  OCLTP3  ♦  DElTQO  -  DUTRO 

BCslTClAl  ■  BCcLTCIA)  ♦  r,£LrPB  ♦  OElTOO  ♦  ocltro 

c 

C«*ACTUATCR  f.VNAMlCS 
CC  20  1*1. A 

ccFiTcu>«i.*i:cLrcm‘pr 

*■  U»l  MA  *1 

ectTcci  i  >»cuia?»ki 
cc e l tc 1 1 1  *  ecstrpi  i » 

ceil  Oil  l*-t?/4l»0Cl.Tcn(  1  )-A3/Al»BOElT01  I  ) -A  A/A  1  •  0OEII  U  I 
•  ♦dCFLTCI  1  l/Al»Bl*FPH  I  )  /  A  i 

c 

C**SU«FAC6  FCSITICNllMIER 

if  ucsteccinm  -lt.  st  i  go  to  30 
eocuui  •  siGNiei  .oofuTinj 

J  •  Nf-LELT  <  1 

VARtJ!  -  COElTin 

If  1SICMI..  BCCCTPim  .N£.  5  IGF*  I  l  •  ^  OOEUTIlll)  GO  to  30 
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tCtlTCI  U  •  0. 

ICIlTFI I)  •  0. 

J  •  J  ♦  * 

•  DCeiTMl) 

30  CCNUNUC 

c 

esopf i-ecut  1 1  i/ft 

8Sl*F2*teUT  m/P  T 
Bbtar  )*i!CCi  r  13)/PT 
eso^Fs^cti r  iaj/pi 
ci  i  ion  •  at  Ft  in) 

CM  10/  I  -  KCfiUZ) 

CUIUI  •  UCELII  31 

c <  iii5)  •  cceiti*) 

C< l l*0»»CClTC3l  11 
CMU<.)*3tirC)l2l 

ciiifii «ccl mi  3 1 

C(  1152MCCLTCJI  A) 

C 

c •  •Cut put  ce-mvatives  of  state  variables  to  integration 
ci  l  loo  »  *  iK ft  r c 1 1 1 
Cl  1 10* I  •  BCEITCI  2  1 
Cl  l  lOtt  »  •  8CEITCI 3  l 
Cl  11 12  I  •  BCEITC14  I 
C.  1  12/MCM  100  ) 

ciiim«ciin*i 
ci ll)5t*Clll0ft> 

Cl  1131  )*CM104) 

RETURN 

EKC 
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APPENDIX  III 

HIGH  FREQUENCY  AUTOPILOT  PROGRAM  LISTING 
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C«*  T!Gf*  #uTO#!lOt  IMT|*Ui>MOX  NOOUlf 

MtCuf*CT  rocfl**«**« 

u:*cvum  cu 

CC**CN  14310» 

Olf’MKN  If  LI  1001 
equivalence  i c <  bjd.cas  ) 
ICllVAiC.NCi  tC  I  8*!I»IYS  I 
CCOtVAlENCC  ICI  BB3)»f/SS  > 
ECUlVAiEsCE  ICC  B<j?),fcYSS  ) 
EQUIVALENCE  I L  (  66*  >*CCA*YS  ) 
ECulVAitNCE  I C  I  6«!  ) 

ECuJvaienu  IC I  1%04 1 1  r^f  ) 
ECuIvalEsCE  (C<  Cft6BI,N*$UK  » 
ECoIvutf.CE  ICICBStl.Gv  ) 
EQUIVALENCE  IC(256U,N  I 

LwLUAlEnCE  <  C  c  2562),  I  Pi  ) 

c 

Kf-SUf*  •  h 
I  f  L  IN  I  •  6)0 
I  Pi  IN* 1  I  •  «06 
UilM2)  •  €08 
IKIN«)I  •  812 
IKIMO  •  6  16 
IPl  |N*3  >  *  u 20 
i pi iMt >  -  e 2 <• 

I f i (K« 7 )  ■  e?8 

I f L  Ift*  a  I  •  6  32 
IHINH)  *  336 
1  f  l  I  N  •  10  I  '  8Q0 
IHIMIII  •  880 
|fllN*l2»  -  866 

I  FL  IN* l  H  «Y12 
: Ft  IN* l*  )-*9CB 
t  FUN*  16  )  *890 
1  Pi  IN* 1  t ) *096 
I Fl IN* 18  I  *9C0 
!  ft  lN*l'3)«9C6 
N  ■  N *  20 

21  E  2  S  >  CCanvS 
E  v  S  •  CCanhS 

22  E2SS  >  E2S 
EtSS  •  EtS 
Cl  603 )  •  0. 

Cl  6071  •  0. 

Cl  Blli  •  0. 

Cl  8153  «  0. 

Cl  619)  •  0. 

Cl  823)  •  0- 

C I  827 )  •  0. 

Cl  Bill  •  0 . 

Cl  839)  •  0, 

»E  RUN 
ENC 

€••  IICfR  ACTCMliOT  MQCUCE 
C • • • • • «h  I  O  F  9  £  CU  t  NCY  N0CEL«»»»*» 

SUJ3CCT INE  Cl 
Cl>*CN  C  (  6 3  10 ) 

CIKENSICN  BCELTCl6),VAP| 10  l ) 

C 

C •  •  INPUT  CAT* 

ECU  IV At ENCE  ICI0B5C),hC  |M0  ) 
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1C*.  IV  Al  *  Kff  1C  lOMII.K  |4(  i 
CQUIVAL‘M.C  «C  <  e  5  7  > ,  u I  a  $  » 
ecu  hale -.a  i*:*  ii;ii«ri|M^*i 
EvUVAlcNCt  ICI  ) 

EQUIVALENCE  ICI  Uv».U*5».  1C  I  *16  l.uASOt  ,  1C  00  81  .OOELPOI 
f  QLlVSL  t*u  I  Cl  Q  15  » ,uOS  )  I  tl  <  M?  IfUdM;)  •  1C  <<M1I  BbCf  LfC) 
t  quIval  *n;e  to  •«>  n  ,«ucca  i,  in  as?  i  * eococ  i . ici*iQi)  ,evnca  i 
tQUlVAl  INC.  <C  I  ->C7  l  ,IVM  l  It  (CI4S0  I  »i  ictau  ICievi  I  »E 100C  I 
EUUlVAl  VOOI.eVMi*  >,  (C(O04),EVNUll 

EQUIVALENCE  ICIC3*‘I»C/  > 

C  Cl  OS?)  1HU  Cl  »6H|  A«E  USCO  QV  ECHUUM 

EQUIVALENCE  4 C  <  0*43)  ,  M.\jI  ) 

EClIvalCNCL  1CIC864), Tauv  » 

EQUlVAuNCE  1CI  H65|,T(Y1  ) 

ECLlVALCr.Cc  1 C  «  8*CI,TCY?  ) 

tQUlVALLNCE  <  C I  07M,MU  ) 

c 

c«*iNF«jrs  f?Cf  ctler  pcclles 

EQUIVALENCE  1C  1 C  35? I , QPh l  I 

EQUIVALENCE  1 C  I  0355  )  ,  B  PM  D  » 

ECt.WAlt*Ct  1C(C*.0!),EC  I 

ECLtVALff.Cc  ICU'.OI.ET  I 

EQUIVALENCE  (Cl COAfl I.NESVN  i 

eclival c *. c £  icun9).wp  i 

EQUIVALENCE  *  C  1  1  KO  I  .  «  Q  l  I 

ECLIVAL  ESce  « C  1  174*1. V.C  I 

EQUIVALENCE  1  C  t  l  7 «.  N  )  ,  WS  C  I 

EQUIVALENCE  ICll74?l,WR  I 

EQUIVALENCE  I C  I  l  751 1 *CR AC  I 

C 

t • ■ I N  Pul  5  FUC?*  KAIN  P&CCRAM 

*  n.  f  y  m_  fsrt  K  *  7C’}r  )  t  r  ) 

EQUIVALENCE  1C i 2965  I •  VAR  ) 

EQUIVALENCE  (CI2664),tfR  I 

C 

C««  STATE  VAR  l  A  9  L  E  OUTPUTS 

EQUIVALENCE  1C(  oori.e**HSt) 

EQUIVALENCE  (Cl  QOJI.Ci'MS  I 
ECLlVALtsCE  I c  1  B04),*CSCu  > 

EQUIVALENCE  (Cl  8071, wCSP  > 

EQUIVALENCE  I C  (  bUfc  )  »  W.SC  ) 

EQUIVALENCE  IC(  Hlll.WQS  ) 

EQUIVALENCE  (Cl  9l2),w*StO  I 
EQLIVAl  ENul  I C  I  tilSI.WTSP  I 
EQUIVALENCE  ( C I  8UI.WRSC  ) 

EQUIVALENCE  (C(  9151,  VMS  1 
EQUIVALENCE  1CI  820 > , CSuKJC > 

EQUIVALENCE  4  C I  *J23l,tSUMO  I 
EQUIVALENCE  4  C  «  624  ) , CSU*EC  ) 

EQUIVALENCE  (C(  62  7  1  ,  E  5U.ut  I 
EQUIVALENCE  (Cl  9?6)»t/SCC  I 
EQu1vtLfr.CE  (C(  6311, E/SP  I 
L^ulvALfNCE  €  C I  8321. E/SC  I 
EQUIVALENCE  (C(  8351. E/S  » 

EQUIVALENCE  tC(  03/  l.frSCO  I 
EQUIVALENCE  (C(  UTOI.CYSP  ) 

equivalence  (ci  ascifEYsu  > 

EQUIVALENCE  I C  4  052), tYS  I 
E  Cu I  V  A  L tNCE  ( C  <  ft'*‘  E/SSC  I 
EQUIVALENCE  1CI  b  E/SS  > 

equivalence  ici  aasi.EYssc  ) 

EQUIVALENCE  ICI  ean.EYSS  I 
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c 

C**CuTPuTS 

tctl VAtfNCE 

c 

c«*cther  cl t puts 

£Co  I  v* LEKCE 
cvl I  VALENCE 
EUlVAlENtE 
ELt 1  V  AL  EKCC 

c 


ICt  fl  5  T  )  • BCElTC  I 


1C  4 087t ),E*ftR  I 
101 0U6S ) . EY*R  I 
ICI  «iS).wCC  I 
ICI  87CI.WAC  ) 


C«*GLICakCE  SIGNAL  $HAP|*C 
€  lCCa»  SCCCC 
Evm:m«Evkco 


€2sc  •  ezsp 

tvsc  *  EY5P 

fc/SCC  •  Tau2 #M AUZ •  ! Cl ■  £*  -  EZS) 
fvscr  -  tauV'I tauy*igy*€y  -  £▼$) 
f/SSC  *  TauE«I CiSC/TAUL  ♦  6 zs  • 
EySSC  •  tAUY.IE  VSC/TAUL  *  6  Y  s  - 

c 

C**G«AVi  TV  A  K  C  RATE  9  I  AS 
*;c  •  f  MS  4  CO  I  AS 
W*C  ■  kYSS  ♦.  «Q|AS 
C 
C 

C««OCCY  P«TE  SrAfMHG  ANL  CYRQ  DYNAMICS 
M.SC  •  wCSP 
wasc  *  nrfSP 


2  .  •  E  2  SO  I 
2 . *E YSO I 
EMSI 
EYSSI 


wocn  •  <J4.2M0A.2MWQ  -  MI'S  I  -  2.AW050) 
hASCC  ■  <J4. 2 M 94. 2M WR  -  WAS)  -  2.»WR50) 
If  CA8SUCSI  « l  E  •  ?0.)  CO  10  3C 
*r.S  *  SICMJC.,  V.CS  I 


VA*  INMS'jM*  J>  -  WGS 

If  <StGMl.»  WQSPI  .*E.  SICNIU,  wQSK  CO  TO  3C 
US*  *  C. 

*CSC  •  0, 

vtfi  (NhS0«*2  J  •  WCSP 
30  If  laesiwpsi  .IE.  ?0.I  CO  TO  32 
-«S  -  SICM  JO..  HP S 1 
VAA  IM-$LP*S  »  •  V.  R  S 

If  I S  I  CM  |  • .  WRSPI  .ME  .  SICNIU.  WASH  CO  TO  32 

MS*  ■  0. 

WRSC  »  0. 

VA«  (NfSLR***  1  ■  HR$P 
32  CCMIKUE 

C 

If  t  AtfSftCSI  .CT  .30.  )UCS-wQSi;/RO.*WOS 
If IA9SILPSI .CT ,30. IL«S«wR5C/90.4WRS 
LLSC-I (t»CSD/80.*WQ5J-UCS)-8C0. 

0*SCM  I »K  SG/00 .  **«*  S  I  -OR  S  )  •  a  CO  - 

t 

C«*SUP^ATICN  CF  RATE  Camping  AND  GUIDANCE  SIGNALS  AND  THE  I A  OcMVATlVES 
r'2RR  •  UCS  -  wgC 

EYAR  *  OR S  *  WRC 

C 

OKR  •  *95 

If  IT  .IT.  TCY  ?  I  UXR  •  A. 25 
If  IT  .If.  TCVi  J  UK«  •  0. 

C 

ISL*CC  <  UKAME2RR  -  EYfifl) 

CStMC  4  U4RME2RR  ♦  £YRR| 

IFUasa5UMEC).GT.9  0.  IESUME0»S1CN(  60..ESUME0) 
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l*  ues i«iu*cc  j  .cr.cc.  >ESu**cc-siLNi4a.,ESL>«CD) 

c 

C»«TCT/l  r.U«CANCe  SIGNAL  $»-APlSG  «‘iO  LIMITING 
tltiC  »lCO,*H  £M/“CC/9  .  ♦*  SL*0  I  -  E  10'.  rt  ) 

£v*C 1  *  I COO -  * ( I tSyxfcC/d  .  *t it«c i-£VND* I 
!F  UBS  IEVMJ1  I.GT.ICCQ.  )£VNl«»C'GNF  *cSum£0/8. 

I  -  (  N  (  l  •  ,  t  V  N  C  1  ) 

if u j* ceucc  i.ot.icoo.  )Eucn*esuHO*esuMtjo/o. 
i-sicm  i.,eiccc » 

If  I  /f»s  (  6  VAC  1  )  -GT.IOCO .  )  EVNC  l«s  ICNf  ICOO.tCVNDl  I 
I  C  { -s  OS  I  6  XCCC  I.GT.1CT0.  )  E  1 C  C  C  r  S  I  CM  ICCO..EIOOC) 

EVACC*  SVhCR 

trccc- 1  icc« 

if  ubsucucai  ter.  hlihqj  ccocr  •  signihlimo,  eodcai 
if  uesitvsc.n  .ot.  hlihei  evncr  ■  signihlihe.  evncai 
c 

C ■ • RC l l  SICNAl  SHAPING 
U*P  »  .33 

IF  (  T  .LT,  TCY2  I  UK  P  -  1.65 
If  IT  .LT.  f  CYI  )  L'<P  «  0. 

UFHS  t'KPa  I  OPHC/I?.  *  6  PH  t  I 
PPMSC06.*  (UPHlS-'.-PHlS  ) 

RLLlPf  *-IL’  Pf  ISG/1QP.  ♦  0  P  ►  I  5  ) 

CCElPL-lCCO.al  bULPT-BCCLPC  ) 

rc*er m x/ia. 

If  I  Ah*  IPPHI  3  )  .GT.CC  mPh  IS-SIGNt  BO.  APHIS) 

IF  U£i  C  eCtl  FT  )  •  CT.SC'*.  tOCE  I.PC  •CCLL  P  f-S  IGNI  .5.OUCLP0 I 
IFUlTilECELPC  J  .CT.5CC.  IbCtlPD- 5  IGN(  5C0  .  .  OOtL  PU  I 
IH  ACSICCtLPCl  .f F.CC IBCcLPC •SIGNl3n,30ELPC I 
!»  I  I  I  .  C  I  .  .  9  )  .  AMI .  (  T  .1  f  .  .  )  I  I 

1-MKI6.6II  I ,  C.  I  I  I,  Cl  IH  I.  Li  I  *2  su  r  Ol  (COM  htthSi  fCI  I *6 U l -l .3* 
* iu . 7  * 

6  FCSPAMlH/t  15. 1P7E1S.YM 
C 

C**AUTCPllCT  CUTPUT  CL»A£N?S  TO  £ACh  ACTUATUR  IFROM  SUMMATION  AAPSI 
BCflTCl  l)  -  EUCCR  -  BCELPC 

8C£lTCI2)  ■  £VnC*  -  BCELPC 

CCtlTCO)  ■  £Gt  C«  ♦  BCELPC 

OCEUCU)  •  CVNCR  ♦  BCELPC 

fltKRN 
EM) 


132 


Y7^ 


J 


>.;*  l?5i  /3y 


UNCLASSIFIED 


Sccunl^latiuficiilio* 


i  DOCUMENT  CONTROL  DATA  RAD  1 

1  (Sacwifr  rUiidkftltaM  of  tflM.  Mf  aAtnaci  amt  tmdmwk 14  amlarion  MU  Me  wntwrwd  ■A—  Me  omclf  re pmt  is  rlaeitfMj  | 

Department  of  Mechanical  Engineering 
University  of  Florida 

Gainesville,  Florida 

_ 1 

1  *CPO*T  TITC* 

CLOSE  AIR  SUPPORT  MISSILE  GUIDANCE  AND  CONTROL  STUDY 

Volume  I.  Six-Degree-of-Freedom  Simulation 

*  OUCNieiivt  motes  (Typ*  *> inmmt  *md  IncfMlM  diM«) 

Final  Report  (9  December  1970  to  9  December  1971) 

%  Au  »  HORItl  (Finl  naw,  Mdtf,  Bmmi  MM) 

J .  Mahig 

«  MCPOM  T  OA  TC 

December  1971 

T.  TOTAL  NO-  OF  P*6EJ  7b.  NO  OF  REFT 

142  None 

M  C  0*1  TNAC  T  ON  «R*NT  MO 

FO8635-71-C-0073 

670b 

Task  No.  01 

«.  Work  Unit  No.  010 

M.  ONICIM  ATOM'S  REPORT  MUMBF  NISI 

•6.  OTHER  REPORT  MOIST  (Aftf  r-Ttiet  11 1  *1 - Ruf  M  Mliflfd 

fftJl  T9por1) 

AFATL-TR-71-169 

19  DtlTMISUTION  STATCMEM'  "  W 

Distribution  limited  to  U.  S.  Govexro&ent  agencies  only; 

a&tdf*  distribution  limitation  applied  December  1971.  Other  re¬ 
quests  for  this  document  must  be  referred  to  the  Air  Force  Arma¬ 
ment  Laboratory  (DLWG) ,  Eglin  Air  Force  Base,  Florida  32542. 

i«  !uppum{hi«rt  moth 

Available  in  DDC 

12-  SRONIORlNC  MILITARY  ACTIVITY 

Air  Force  Armament  Laboratory 

Air  Force  Systems  Command 

Eglin  Air  Force  Base,  Florida  32541 

This  report  describes  in  detail  a  six-degree-of-f reedom  pro¬ 
gram  which  can  be  used  to  determine  the  trajectory  and  miss  dis¬ 
tance  of  a  missile  system.  The  options  for  the  program  are  such 
as  to  permit  variation  of  the  aerodynamics,  seeker,  autopilot 
actuator  and  missile  motor  performance  for  the  purpose  of  accu¬ 
rately  simulating  a  given  missile  design  and  evaluating  the  ef¬ 
fects  of  changes  in  system  parameters.  Sufficient  detail  has 
been  included  in  the  text  in  order  to  minimize  the  users'  effort 
needed  to  know  how  to  update  or  modify  the  program  for  his  pur¬ 
poses  . 


DD 


FORM 

I  MO*  •  • 


1473 


UNCLASSIFIED 


Security  Cltttiftrttion 


